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READER’S  REFERENCE 

This  report  is  Volume  IV  of  IDA  Study  S-341  "Determination  of 
Winds  and  Other  Atmospheric  Parameters  by  Satellite  Techniques." 

The  volumes  describing  the  study  are  the  following: 

Volume  I  Summary 

II  Potential  Needs  for  Determination  of 
Atmospheric  Parameters 
III  Techniques  for  Determining  Winds  and 
Other  Atmospheric  Parameters 
IV  Physics  of  the  Atmosphere 
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PREFACE 


The  work  reported  in  this  volume  was  accomplished  under  IDA’s 
Independent  Research  Program.  It  is  included  here  as  part  of  S-341 
because  of  its  relevance  as  a  background  to  tne  measurements  of 
winds  and  densities  from  satellites  and  to  the  computational  simula¬ 
tion  of  winds. 
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ABSTRACT 

The  physical  characteristics  of  the  terrestrial  atmosphere 
are  reviewed  for  the  purpose  of  providing  a  physical  background 
for  a  study  of  requirements  and  techniques  for  the  determination 
from  satellites  of  winds  and  other  atmospheric  parameters.  The 
nature  of  the  atmosphere,  its  physical  processes,  and  the  vari¬ 
ability  of  its  parameters  are  described.  Topics  treated,  in  turn, 
are  solar  radiation  as  its  affects  the  atmosphere,  the  troposphere, 
stratosphere  and  mesosphere,  the  thermosphere  and  the  ionosphere. 
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INTRODUCTION 


Physics  of  the  Atmosphere,  volume  IV  of  a  study  of  requirements 
and  techniques  for  the  determination  of  winds  and  other  atmospheric 
parameters  by  measurements  made  from  earth  orbiting  satellites,  is 
a  review  of  the  physical  characteristics  of  the  terrestrial  atmosphere 
which  is  to  be  observed.  The  nature  of  that  atmosphere,  its  physical 
processes  and  the  variability  of  the  parameters  by  which  they  may  be 
described  are  essential  knowledge  required  for  the  understanding  and 
evaluation  of  measurement  data,  for  the  design  of  instrumentation 
systems,  and  for  the  computational  simulation  of  the  atmosphere. 

Much  of  the  material  contained  in  this  description  may  seem  either 
unnecessarily  complicated  for  the  needs  of  the  pragmatic  reader,  or 
roo  elementarily  tutorial  for  the  knowledgeful  specialist.  The 
meteorologist,  for  ii. stance,  will  find  little  new  in  the  discussion 
of  the  troposphere.  The  radio  communications  physicist  will  not  be 
particularly  enlightened  by  the  discussion  of  the  ionosphere.  How¬ 
ever,  in  this  volume  is  the  background  essential  for  a  proper  appre¬ 
ciation  of  the  requirements  and  of  the  measurement  techniques  discussed 
in  Volumes  II  and  III-  It  is  hoped  that  even  the  specialists  will 
find  some  value  in  this  review  of  the  total  picture. 

In  the  following  pages,  the  influence  of  the  solar  radiation 
on  the  terrestrial  atmosphere  is  described  first.  The  sun  provides 
the  driving  energy  of  the  atmosphere  and  causes  its  large  variability. 
Most  of  the  sun’s  energy  arriving  at  the  earth  is  absorbed  by  the 
solid  earth,  the  oceans  and  by  the  lower  atmosphere,  the  troposphere. 

A  very  small  fraction  of  this  energy  is  deposited  in  the  upper  atmos¬ 
phere;  however,  its  effect  is  very  large. 


With  the  energy  derived  from  insolation  as  the  primary  cause  of 
variability,  the  terrestrial  atmosphere  varies  most  strongly  as  a 
function  of  altitude.  In  trying  to  characterize  the  variation  with 
altitude,  a  number  of  models  of  the  atmosphere  have  been  described 
as  "standards"  for  use  in  detailed  studies.  Of  these,  a  recent 
example  commonly  cited  is  the  U.  S.  Standard  Atmosphere  of  1962. 
Shown  in  Fig.  1  is  the  variation  of  temperature,  pressure,  density, 
and  molecular  weight,  with  altitude,  according  to  the  U.  S.  Standard 
Atmosphere,  1962.  Also,  in  Fig.  2  is  the  variation  with  altitude  of 
other  atmospheric  parameters:  number  density,  collision  frequency, 
mean  free  path  and  particle  speed,  according  to  the  same  standard. 


TEMPERATURE  (°K) 


PRESSURE  (mb)  OR  DENSITY  (g  m-3) 


FIGURE  1.  Temperature,  Pressure,  Density,  Molecular  Weight  (U.S. 
Standard  Atmosphere,  1962)(  After  Valley,  1965) 
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FIGURE  4 


FIGURE  3.  Moisture  in  Lower  Troposphere  (After  Valley,  1965) 
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1966) 


of  excursions  of  two  of  the  parameters :  mass  density  and  tempera¬ 
ture.  Mass  density  may  vary  from  plus  60  to  minus  40  percent  of 
standard  magnitude  at  altitudes  below  120  km,  and  by  a  factor  of  100 
at  altitudes  of  600  km.  Temperature  may  vary  over  a  range  of  10  to 
25  percent  of  standard  value  at  altitudes  below  120  km,  and  by  a 
factor  of  three  at  an  altitude  of  600  km. 

The  variation  of  the  atmospheric  temperature  profile  has  been 
the  basis  for  the  nomenclature  of  the  atmospheric  regions,  which  is 
currently  common  and  which  will  be  used  in  this  report  hereafter  and 
is  given  in  Fig.  5.  The  extrema  A  the  temperature  profile  are 
called  the  tropopause,  the  stratopause  and  the  mesopause,  and  occur 
at  altitudes  of  about  20,  50,  and  80  km,  respectively.  Below  the 
tropopause  is  the  troposphere ;  between  tropopause  and  stratopause 
is  the  stratosphere;  between  the  stratopause  and  the  mesopause  is  the 
mesosphere,  and  above  the  mesopause  is  the  thermosphere.  The  physics 
of  each  of  the  regions  is  distinctive. 

The  review  of  the  physics  of  the  troposphere  is  concerned  mainly 
with  estimating  the  characteristic  dimensions  of  the  large-scale 
phenomena  which  are  of  importance  to  numerical  weather  prediction,  and 
with  a  treatment  of  the  theory  and  examples  of  the  atmospheric  wave 
motion  largely  a  paraphrase  of  Eliassen  and  Kleinschmidt  (1957), 

Mintz  (1961),  Lorenz  (1965),  and  Starr  (1968).  Such  wave  motions, 
including  both  impulsive  and  quasi -continuous  phenomena,  have  a  highly 
complicated  theory,  which  is  however  a  powerful  tool  for  understanding 
the  dynamic  processes  of  the  atmosphere. 

In  the  review  of  the  physics  of  the  stratosphere  and  the  meso¬ 
sphere,  some  distinctive  characteristics  of  this  stable  altitude  re¬ 
gion  are  given.  The  region  is  quite  unlike  either  the  thermosphere 
above  it,  or  the  troposphere  below.  To  a  degree,  it  behaves  as  an 
energy  trap,  it  inhibits  the  transmission  of  several  forms  of  energy 
flux  from  outer  space  to  the  earth,  and  from  the  earth  outward.  Its 
importance  lies  in  the  fact  that  it  uossibly  serves  as  an  upper  bound 
to  the  troposphere  (and  thereby  affects  our  weather)  and  certainly 
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of  excursions  of  two  of  the  parameters:  mass  density  and  tempera¬ 
ture.  Mass  density  may  vary  from  plus  60  to  minus  40  percent  of 
standard  magnitude  at  altitudes  below  120  km,  and  by  a  factor  of  100 
at  altitudes  of  600  km.  Temperature  may  vary  over  a  range  of  10  to 
25  percent  of  standard  value  at  altitudes  below  120  km,  arid  by  a 
factor  of  three  at  an  altitude  of  600  km. 

The  variation  of  the  atmospheric  temperature  profile  has  been 
the  basis  for  the  nomenclature  of  the  atmospheric  regions ,  which  is 
currently  common  and  which  will  be  used  in  this  report  hereafter  and 
is  given  in  Fig.  5.  The  extrema  A  the  temperature  profile  are 
called  the  tropopause,  the  stratopause  and  the  mesopause,  and  occur 
at  altitudes  of  about  20,  50,  and  80  km,  respectively.  Below  the 
tropopause  is  the  troposphere;  between  tropopause  and  stratopause 
is  the  stratosphere;  between  the  stratopause  and  the  mesopause  is  the 
mesosphere,  and  above  the  mesopause  is  the  thermosphere.'  The  physics 
of  each  of  the  regions  is  distinctive. 

The  review  of  the  physics  of  the  troposphere  is  concerned  mainly 
with  estimating  the  characteristic  dimensions  of  the  large-scale 
phenomena  which  are  of  importance  to  numerical  weather  prediction,  and 
with  a  treatment  of  the  theory  and  examples  of  the  atmospheric  wave 
motion  largely  a  paraphrase  of  Eliassen  and  Kleinschmidt  (1957), 

Mintz  (1961),  Lorenz  (1965),  and  Starr  (1968).  Such  wave  motions, 
including  both  impulsive  and  quasi-coritinuous  phenomena,  have  a  highly 
complicated  theory,  which  is  however  a  powerful  tool  for  understanding 
the  dynamic  processes  of  the  atmosphere. 

In  the  review  of  the  physics  of  the  stratosphere  and  the  meso¬ 
sphere,  some  distinctive  characteristics  of  this  stable  altitude  re¬ 
gion  are  given.  The  region  is  quite  unlike  either  the  thermosphere 
above  it,  or  the  troposphere  below.  To  a  degree,  it  behaves  as  an 
energy  trap,  it  inhibits  the  transmission  of  several  forms  of  energy 
flux  from  outer  space  to  the  earth,  and  from  the  earth  outward.  Its 
importance  lies  in  the  fact  that  it  nossibly  serves  as  an  upper  bound 
to  the  troposphere  (and  thereby  affects  our  weather)  and  certainly 
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is  involved  in  the  full  exploitation  of  radio  communications  and  in 
the  military  problem  of  targeting  of  weapons  from  higher  regions  in 
the  atmosphere  and  outer  space. 


FIGURE  5.  Atmospheric  Temperature  Profile  and  Different  Regions 
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In  the  review  of  the  physics  of  the  thermosphere  and  ionosphere 
are  given  some  features  of  a  quite  distinctively  different  region 
which  affects  radio  communications,  targeting  of  ballistic  weapons 
and  the  predictability  of  low  orbits  of  earth  satellites . 
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I.  INFLUENCE  OF  SOLAR  RADIATION  IN  THE  ATMOSPHERE 


The  driving  source  in  producing  the  variability  of  the  terres¬ 
trial  atmosphere  is  the  sun  which  is,  itself,  a  very  complex  system. 
The  regions  of  the  visible  sun  have  been  named  in  a  manner  illustrated 
in  Fig.  6. 


FIGURE  6.  Regions  of  the  Sun  ( After  Friedman,  1965) 

mchuhI  PAK  SUNK 


The  interior  of  the  sun  consists  mainly  of  hydrogen  which  by  exo¬ 
thermic  nuclear  transformation  produces  heavier  elements.  The  heat  of 
the  exothermic  process  is  transferred  from  the  outer  shell  by  a  com¬ 
bination  of  convection  and  radiation  processes  and  from  the  inner  core 
by  radiation  alone.  Only  a  small  fraction  of  the  exothermal  energy 
escapes.  It  is  thought  that  the  sun  has  been  in  essentially  its  pres¬ 
ent  form  for  about  seven  billion  years,  and  will  remain  essentially 
in  its  present  form  for  another  three  billion  years  ;  then  it  will  be¬ 
come  a  "red  giant."  As  a  "red  giant,"  it  will  expand  to  a  diameter 
which  will  extend  to  the  present  orbit  of  Mars.  It  will  consume  the 
interior  planets  (including  the  earth). 

Be  that  as  it  may,  the  present  sun  is  for  the  time  of  recorded 
observation  in  an  essentially  stationary,  although  complex,  state. 

The  most  obvious  visual  feature  of  the  sun  is  the  photosphere,  which 
we  may  see  with  the  naked  eye.  It  subtends  about  a  half-degree.  Its 
minimum  surface  temperature  is  about  5700  to  400D°C.  The  sun  suffers 
small  temporal  variations  of  apparent  size  and  temperature  that  are 
due  principally  to  the  variations  of  orbital  geometry-  and  to  the  ab¬ 
sorption  in  so  called  sunspots  by  cooler  gases  of  a  small  fraction  of 
its  total  radiation.  The  chromosphere  and  corona,  on  the  other  hand, 
are  at  much  higher  temperatures  and  are  much  more  variable. 

The  chromosphere  extends  from  the  surface  of  the  photosphere  out 
about  10,000  km  and  is  extremely  variable  in  thickness.  It  is  made 
up  of  cool  spicules  which  penetrate  into  the  corona  and  of  hotter 
inter-spicules.  Its  temperature  and  density,  illustrated  in  Fig.  7, 
range  from  about  440C  °K  and  (10)19  cm  9  near  the  photosphere  to 
about  40,000  °K  and  (10)9  cm-3  at  a  distance  10,000  km  from  the 
photosphere . 

Shown  in  Fig.  8  is  a  model  of  chromospheric  and  coronal  temper¬ 
ature  and  density.  One  may  see  that  at  about  10,000  km  from  the 
photosphere,  there  is  a  sudden  transition  to  the  much  higher  temper¬ 
atures  of  the  near  corona  (millions  of  degrees  Kelvin).  The  density 
distribution  and  the  temperature  are  strongly  variable  according  to 
magnetic  fields  in  the  sun. 
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FIGURE  7.  Density  and  Temperature  in  Chromosphere 


H.  W.  Babcock  (1961)  has  given  the  most  plausible  explanation 
of  the  principal  mechanism  that  produces  the  variability  of  the 
visible  sun.  His  ideas  are  illustrated  in  Fig.  9.  The  magnetic 
field  of  the  sun,  thought  to  have  an  exterior  expression  like  that  of 
a  in  Stage  1  of  Fig.  9  is  dragged  at  different  speeds  in  the  interior 
by  the  solar  plasma.  The  solar  plasma,  located  within  the  photosphere, 
is  known  to  rotate  more  rapidly  at  the  solar  equator  than  at  the  solar 
poles.  The  distribution  of  the  angular  velocity  (w)  is  as  shown  in 
b  in  Stage  1  of  Fig.  9.  Moreover,  the  magnetic  field  lines  do  not 
penetrate  to  a  great  depth  within  the  sun,  but  rather  as  illustrated 
by  H  in  b.  As  a  consequence  of  the  shallow  penetration,  and  of  the 
dragging,  the  magnetic  field  lines  are  wound  around  the  solar  sphere, 
as  illustrated  in  Stage  2  of  Fig.  9.  With  more  and  more  lines  pressed 
ever  more  tightly  together  with  successive  rotations  of  the  solar 
equator  with  respect  to  the  solar  interior,  the  consequent  lateral 
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FIGURE  8.  Model  of  Chromospheric  and  Coronal  Temperature  and  Density 
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pressure  may  cause  one  of  the  lines  to  be  extruded  as  a  small  loop, 
as  illustrated  in  Stage  3  of  Fig.  9.  The  small  loop  of  magnetic 
field  is  the  site  of  the  sunspot  as  cooler  gases  from  below  are 
entrained  in  the  magnetic  field  to  flow  along  the  loop  into  a  higher 
region  of  hot  gas .  The  sunspot  itself  is  the  region  of  the  sun  in 
which  radiation  normally  visible  from  the  hot  gas  is  self-absorbed 
by  the  cooler  entrained  gas.  Seen  from  the  side,  the  hot  gases  in  the 
magnetic  loops  are  the  "solar  prominences"  which  project  out  in  the 
corona  where  the  temperature  is  the  maximum. 


FIGURE  9.  Development  of  Bipolar  Field  (After  H.W.  Babcock,  1961) 
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We  now  begin  to  have  a  basis  for  understanding  some  reasons  for 
the  variability  in  character  of  the  solar  sources  of  radiation  that 
affect  the  terrestrial  atmosphere.  Some  cf  the  ideas  are  outlined 
in  Fig.  10. 
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FIGURE  10.  Source  Location  of  Solar  Emission 


The  atomic  processes  that  give  rise  to  the  radiation  may  be 
classified  as  transitions  from  bound  to  bound  states,  from  bound 
to  free  (and  vice  versa)  states,  and  from  free  to  free  states.  The 
energy  relations  characteristic  of  the  three  types  of  transition  are 
listed  in  Fig.  10  as  functions  of  density,  temperature,  and  frequency 
according  to  Aller’s  classification  (1963).  By  using  the  energy  rela¬ 
tions,  and  the  temperature  and  density  distributions  of  Fig.  8,  one 
perceives  that  the  visible  spectrum  originates  within  the  cool  photo¬ 
sphere,  that  the  Hydrogen  Lyman-alpha  line  originates  within  an  area 
of  about  a  1000  km  above  the  photosphere,  that  the  Schumann-Range 
continuum  originates  in  the  mid-chromosphere  about  6000  km  above  the 
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photosphere,  and  that  the  X-radiation  of  wavelength  greater  than  a  few 
Angstroms  originates  in  the  near  (hot)  corona.  In  general,  the  radio 
frequency  radiation  must  originate  in  the  far  corona  where  the  smaller 
densities  make  possible  the  radiation  to  earth  without  self -absorption. 
Although  the  mechanism  is  not  clearly  understood,  particles  radiated 
to  the  earth  from  "solar  sources"  do  not  originate  within  the  sun,  but 
also  come  from  a  location  in  the  far  corona.  Apparently  they,  like 
the  X-radiation,  arise  from  the  interaction  of  jets  of  magnetically 
confined  cooler  gases  from  the  sun  into  the  hot  solar  corona.  Travel¬ 
ling  at  slower  speeds  than  the  electromagnetic  radiation,  the  par¬ 
ticles  travel  from  the  source  to  the  earth  by  a  different  route,  and 
once  within  the  magnetosphere  are  further  constrained  to  travel  along 
the  lines  of  the  earth's  magnetic  field. 

The  spectrum  of  the  solar  electromagnetic  radiation  in  the  visible 
range  may  be  roughly  characterized  as  a  group  of  lines  and  continue  in 
an  envelope  like  the  black -body  characteristic  of  6000  °K.  For  shorter 
wavelengths  (less  than  2000  A)  the  radiative  spectrum  is  as  shown  in 
Fig.  11.  For  wavelengths  1000  to  2000  A,  the  radiation  appears  in 
general  to  emanate  from  a  black-body  of  temperature  4500  to  5000  °K. 
Notably  differing  from  this  generalization  are  lines  of  helium  and 
hydrogen.  Since  these  gases,  however,  are  the  principal  constituents 
of  interplanetary  space,  the  lines  represent  an  accumulation,  in  de¬ 
graded  form,  of  energy  originating  in  much  more  energetic  form,  but 
not  reduced  to  the  lowest  energy  levels  (HLyA,  He  I,  He  II)  charac¬ 
teristic  of  the  two  species. 

In  contrast  to  the  radiation  at  wavelengths  greater  than  1000  A, 
there  is  the  radiation  at  wavelengths  less  than  200  A.  This  appears 
(as  we  have  already  shown  to  be  plausible)  to  originate  from  different 
sources;  a  black-body  of  highly  variable  area  (peak  energy  varying  by 
factor  ten  between  minimum  and  maximum  of  eleven-year  solar  sunspot 
cycle)  and  of  temperature  about  5 JO, 000  °K  seems  to  be  a  source  of 
radiation  greater  than  20  A,  and  a  different  black -body  of  even  more 
variable  area  (the  shaded  zones  in  Fig.  11  which  vary  according  to 
solar  flare  activity)  and  of  tempt' ature  of  about  two  billion  °K.  The 
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FIGURE  1 1 .  Solar  Spectral  Energy  Below  2000  A 
(After  Friedman,  1962) 


dotted  lines  of  Fig.  11  show  an  energy  flux  curve  expected  of  an  area 
like  that  of  the  photosphere  at  various  black -body  temperatures.  One 
may  see  by  comparing  the  actual  flux  curves  with  the  hypothetical 
black -body  curves  (dashed  lines),  that  the  solar  radiation  at  wave¬ 
lengths  greater  than  1000  °K  seems  to  emanate  from  a  source  area  ap¬ 
proximately  that  of  the  photosphere,  but  that  solar  radiation  at  wave¬ 
lengths  less  than  200  A  emanates  from  a  source  area  that  is  orders  of 
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magnitude  smaller  than  the  photosphere.  Photographs  of  the  sun,  made 
with  X-radiation,  confirm  that  the  solar  sources  are  indeed  spots  on 
the  sun  in  locations  roughly  coinciding  with  solar  sunspot  activity 
and  of  area  of  the  order  of  one  millionth  of  the  photosphere  area. 

Seen  in  profile,  the  activity  which  causes  the  X-radiation  is  at  the 
extremes  of  the  solar  prominences  where  the  chromosphere  penetrates 
the  corona. 

The  measured  variability  of  the  solar  X-ray  flux  is  shown  in  Fig. 

12  in  which  is  compiled  measurements  by  a  number  of  observers.  Radia¬ 
tion  in  the  band  44-60  A  varies  over  a  factor  of  about  ten  between 
solar  minimum  and  the  maximum.  Radiation  in  the  band  2-8  A,  on  the 
other  hand,  is  seen  to  vary  over  five  orders  of  magnitude  between 
solar  minimum  quiet  conditions  and  solar  maximum  (2  +  flare)  conditions. 


FIGURE  12.  X-Ray  Solar  Flux  Measurements 
(After  Kreplin,  1963) 


The  occurrence  of  radio  bursts  and  flares  has  been  described  by 
H.  W.  Dodson'  (1959),  as  shown  in  Fig.  13.  Also  shown  in  the  figure 
are  the  minimum  distances  from  the  photosphere  of  the  sources  of  the 
radiation  perceived  at  the  earth,  as  indicated  by  the  limiting  plasma 
frequency 
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where  nQ  is  the  electron  number  density  (it  is  assumed  to  be  equiva¬ 
lent  to  the  number  density  of  the  model  of  Fig.  8).  Of  particular 
note  in  Fig.  13  is  the  3000  MHz  radiation  which  roughly  correlates 
in  time  with  solar  ultraviolet  radiation  and  is  indicated  as  possibly 
originating  similarly  in  the  chromosphere  about  2000  km  from  the 
photosphere.  Radiation  at  3000  MHz  is  taken  as  an  index  of  solar 
activity  as  a  parameter  differentiating  between  CIRA  1S65  (COSPAR 
International  Reference  Atmosphere  of  1965 )  models  of  the  thermo¬ 
sphere.  Such  an  index,  while  a  feeble  measure  by  which  to  character¬ 
ize  variability  of  the  thermosphere,  is  better  than  none  at  all. 
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FIGURE  13.  Radio  Bursts  and  Flares 

(Adapted  from  H.  W.  Dodson,  1959) 
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The  delayed  arrival  of  particulate  radiation  on  earth  is  due  to  1 

> 

transit  time.  An  example  is  the  cosmic  ray  burst  of  February  23,  j 

1956,  pictured  in  Fig.  14.  The  originating  flare  observed  in  He* 
apparently  started  at  03  hr  31  min  ±  1  min  UT.  The  first  cosmic  ray 
particles  arrived  at  earth’s  nightside  at  03  hr  45  min  ±  1  min.  The 
time  distribution  of  relative  increases  as  observed  at  Stockholm, 

Chicago  and  Berkeley  are  shown  in  the  figure.  The  earth’s  magnetic 
field  acts  as  a  spectrograph  separating  particles  of  different 
energies  according  to  geomagnetic  latitude.  The  energy  may  be 
parameterized  as  rigidity,  Z,  defined  to  be 

Z  =  *L£  1.2 

q  e 

wherein  qe  is  the  electric  charge  of  the  particle,  p  is  momentum,  and 

c  is  velocity  of  light  in  a  vacuum.  Z  is  usually  measured  in  units 

g 

of  (10)  eV  (i.e.,  GeV).  The  vertical  cutoff  rigidity  at  the  equator 
is  14.9  GeV;  only  particles  with  energies  greater  than  60  GeV  can 
approach  rhe  earth  from  any  direction.  Stormer’s  theory  has  shown 
that  the  cutoff  rigidity  for  vertical  incidence  is 

4 

Zc  =  14. y  cos  e  1.3 

where  e  is  the  geomagnetic  latitude.  Particles  with  a  rigidity  of 
somewhat  less  than  ZQ  can  reach  the  earth  at  latitude  e  at  large 
zenith  angles  and  at  certain  azimuths.  For  positively  charged  particles 
in  the  northern  hemisphere ,  the  minimum  rigidity  is  realized  for 
particles  arriving  from  the  western  horizon. 

The  path  traveled  by  a  particle  from  the  sun  to  the  earth  is 

not  radial.  One  hypothesis-;  pr^p^ by  Parker  tha*- _ 

the  ejected  particles,  having  been  ionized,  travel  along  magnetic 
field  lines  that  are  anchored  to  the  emissive  region  on  the  sun. 

Since  the  sun  is  rotating  counterclockwise  (see  Fig.  15),  its  rota¬ 
tion  produces  a  curvature  of  the  field  lines  which  intercepts  the 
earth’s  orbit  along  which  the  earth  moves  counterclockwise  (see 
Fig.  15).  The  angle  between  the  magnetic  field  path  of  particles 
and  the  radial  path  of  the  electromagnetic  radiation  can  be  computed 
in  the  following  manner 


19 


UT  (hrs) 

1111 _ 1 _ L_ 

100  200  500 

'(MIN) 


y  Burst  (After  Ehmert  and  Pfetzer,  1956) 
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6  =  tan 


-1  (cu  r) 


wherein  cu  is  the  angular  rotation  velocity  of  the  sun,  r  is  distance 
from  the  sun,  and  v  is  the  velocity  of  the  particle. 
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FIGURE  15.  Solar  Magnetic  Field  lines  as  Viewed  from  Sun’s  North  Pole 
(After  Parker,  1958) 


A  spectrum  of  solar  particles  energies  is  shown  in  Fig.  16 
which  also  indicates  the  integrated  intensity  of  particles  and  cosmic 
l-ays  of  galactic  (non-solar)  origin. 
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FIGURE  16-  Integrated  Energy  Spectrum  of  Solar  Particles  Associated 
with  Solar  Flares  (After  Obayashi  and  Hakura,  1960) 


ATMOSPHERIC  ENERGY  SOURCES 

The  sources  of  atmospheric  energy  are  the  difference  between  the 
solar-originated  energy  absorbed  by  the  atmosphere  and  the  energy 
emitted  by  the  atmosphere. 

At  the  distance  earth  is  from  the  sun  the  solar  flux  is  given  by 

■— 2  “  1  6 

the  so-called  "solar  constant"  (1.99  cal  cm  min  ,  or  1.388  (10) 

—  9  —1 

erg  cm  sec  )(Allen,  1963).  Of  this  energy  a  fraction  of  about  0.99 
is  absorbed  at  the  surface  and  by  carbon  dioxide  and  water  vapor  in 
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the  troposphere,  primarily  in  the  visible  range  of  the  spectrum.  The 
distribution  of  absorbing  water  vapor  and  of  reflecting  water  particles 
in  the  troposphere  is  strongly  variable  in  time,  latitude,  longitude, 
and  altitude. 

Above  the  troposphere,  the  distribution  of  absorbers  is  not  so 
strongly  variable.  A  measure  of  the  absorption  capabilities  is  the 
absorption  cross  section  a;  the  fraction  of  incident  flux  that  is 
absorbed  in  unit  distance  (dx)  of  a  radiation  path  in  which  the 
number  density  is  n  is  given  by 

dE  .  ,  _ 

— g  —  on  dx  1.5 

The  principal  constituents  of  the  atmosphere  below  200  km  are  atomic 
oxygen  (0),  molecular  oxygen  (02),  and  molecular  nitrogen  (N2). 

These  constituents  have  the  absorption  cross  sections  shown  in  Fig.  17. 

As  we  have  shown  in  Fig.  6,  the  hydrogen  Lyman -alpha  line  is  a 
very  intense  one.  It  penetrates  to  the  mesosphere  because  it  coincides 
with  a  minimum  in  the  wavelength  versus  absorption  cross  section  plot 
of  the  principal  absorber  molecular  oxygen.  This  minimum  appears  in 
the  illustration  of  the  absorption  cross  section  for  09 ,  shown  in 
Fig.  18- 

A  principal  absorber  in  the  stratosphere  is  ozone  O^.  (Its 
absorption  cross  section  is  shown  in  Fig.  19.)  It  may  be  noted  that 
the  absorption  cross  section  of  0^  at  the  wavelength  of  hydrogen 
Lyman-alpha  (1215.7  A)  is  larger  than  that  of  02 .  However,  in 
the  region  where  Ly-A  is  absorbed  (70  -  90  km),  the  number  density 
of  0^  is  not  large  enough  to  match  the  absorption  by  02  in  spite  of 
the  larger  absorption  cross  section. 

In  summary,  then,  the  solar  heat  sources  in  the  atmosphere  above 
the  troposphere  are  as  indicated  in  Table  1  which  shows  the  magnitude 
of  the  solar  flux  absorbed  in  the  altitude  regime,  and  the  spectral 
band  of  the  solar  radiation  from  which  the  heat  is  derived. 

In  Fig.  20  is  shown  the  rate  of  energy  absolution  of  X-rays  — 

220  to  900  a  and  1300  to  1800  a  wavelength  bands. 
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ABSORPTION  CROSS  SECTION  ( CM  ) 


ABSORPTION  COEFFICIENT  (CM  NTP) 


TABLE  1.  SOLAR  HEAT  SOURCES  IN  THE  ATMOSPHERE 
(After  Craig,  1965) 


REGION 

ION  REG 

ALTITUDE 

(KM) 

SOLAR  HEAT  SOURCE 

FLUX 

-2  -1 
(ergs  cm  sec  ) 

Thermosphere 

Upper 

F 

ISO  -  300 

Conduction  From  Corona 

0.25 

Lower 

E 

80  -  150 

Schumann-Runge  Cont.  1300-1750A 

EUV  (X<1000A)  From  Chromosphere 

30.0 

Mesosphere 

D 

50  -  80 

Schumann  Runge  Bands  17S0-1950A 

From  Chromosphere 

100. 0 

Stratosphere 

30  -  50 

Hartley  Continuum  2000-3000A 

From  Lower  Chromosphere 

1.8(10)4 

Troposphere 

0-30 

Absorption  From  Photosphere  And 
Re-Emission  From  C02»  H^O,  Surface 

1.37(10)S 

_1  _o 

RATE  OF  ENERGY  ABSORPTION  (erg  sec  cm  °) 

ST 11 -20-68-21 


FIGURE  20.  Rate  of  Energy  Absorption 
(After  Johnson,  1958) 
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In  addition  to  the  sources  of  heat  in  the  atmosphere,  it  is 
necessary  to  consider  the  sinks  of  energy.  Chief  among  these  is  the 
radiative  emission  of  energy  out  of  the  atmosphere.  In  Fig.  21  is 
shown  a  composite  plot  of  vertical  radiance,  as  measured  by  a  number 
of  observers  from  balloons,  aircraft,  rockets  and  satellites  in  a 
variety  of  observation  zenith  angles.  The  daylight  radiance  appears 
to  emanate  from  a  source  like  a  black-body  at  about  6000°  X,  and 
the  nighttime  radiance  as  from  a  black-body  of  about  300°  K.  An 
exception  to  the  latter  generalization  is  the  nighttime  Ly-A,  f-.:r 
which  there  may  be  several  explanations. 
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FIGURE  21 .  Vertical  Radiance  (Dayglow,  Nightglow) 
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In  Fig.  22  is  shown  a  schematic  representation  of  the  infrared 
spectrum  of  air.  The  curve  shows  black-body  emittance  at  300°  K  in 
relative  units.  The  hatched  areas  represent  the  spectral  regions  of 
appreciable  absorption  and  emission  due  to  CC>2,  H20,  and  0^.  In 
Fig.  23  are  given  the  cooling  rates  due  to  the  15  micron  C02  band 
and  to  the  9.6  ozone  band,  according  to  calculations  by  Plass  (1956) 
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FIGURE  22.  Infrared  Spectrum  of  Air  and  Spectral  Regions  of  Appreciable 
Absorption  Due  to  CC>2,  ^O,  and  O3  (After  Craig,  1965) 
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FIGURE  23.  Cooling  Rates  Due  to  15— CO2  Band  and  9.6-fl  O  Band 
(After  PI  ass,  1956)  3 


II. 


PHYSICS  OF  THE  TROPOSPHERE 


The  following  review  of  the  physics  of  the  troposphere  is  an 
attempt  to  summarize,  from  a  very  large  body  of  literature,  the 
features  of  the  troposphere  which  are  pertinent  to  observations 
made  from  satellites.  A  particular  need  is  to  relare  the  features 
by  some  sort  of  generalization  which  makes  possible  an  appreciation 
of  the  techniques  for  determination  of  tropospheric  parameters. 

Treated  first  is  an  effort  to  appraise  in  general  the  magnitude  and 
range  of  variation  ("scale")  of  important  atmospheric  parameters  in 
terms  of  spectral  distribution  in  frequency,  lifetime  and  dimension 
and  in  particular ,  the  scale  of  those  parameters  treated  in  numerical 
weather  prediction.  By  "numerical  weather  prediction"  is  meant  the 
forecasting  of  weather  for  more  than  24  hours  ahead.  This,  as  well  as 
the  analysis  of  the  current  weather  situation,  is  accomplished  by  the 
use  of  high-speed  electronic  computing  machines. 

Following  this  general  description  of  particular  features  of  the 
troposphere  is  a  review  of  some  important  features  of  the  hydrodynamic 
theory  of  wave  motion  in  the  troposphere.  Finally,  some  of  the  gross 
features  of  the  general  circulation  of  the  atmosphere,  involving  both 
friction  and  heat  sources  and  sinks,  are  described. 

GENERAL  DESCRIPTION  OF  PARTICULAR  FEATURES  OF  THE  TROPOSPHERE 

The  climate  that  we  experience  is  affected  by  a  variety  of  in¬ 
fluences.  Of  these  influences  on  climate  (some  of  which  have  been 
listed  by  Landsberg  [1945],  and  are  shown  in  Fig.  24),  a  dominant  one 
is  insolation  which  has  effects  listed  in  Tables  2  to  4.  At  the 
surface  of  the  earth  some  of  these  parameters  of  the  weather  have  had 
extreme  values  (see  Table  5).  The  scales  of  various  phenomena 

mCfDIHB  m  BUNK 
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TABLE  2.  ALBEDO  OF  SOLID  AMD  SMOOTH  SEA  SURFACES 


I  SOLID  SURFACES  | 

CULTIVATED  SOIL 

SAND 

■r wm 

OLD  SNOW 

TYPE  OF  SURFACE 

AND  VEGETATION 

AREAS 

AND  SEA  ICE 

PERCENT  REFLECTION 

7-9 

13-18 

1  gj  | 

50-60 

|  SMOOTH  SEA  SURFACE  ! 

SOLAR  ELEVATION,  DEGREE 

— 

5 

10 

20 

30 

— 

40 

50-90 

PERCENT  REFLECTION 

40 

25 

12 

6 

4 

3 

TABLE  3.  MEAN  DAILY  HEAT  SUM  RECEIVED  ON  CLEAR  DAYS 
IN  VARIOUS  LATITUDES  OF  THE  NORTHERN  HEMISPHERE 


LATITUDE,  DEGREE 

0 

15 

30 

45 

60 

75 

HEAT  SUM,  CAL/CM2/DAY 

510 

510 

470 

380 

300 

220 

TABLE  4.  HEAT  BUDGET  OF  THE  EARTH  AS  A  WHOLE  FOR  AVERAGE 

DAY  OF  THE  YEAR 


HEAT  RECEIVED  BY 

cal/cm2/day 

HEAT  LOST  BY 

cal/cm2/day 

DIRECT  SOLAR  RADIATION 

170 

OUTGOING  RADIATION 

790 

DIFFUSE  RADIATION 

140 

EVAPORATION 

120 

BACK  RADIATION  FROM  ATMOSPHERE  ' 

600 

TOTAL 

910 

total 

910 

AFTER  H.  E.  LANDSBERG,  IN  F.  E.  BETxY  ET  AL,  HANDBOOK  OF  METEOROLOGY  1945 
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TABLE  5. 


EXTREMES  OF  SOME  CLIMATIC  ELEMENTS  OBSERVED  ON  EARTH* 


ELEMENT 

VALUE 

LOCATION 

IN  THE  UNITED  STATES 

VALUE 

LOCATION 

ABSOLUTE  MAXIMUM  TEMPERA¬ 
TURE 

136°F 

AZIZIYA,  TRIPOLITANIA 

134°F 

GREENLAND  RANCH, 

DEATH  VALLEY,  CALIF. 

ABSOLUTE  MINIMUM  TEMPERA¬ 
TURE 

-I08°F 

OIMEKON,  SIBERIA,  US¬ 
SR 

-66°F 

RIVERSIDE,  YELLOWSTONE 
PARK,  WYO. 

ABSOLUTE  MAXIMUM  PRESSURE 

1,078  MB 

WESTER1 1  SIBERIA 

1,060  MB 

LANDER,  WYO. 

ABSOLUTE  MINIMUM  PRESSURE 

887  MB 

PACIFIC  OCEAN,  OFF 
PHILIPPINES 

892  MB 

LOWER  MATECUMBE 

KEY,  FLA. 

ABSOLUTE  MAXIMUM  RAINFALL  IN 
24HR 

46  IN 

BAGUIO,  PHILIPPINES 

25.83  IN 

CAMP  LEROY,  CALIF. 

ABSOLUTE  MAXIMUM  RAINFALL  IN 

I  MONTH 

241  IN 

CHERRAPUNJI,  INDIA 

HIGHEST  MEAN  ANNUAL  RAINFALL 

476  IN 

MT.  WAIALSALE,  T.H.** 

150  IN 

WYNOOCHEE  OXBOW, 
WASH. 

ABSOLUTE  MAXIMUM  RAINFALL 
MEASURED  IN  1  MIN. 

1.02  IN 

OPID'S  CAMP,  CALIF. 

ABSOLUTE  MAXIMUM  WIND  SPEED, 
ACCURATELY  MEASURED 

225  MFH 

MOUNT  WASHINGTON, 
N.H. 

MAXIMUM  MEAN  WIND  SPEED, 
ACCURATELY  MEASURED  FOR  5 
MIN. 

**• 

188  MPH 

MOUNT  WASHINGTON, 
N.H. 

MAXIMUM  MEAN  WIND  SPEED, 
ACCURATELY  MEASURED  FOR  1  HR 

173  MPH 

MOUNT  WASHINGTON, 
N.H. 

MAXIMUM  MEAN  WIND  SPEED, 
ACCURATELY  MEASURED  FOR  24 

HR 

129  MPH 

MOUNT  WASHINGTON, 
N.H. 

*  THE  VALUES  GIVEN  DO  MCI  INCLUDE  MEASUREMENTS  IN  THE  FREE  ATMOSPHERE. 

**  VALUE  BASED  ON  7  YEARS  OF  INTERRUPTED  RECORD.  A  COMPUTATION  BY  USING  THE  NORMAL  VALUE 
AT  A  BASE  STATION  WITH  LONG  RECORD  (  FORMULA  ON  P.  958  )  YIELDS  462  IN.  BOTH  VALUES  ARE 
DEFINITELY  ABOVE  THE  LONG  RECC7.D  VALUE  OF  428  IN.  AT  CHERRAPUNJI,  INDIA. 

•**  FOR  WINDS  THE  VALUES  GIVEN  FOR  MOUNT  WASHINGTON  IN  THE  COLUMN  FOR  THE  U.S.  EXTREMES 
ARE  THE  HIGHEST  THAT  ARE  ACCURATELY  MEASURED  AT  STATIONS  ON  THE  EARTH'S  SURFACE.  THEY 
ARE  KNOWN  TO  HAVE  BEEN  EXCEEDED  IN  THE  FREE  ATMOSPHERE. 


AFTER  H.  E.  LANDSBERG,  IN  F.  E.  BERRY  ET  AL,  HANDBOOK  OF  METEOROLOGY  1945 
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characteristics  of  tne  troposphere  are  essential  features  to  be  borne 
in  mind  in  any  analysis.  These  are  the  scales  of  the  weather  as  we 
know  it,  to  be  sure,  but  the  dimensions  of  the  phenomena  are  not 
always  clearly  understood. 


FREQUENCY  SCALES  OF  TROPOSPHERIC  MOTIONS 

The  characteristic  spectrum  of  atmospheric  motions  is  not  well 
known.  Its  range  of  frequency  may  be  comparable  to  its  range  of 
characteristic  lengths,  which  ranges  from  the  equatcr-to-pole  distance 
(about  10 1  m)  to  atomic  dimensions  (about  10  ±u  m),  a  total  range  of 
17  orders  of  magnitude.  Some  idea  of  wind  speed  spectra,  made  by 
measurement  of  wind  speeds  at  Brookhaven,  Long  Island,  and  Caribou, 

Maine,  from  towers  about  100  and  10  m  high,  respectively,  and  covering 
many  orders  of  magnitude  of  frequency,  has  been  described  by  van  der 
Hoven  (1357)  and  Oort  and  Taylor  (1969),  as  shown  in  Fig.  25-  The  peak 
of  power  at  periods  of  about  four  days  corresponds  perhaps  to  planetary 
wave  propagation  represented,  for  example,  by  the  propagation  of 
cyclones  and  anti-cyclones-  The  lesser  peaks  near  the  12  and  24  hour 
periods  correspond  to  the  semidiurnal  and  diurnal  tides.  The  peak  near 
one  minute  seems  to  correspond  to  the  motion  associated  with  small- 
scale  turbulence,  familiar  as  gusts  and  lulls  of  wind.  Oort  and  Taylor 
(1969)  have  also  shown  effects  of  aliasing  from  periods  of  one  minute 
(the  basic  averaging  period  of  rhe  wind  speed  of  rhe  wind  reports)  and 
the  two  hour  Nyquist  frequency  of  their  own  study,  with  the  assumption 
that  the  true  spectrum  at  high  frequency  is  the  dotted  curve  which  has 
been  drawn  to  correspond  with  van  der  Hoven  (1957)  data.  The  difference 
between  the  true  data  (the  dotted  curve)  and  the  apparent  but  erroneous 
spectrum  (solid  line)  is  a  dramatic  demonstration  of  possible  error  of 
using  one  minute  averaged  observational  data  as  inputs  at  two  or  three 
hour  intervals  for  purposes  of  numerical  weather  analysis  and  prediction. 


TIME  SCALES  FOR  ATMOSPHERIC  MOTIONS 

In  numerical  weather  analysis  and  prediction  it  is  important  to 
distinguish  between  large  and  small  length-scales  of  weather  phenomena. 


35 


In  the  past,  the  effort  has  been  chiefly  concerned  with  weather 
phenomena  of  large  dimension  on  the  supposition  that  the  small-scale 
activity  is  embedded  within  the  large-scale  effects  and  is  carried 
along  with  it. 


SPECTRUM  WIND  SPEED  AT  CARIBOU,  MAINE 
(AFTER  OORT  AND  TAYLOR,  1969) 


FIGURE  25.  Horizontal  Wind  Speed  Spectra 
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G.  D.  Robinson  (1967)  has  suggested  one  basis  for  identifying 
specific  time  scales  with  specific  space  scales.  He  notes  that  the 
Navier- Stokes  equation,  the  basic  equation  of  hydrogynamics,  is  in 
the  form 

^rj^s  — :  +  u  etcj  =  forces  +  vv2u,  etc.  (2.1) 

where  v  is  a  coefficient  of  viscosity  and  u  is  a  velocity  component. 
It  is  concerned  with  the  equation  of  motion  of  a  fluid. 

The  lifetime  of  the  "forms"  of  the  fluid  variations  may  be 
estimated  to  be  of  order 

L2D  -  L2/v  (2.2) 

if  L  is  a  characteristic  dimension  and  D  is  the  coefficient  of 
diffusion.  Then,  if  the  bulk  velocity  is  U,  the  time  for  a  form  to 
travel  distance  L  is  L/U.  A  form  cannot  travel  more  than  this 
distance  as  a  recognizable  entity  if 

2 

L/U  >  L  D  (prediction  impracticable) 


or  if 


LU  D  =-  L-n/v  <  1 

In  turbulent  conditions,  the  analog  of  Eq.  2.3  is 

O 

L/U  >  If /K.  (prediction  impracticable) 


(2.3) 


or 


LU/K  <  1  (2,a) 

wherein  K  is  a  lateral  eddy  diffusion  coefficient  for  momentum 
appropriate  to  the  length  scale. 
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The  scales  L,  K  and  v  are  not  independent  of  time;  however,  as 
an  example  by  which  the  predictability  time  may  be  estimated,  Robinson 
considers  the  following  model:  a  region  of  diameter  L  with  parabolic 
velocity  distribution,  zero  at  the  boundary,  maximum  U.  He  has  also 
taken : 


=  !u> 


“Hi 


(2.5) 


Taking  a  particle  to  be  accelerated  to  its  maximum  velocity  by 
external  forces  and  thereafter  retarded  by  lateral  diffusion  of 
momentum,  he  equates  average  inertial  and  dissipative  terms 

I  u  |Jl  (2-6) 
ay 


and  finds  that 


K=-if  (2.7) 

The  rate  of  dissipation  of  the  energy  of  the  large  scale  motion  is 


ay 


0  3 

16  U“  2  ir 

3  L2  3  L 


(2-8) 


By  deriving  U  from  (2.8)  in  terms  of  6,  and  by  substituting  into 
(2.7),  we  find 


1/3  1/3  4/3 
€  L 


(2.9) 


This  is  the  Richardson’s  law  (1926),  which  predicts  that  the  diffusion 
of  a  pair  of  particles  separated  by  a  distance  I.  has  a  dif fusivity 
K  proportional  to  L4//3.  The  numerical  coefficient  =  0.14 

obtained  here  is,  however,  at  varian,  e  with  that  proposed  by 
Richardson. 


The  limiting  values  of  the  frictional  dissipation  of  kinetic 

energy  tabulated  by  Robinson  are  1  <6  <  1G  cmc  sec  3 ,  established 

from  the  observed  radiation  field  without  regarr.  to  atmospheric  motion. 

3  2  -3 

Actually  other  observers  have  reported  values  or  6  ~  10  cm  sec  , 


7  r-\ 

30 


m, 


using  Eq.  2.9,  Robinson  found,  with  his  assumptions,  the  limits  of 
K  are 

0.15L4//3  <K  <  0 . 3L4^3  (2.10) 

Robinson's  values  of  the  eddy  diffusion  coefficient  K.,  bulk  velocity 

9 

U  and  predictability  time  T  -  L/U  =  L  /3  K  for  the  various  scales  of 
atmospheric  motion  are  given  in  Table  6.  For  values  of  dissipation 
rate  greater  than  those  assumed  by  Robinson,  predictability  times  are 
even  shorter  than  those  of  Table  6. 

Robinson’s  suggested  relationship  assumes  that  predictability 
time  (or  time  interval  over  which  a  state  is  predictable)  is  the  same 
as  the  lifetime  (or  the  time  over  which  a  state  may  be  expected  to 
exist).  That  they  are  the  same  is  controversial-  Also,  neither  is 
likely  to  be  the  same  as  the  scale  time,  which  for  phenomena  of 

c; 

interest  (length  scale  of  1000  km),  is  more  than  10'  sec,  say  about 
four  days,  as  described  in  later  paragraphs. 

SCALE  OF  ATMOSPHERIC  PARAMETERS  FOR  NUMERICAL  WEATHER  PREDICTION 

The  need  for  exact  knowledge  of  atmospheric  parameters  for  the 
purpose  of  numerical  prediction  of  weather  for  more  chan  24  hours  in 
advance  has  been  reviewed  by  Phillips  (1960).  In  his  review,  Phillips 
notes  that  0.  G.  Charney  (1948)  made  useful  estimates  of  space  and 
time  scales  of  the  meteorologically  important  motions  of  the  tropos¬ 
phere  which  may  be  numerically  simulated  with  electronic  computer 
computations. 

Charney  (1948)  and  Phillips  (196P)  have  been  guided  by  the 
experience  of  synoptic  meteorologists  who  have  found  that  the  weather 
producing  motions  of  the  free  atmosphere  can  be  characterized  as 
quasi-hydrcstatic,  quasi-aaiabatic,  quasi-horizoncal  and  quasi- 
geostrophic.  The  theory  is  based  on  a  kind  of  dimensional  analysis 
similar  to  that  used  in  th-  boundary  layer  theory  of  aerodynamics, 
in  which  the  motion  is  eft; s • terized  by  a  length  parameter  and  a 
velocity  parameter,  in  terms  of  whicn  the  orders  of  magnitude  of  the 
individual  terms  of  the  equations  of  motion  are  evaluated.  In  the 
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Charney  Treatment,  the  parameters  characterize  the  horizontal  and 
vertical  scales  of  motion,  the  speed  of  propagation  of  the  streamline 
pattern,  the  horizontal  particle  speed,  and  the  internal  static 
stability. 

Since  the  governing  equations  are  non-linear  and  the  small  scale 
motions  are  not  superposable  on  the  large  scale  motions,  it  is  not 
strictly  proper  to  regard  the  two  scales  of  motions  as  independent. 
However,  in  the  atmosphere,  since  the  bulk  of  the  energy  is  asso¬ 
ciated  with  the  large-scale  systems,  the  small-scale  motions  may  be 
regarded  as  turbulent  fluctuations  giving  rise  to  small  Reynolds’ 
stresses  and  heat  transports  which,  in  the  first  approximation,  may 
be  ignored. 

Charney ’s  analysis  begins  by  assigning  the  following  orders  of 
magnitude  for  the  large-scale  motions,  which  are  assumed  to  be  wave¬ 
like  in  shape: 


S 

H 


V 


C 

K 

K 


g 

f 


horizontal  scales  mean  horizontal  distance  between 
velocity  extremes  (~106  m) 

vertical  distance  scale,  mean  vertical  distance  between 

A 

velocity  extremes  (~10  m) 

mean  horizontal  velocity  (~10  m  sec-1) 

mean  speed  of  propagation  along  horizontal  streamlines 

(~10  m  sec  1) 

non-dimensional  parameter  characterizing  the  static 
stability  of  the  atmosphere 


(vd  -  v) 


where  =  dry  adiabatic  lapse  rate 


^  and 


y  is  the  actual  lapse  rate  in  the  atmosphere  (~10  -) 

_2 

acceleration  due  to  gravity  (~10  m  sec  ) 

the  z-component  of  the  earth’s  vorticity,  f  =  2-nQsina? , 

where  Q  is  the  angular  speed  of  the  earth’s  rotation 

-4  -1 

about  its  axis  and  $  is  the  latitude  (~10  sec  )' 
the  y-component  of  the  earth’s  vorticity,  j  =  gnncosl , 
(~10  4  sec  ■*■). 
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The  first  four  of  these  are  determined  by  observation  of  weather  naps ; 
values  of  K  are  given  for  normally  observed  lapse  rates  in  the  tropos¬ 
phere;  values  of  f,  j  hold  for  latitudes  between  0.5°  and  75°. 

The  basic  equations  for  the  motions  in  the  atmosphere  are  given 
by  the  equation  of  motion, 


dv  Bv 
dt  ~  dt 


+  v  .  Vv  =  -20  xv - gk+F 


(2.11) 


for  which  the  components  are  written  if  we  neglect  external  forces 
F,  and  viscous  forces  vt 


fv  +  jw  =  "i  |£ 


-  +  i  ^ )  -  fv  +  jw  =  -i 

:  V  3X  iy  3Z/  P 

l  +  (  W  )  f „  _  -1  3£ 

:  \  BX  +  By  Bz/  p  By 

i  +  (  +iu  +  .  =  il 

:  y  3x+  By  az/+]U  g  P  bz 


Bv  ,  I  uBv  ,  vav 

at  \  ax  ay 


The  equation  of  continuity  of  density  p 


(2.12) 


M  =  -v  .  vp 

and  the  equation  for  the  conservation  of  entropy  s 

as  _  q 

at  +  ~  *  vs  -  ^ 

In  Eq.  2.14,  the  atmosphere  is  considered  as  a  perfect  gas, 
s  =  cv  In  p  -  Cp  in  p  =  0^  tn  T  -  R  ^  p 


and  the  equation  of  state  is 

T  =  — £ 

Rp 


(2.13) 


(2.14) 


(2.15) 


(2.16) 


In  these  equations,  v  is  the  velocity  (three  dimensional  in  u,  v,  w) 
relative  to  the  rotating  earth,  k  is  the  vertical  unit  vector,  f  is 
the  c- component  of  earth’s  vorticity  2nsin|,  j  is  the  corresponding 
South  versus  North  (y  axis)  component  (2Qcoscp),  a  is  the  angular 
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speed  of  earth  rotation ,  cp  is  the  latitude ,  F  is  the  frictional  force 
per  unit  mass,  p  is  the  density,  p  is  the  pressure,  s  is  the  specific 
entropy,  q  is  the  rate  of  heating  per  unit  mass,  T  is  the  temperature, 
c„ 

>  v 

c 


and  c^  are  the  specific  heats  at  constant  pressure  and  volume, 
and  R  =  cp  _y 

Additional  equations  are  required  to  describe  the  flux  of  water  vapor 
and  the  condensation  process. 

A  convenient  measure  of  the  entropy  of  dry  air  is  the  potential 
temperature  which  may  be  defined  by  s  =  Cp-tn0  +  constant,  from  which 
R/c 

9  =  T  (2.17) 


•« 


If  the  right  side  of  the  Eqs.  2.11  to  2.14  is  expanded,  operators 
9  _ _ . _ 9 


like  — appear  in  addition  to  the  ~r  operator 
ox  oy  oz  ou 

assigned  the  order  cf  magnitude 


These  may  be 


a_ 

dX  ~ 

a 

-  3y  ~ 

|  =  TO"6  m  -1 

(2.18) 

3-. 

9z 

i 

J  H 

10"4  m  _1 

(2.19) 

9_ 

at ~ 

r>  a 

J  c  ax 

•  =  10*"^  sec-1 

(2-20) 

Given  these  orders  of  magnitude,  the  following  statements  may  be 
made: 

1.  The  continuity  Eq.  2.13  sets  an  upper  limit  on  the  value  of 
the  vertical  velocity,  which  may  also  by  the  geostrophic 
approximation  to  have  the  value 


W  ■ 


(HV  ,„-l  -l\ 

—  =  10  m  sec  J 


10  ^  m  sec-^ 


(2.21) 


2. 


Introduction  of  the  Inequality  relation  of  Eq.  2.21  into  the 
vertical  component'  of  the  equation  of  motion  (Eq.  2.11)  shows 
that  all  terms  are  smaller  by  a  factor  of  10~4  than  gravity 
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A 


(g)  and  the  vertical  pressure  force.  This  justifies  the 
use  of  the  hydrostatic  equation 

|f  =  -pg  (2.22) 

in  place  of  the  vertical  component  of  Eq.  2.11  in  which  the 
second  and  third  term  of  the  left  hand  side  effectively 
cancel  each  other. 

3.  Introduction  of  Charney’s  original  scales  and  the  inequality 
of  Eq.  2.20  into  the  horizontal  components  of  Eq.  2.11  shows 
that  the  horizontal  components  of  the  Coriolis  acceleration 
(~  fV)  and  the  pressure  forces  are  about  ten  times  as  large 
as  the  horizontal  acceleration  terms  C~  V^S).  In  this 
way;  the  geostrophic  relations  are  justified: 


(2.23) 


u  =  -  2£  ,  v=^^ 

pf  ay  *  pf  BX 

Here  the  differentiation  of  p  is  carried  out  on  a  surface 

of  constant  z,  which  yields  the  same  result  if  the 

differentiation  of  p  takes  place  on  a  constant  pressure 

surface  whenever  the  hydrostatic  relation  Eq.  2.21  is  valid. 

Phillips  notes  that  the  geostrophic  approximation  is  not 

as  good  as  the  hydrostatic  approximation  and  becomes  worse 

at  higher  velocities  (when  V/S  increases)  and  at  low 

latitudes  (when  f  decreases). 


4.  Introduction  of  the  hydrostatic  Eq.  2.21  and  the  geostrophic 
approximations  2.12  into  the  first  law  of  thermodynamics 
Eq.  2.14  and  disregarding  q/T,  as  though  the  process  were 
adiabatic,  leads  to  a  refined  estimate  of  W,  given  by  the 
equality  of  Eq.  2.21.  The  small  value  of  W  so  implied  means 
that  the  vertical  advection  terms  and  on  the  left 
side  of  the  component  Eq.  2.11  are  even  less  important  than 
the  horizontal  advection  terms  by  factor  -=j-g  and  should 
be  neglected  in  any  theory  based  on  Eqs.  2.22  and  2.23.  The 
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acceleration  of  the  horizontal  wind  may  be  computed  as  if 
the  motion  were  purely  horizontal.  A  further  conclusion 
stemming  from  the  inequality  (W/H)  <  (V/S)  is  that  the 
horizontal  divergence 


(2.25) 


Vp  .  (tu  +  jv)  ~  jj  (2.24) 

is  smaller  than  the  vorticity  of  the  horizonral  wind 

A  A  ^  ^ 

£  =  k  .  (  vp  X  (tu  +  jv)  ~  ^  (2.25 

By  a  slightly  different  analysis,  Burger  (1958) 
has  shown  that  the  inequality  (W/H)  <  (V/S)  becomes  an 
equality  when  S  is  as  large  as  10  7m,  comparable  to  the 
radius  of  the  earth.  This  invalidates  the  geostrophic 
theory  for  predicting  motions  of  such  great  horizontal 
extent. 

The  horizontal  increments  of  the  pressure  and  density 
fields  are  of  the  same  magnitude  as  their  mean  values, 
and  vary  in  the  same  horizontal  scale  S  as  do  the  velocity 
increments.  Similarly  the  vertical  increments  of  pressure 
and  density  vary  in  the  same  vertical  scale  H  as  do  the 
velocity  increments,  and  are  of  the  same  order  as  the 
pressure  and  density  itself.  Accordingly,  Charney  notes 
by  using  Eqs.  2.22  and  2.23 


gE  =  10  8  C2.26) 


■Lnp  ~  jj  —  10  5  ■knp  ~  ^  —  10 


(2.2?) 


Differentiating  Eq.  2.17  with  respect  to  x  or  y,  and  sub¬ 
stituting  Eqs.  2.26  and  2.27  yields 


fx  10'8s  iy  «•*>-§  =  1°'8  (2.28) 
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Then  using  Eq.  2.20  yields 


§  /  p  _  \  f  CV  „  -> 

St  — s?  ~  10 


(2.29) 


Applying  the  geostrophic  and  hydrostatic  Eqs.  2.22  and  2-23  and  the 
equality  10  giving  W  into  the  equation 


dv  d  v 
_ n  _  n  n 

dt  ”  dt 


(2.30) 


Charney  obtains 


du  dv  cv  fcv 
dt  ~  dt  ~  s  gHK 


10"4  (1  +  ^ 


(2.31) 


So  it  may  be  seen  that  the  term  in  dv^/dt  involving  w  is  one  order  of 
magnitud-  smaller  than  the  other.  Similarly,  Eqs.  2.22,  2-23,  and  the 
equality  of  10,  taken  with  Eq.  2.26  yields 


d  (top) 
Ilnp)  =  -S-gg- 


__  ..  w  (to?) 


(1  +  i)  =  10"  7  +  10"6j  ~ 


(2.32) 


The  implication  of  Eq.  2.32  is  that  the  individual  change  in  pressure 
is  due  almost  entirely  to  the  vertical  motion,  and  hence, 


Jvh 
dt  5=3 


(2.33) 


By  similar  reasoning,  Charney  shows 


d(tna )  wdtop 

dt  ~  bz 


(2.34) 


Using  Eqs.  2.18,  2.20,  2.21,  and  2.34  in  the  equation  of  continuity 
(Eq.  2,13  written  out  component  by  component)  yields 
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i«Ta»i 


Bu  cv  f CV 
dx  3y  ~  gHK 


10"G  sec"1 


(2.35) 


au 


■  |*  =  10 

ay  s 


5  sec-1  shows  that  the 


Contrasting  this  result  with  ^ 

terms  ~  and  comprising  the  horizontal  divergence  t-'nd  to  compensate 
ox  oy 

each  other,  explaining  perhaps  the  ever  present  difficulties  of  comput¬ 


ing  mutually  consistent  values  of  horizontal  divergence.  In  short, 
because  of  the  compensation,  it  may  be  impractical  to  compute  divergence 
at,  all. 


Charney’s  analysis  leads  from  the  hydrostatic  Eq.  2.22  and  the 
approximation  2.21  to  the  form  of  the  Eulerian  equations  most  often 
used  in  meteorology 


d,  v, 

nr-n  j., 

— r—  +  f  k  x  v 

CiC  Ml 


- 7,  p 

p  h 


i  ap 

P  Bz 


(2.36) 

(2. 37) 


These  scales,  developed  by  Charney,  are  summarized  in  Table  7. 


HYDRODYNAMIC  THEORY  OF  WAVE  MOTION  IN  THE  TROPOSPHERE 

A  most  powerful  tool  aiding  the  understanding  of  the  dynamics  of 
the  atmosphere  below  the  ionosphere  is  the  hydrodynamic  theory  of  wave 
motion.  In  addition  to  obviously  periodic  phenomena,  even  one-time 
only,  short-time  event  may  be  described  and  understood  in  terms  of 
wave  motion  theory.  Since  that  theory  is  rewarding  to  the  student 
of  atmospheric  physics,  it  will  be  discussed  at  some  length  in  the 
following  pages,  reviewing  the  basic,  hydrodynamic  equations,  linear¬ 
ization  of  the  equations  with  some  other  approximations  useful  in  the 
study  of  the  troposphere',  some  special  forms  of  the  equations  of 
limited,  but  useful  applicability,  and  analyses  of  typical  problems 
which  may  be  closely  identified  with  meteorological  phenomena.  This 
section  follows  the  development  of  the  review  by  Eliassen  and 
Kleinschmidt  (1957). 
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TABLE  7.  SCALES  OF  WEATHER  PREDICTION 


PARAMETER 

SCALE 

NOTE 

HORIZONTAL  DISTANCE 

S 

106m 

VERTICAL  DISTANCE 

H 

104  m 

HORIZONTAL  VELOCITY 

V 

10  m  sec  ^ 

HORIZ.  VEL.  GRADIENT 

V/S 

iO  sec  * 

/du  ,&v\  ln-6  -1 

(  — +r — 1~  10  sec 

\ax  ag/ 

HORP.  ACCELERATION 

CV/S 

in*4  ~2 

10  m  sec 

-~j  HORIZONTAL  CORIOLIS  FORCE 

VERTICAL  VELOCITY 

w 

ia-2  -1 

10  m  sec 

VERT.  VEL.  GRADIENT 

W/S 

in-5  -1 

10  sec 

i 

VERT.  ACCELERATION 

cw/s 

lA-8  -2 

10  m  sec 

10“7  ACCELERATION  OF  GRAVITY 

VELOCITY  OF  METEOR.  PROP. 

c 

“1 

10  m  sec 

VEL.  PROP.  OF  GRAV.  WAVES 

in2  -1 

10  m  sec 

VEL.  PROP.  OF  TIDAL  WAVES 

in2 

10  m  sec 

y-1 

LOG  POT.  TEMPERATURE  (8  =  RT 

P 

TIME  RATE  OF  LOG  POT.  TEMP. 

)  K 

fCV/gH 

.0-’ 

io'7sec'‘ 

STATIC  STABILITY 

3z  cv 

R  =  c  -  c 

P  v 

ACCELERATION  OF  GRAVITY 

9 

-2 

10  m  sec 

NEGLIGIBLE  FORAK63km 

HORIZ.  CORIOLIS  FORCE 

fV 

in-3  -2 

10  m  sec 

HORIZ.  EARTH  VORT1CITY  f=2irnsinp 

10  4  sec  1 

VERTICAL  CORIOLIS  FORCE 

jW 

tfr6  -2 

10  m  sec 

IO-4  ACCELERATION  OF  GRAVITY 

VERTICAL  EARTH  VORTICITY  j=2irn«js<p 

IO'4  sec'1 

LATITUDE 

(°-f)rod 

LOG  PRESSURE 

In  p 

1 

HORIZ.  GRAD.  OF  LOG  PRESS 

fV/gH 

Hf8  rr'1 

TIME  RATE  OF  LOG  PRESS  -££(  H 

gn 

b 

10'6sec'’ 

LOG  DENSITY 

In  p 

1 

HORIZ.  GRAD.  OF  LOG  DENSITY 

fV/gH 

10'8  m"1 

fry 

TIME  RATE  OF  LOG  DENSITY  --tH  1 
9!« 

10  3  sec  * 

TIME 

s/c 

105  sec 
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Basic  Hydrodynamic  Equations 


In  Cartesian  coordinates,  the  complete  set  of  equations  governing 
the  atmosphere  are  the  following. 


The  equations  of  motion 


If +  “  if +  v  t? + » it  = 2n  (vsini  -  -  ?  If +  s  F*  <2-3B> 


|^  +  u|~  +  v|^  +  w||-  =  -2Qusin$  -  -  |£  +  ~  P 
ot  ax  3y  az  o  Bv  m  3 


o  ay  m  % 


(2.39) 


+  +  +  =  “u=os5  -  i  |E  +  i  F  -  g  C2.<'» 

ot  ox  oy  a z  p  m  z  57 

and  the  equation  of  continuity  - 


§4  +  U  +  V  14 


at  ax 


»  If  -  -p(fe 


av  +  awl 
ay  az  i 


(2.41) 


These  are  the  four  equations  in  five  dependent  variables,  u,  v,  w,  p, 
and  p.  A  fifth  relationship,  an  "equation  of  state,"  is  usually 
derived  from  some  physical  restriction  on  the  system  based  on  the 
first  law  of  thermodynamics,  for  example: 


Tds  =  de  +  p  d 


t  -i 


=  C  dT  -  -  dp  =  C  dT  +  p 

p  p  v  dt 


(2.42) 


The  thermodynamic,  or  internal  energy,  of  dry  air  may  be 
described  by  the  equation 


§|  -  -P  -g^  +  6  +  jj§  (2.43) 

wherein  6  F,^,  Stokes  viscous  dissipation  function  and  D 

P  ~ 

represents  the  Eulerian  derivative 
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^  -  heat  received  per  unit  mass  and  unit  time 

e  =  CyT  +  constant  =  internal  energy  per  unit  mass 
of  dry  air 

6  =  entropy  per  unit  mass  of  dry  air. 

The  thermodynamic  equation  indicates  that  internal  energy  may  be 
increased  by  compression,  dissipation  and  heat  supply. 

As  long  as  the  effects  of  water  vapor  are  ignored,  the  thermo¬ 
dynamic-hydrodynamic  state  of  the  atmosphere  is  characterized  by  the 
velocity  field  v  and  the  two  scalar  fields  p  and  p .  The  change  with 
time  of  these  fields  is  described  by  the  above  equations  of  motior , 
continuity,  -and  thermodynamic  energy. 

The  effect  of  water  vapor  may  be  taken  into  account  by  intro¬ 
ducing  the  specific  humidity  as  another  dependent  variable,  for 
which  must  be  introduced  the  continuity  equation  for  the  water 
component,  and  a  suitable  energy  equation. 

Denoting  the  partial  pressure  and  density  of  water  vapor  by  pT 
and  p  ’ ,  its  equation  of  state  may  be  written 


wherein  R  is  the  gas  constant  of  dry  air,  and  a  =  0.622  is  the 
relative  density  of  water  vapor  with  respect  to  dry  air.  The 
enthalpy  of  water  vapor  is 

h’  -  Cp*T  +  constant 

and  of  liquid  water  is 

h"  =  c"T  +  constant 

wherein  c^’  is  the  specific  heat  of  water  vapor  at  constant  pressure;, 
c"  is  the  specific  heat  of  liquid  water,  and  the  latent  heat  of 
evaporation  is  given  by 

L  =  h*  -  h"  (2.45) 


50 


The  enthalpy  of  water  vapor  may  also  be  written 


h’  =  c"T  +  L  +  constant 
The  mixing  ratio  is  defined  as 


(2.46) 


(2.47) 


It  is  a  small  quantity,  of  the  order  1%  in  the  lower  troposphere 
and  still  smaller  at  greater  heights.  In  unsaturated  air  it  changes 
only  as  a  result  of  diffusion  and  turbulent  mixing,  a  slow  process. 

In  many  cases  the  mixing  ratio  may  be  considered  a  constant.  In 
this  case,  the  enthalpy  per  unit  mass  of  moist  unsaturated  air 

c  +  p,  Tc”  t 

h  =  -£=- - j—  T  +  - f  L  +  constant  (2.48) 

I  +  p,  JL+p, 


may  be  approximated  by 


h  ca  Cp  T  +  p.’L  +  constant  (2.49) 

and  the  energy  equation  for  unsaturated  moist  air  is 

TDs  =  Dh  -  i  Dp  (2.50) 

as  in  the  case  of  the  dry  air  energy  equation.  Thus,  nonsatui'ated 
air  behaves  in  this  approximation  in  the  same  way  as  dry  air.  In 
particular,  the  isentropic  change  of  state  will  be  very  nearly 
unaffected  by  the  humidity. 

The  case  of  saturated  air  is,  however,  different.  If  a  mass 
of  moist  air  contains  liquid  water,  so  that  one  mass  unit  of  dry 
air  is  associated  with  p,’  mass  units  of  water  vapor  and  p,"  mass  units 
of  liquid  water,  the  vapor  may  be  assumed  to  be  in  equilibrium  with 
the  liquid  phase  so  that 


p’  =  P(T) 


(2.51) 
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wherein  P(T)  denotes  the  saturation  pressure  of  water  vapor  with 
respect  to  liquid  water  with  a  plane  interface.  In  this  case,  the 
.nixing  ratio  is 


li '  -  a 


P(T)  P(T) 

-  P(T)  w  0  p 


(2.52) 


that  is,  the  mixing  ratio  in  saturated  air  is  no  longer  a  small 
constant,  but  is  a  function  of  both  p  and  T.  The  quantity  will 
change  as  a  result  of  condensation  and  evaporation. 

The  total  mass  of  the  water  component  per  unit  mass  of  dry  air 
is 

p,  =  p.  *  +  p  " 

In  the  absence  of  diffusion,  turbulent  mixing  and  precipitation, 
one  has  an  equation  of  continuity, 


=  o 

Dt 


(2.53) 


The  enthalpy  per  unit  mass  of  the  two-phase  system  is 


%  +  V-c  u  * 

h  =  - T  +  Y-tt— 

1  +  p  1  +  p 


U: - ,  L  + 

+  U 


constant 


(2.54) 


and  the  energy  equation  is 


TDs  =  Dh  -  —  Dp  «  (c  +  ~\m  -  (- 

P  \  P  P  dty 


+  c  i|)Dp 

P  / 


(2.55) 


The  difference  between  the  saturated  case  and  the  nonsaturated  case 
is  apparent  from  consideration  of  the  entropy  equations. 

From  the  above  considerations,  one  may  see  that  the  thermodynamic 
properties  of  atmospheric  air  are  such  that  nonsaturated  or  clear  air 
behaves  nearly  like  dry  air,  regards ^ss  of  humidity,  and  saturated  or 
cloudy  air  has  properties  differing  from  those  of  nonsaturated  air 
and  which  are  nearly  independent  of  the  content  of  the  liquid  water 
or  ice.  The  total  content  of  the  water  constituent  is  primarily 
important  as  a  means  of  determining  when  saturation  will  begin  and 
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end,  and  as  means  of  distinguishing  between  the  two  alternate  forms 
of  the  energy  equation. 

Complete  Linearized  Equations 

Since  the  complete  non-linear  equations  of  thermodynamics  and 
hydrodynamics  are  difficult  to  solve,  many  efforts  have  been  made  to 
linearize  the  equations  by  considering  only  small  departures  from  a 
state  of  static  equilibrium. 

The  complete  linearized  equations  of  motion,  derived  by  the  methods 
of  perturbation  theory  which  result  in  omission  of  the  advection  term 
u*v,  are  given  in  a  convenient  form  by  Eckart  (1960)  in  the  form: 

Su  *  o 

~  +  Co(v  +  T£)  P  -  1TQ  +  n  x  U  =  F  (2.56) 

The  equation  of  continuity  of  pressure  is 

pkP 

ft  +  Co  v.U  =  G  (2.57) 

The  equation  for  conservation  of  entropy  is 

|§=C.U=H  (2.58) 

Expressed  in  spherical  coordinates  (r,  9,  X),  where  r  is  the  distance 
from  earth  center,  X  is  longitude  measured  from  west  to  east,  and  0  is 
latitude  measured  from  equator  north,  the  equations  become 


co  || 

St  or 

+  Co  rP  -  n^q  -  ncoseu^ 

=  F 

r 

(2.59) 

5UX+  Co 

-  nsineu  +  ncoseu 

o  a.  y  r 

—  T? 

(2.60) 

St  rcos9 

Fx 

sue  Co  SP 
St  r  se 

+  nsinSU^ 

=  Fe 

(2.61) 

SP  ,  5 Ur 

it  +  U)17 

J  Co  sUX  ,  Co 

rcose  SX  r 

=  G 

(2.62) 

+  Ur 

=  H 

(2.63) 
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wherein  boundary  conditions  are 


Ur  =  0  at  r  =  r 


Ur  =  0  at  r  =  « 


(2.64) 


(2.65) 


and  the  variables  of  the  field  equations  are 


a  =  cp o' co)  4 


r  =  &  (log  pnc  )  + 


P  =  (PoCo^  4 


ffe  Po%] 


l  =  (pqV  k 


PQ(y  -  i) 

G - j- — 

^PoCo>2 

P1  H  -  (p°C°)  q 

1  H  "  /ds _\  -4 


N  =  4 


m  1 


2  o  dz 


(2.66) 


The  chief  error  introduced  by  the  linearization  is  the  elimination 
of  (.v  .  grad  v)  of  the  complete  equations.  In  the  perturbation 
theory  approximation,  in  which  v  =  0,  the  term  falls  out.  However, 
for  many  tropospheric  phenomena  involving  large  scale  motions,  the 
perturbation  theory  approximation  is  not  applicable,  and  the  non¬ 
linearity  must  be  taken  into  account. 

The  subscripts  o  and  1  imply  steady  state  and  first  order 
perturbations ,  respectively,  of  the  variables,  which  may  be  written 
as  functions  of  (r,  6,  X,  t) 


Velocity 


v  =  0 

~o 


Pressure 


Velocity  Change 


v  =  v 
~1  ~1 


P0  =  P0Cz> 


Pressure  Change 


Pi  =  P-j. 


(2.67) 


(2.68) 


-s**er*- 


Density 

Po 

Force 


f  =  ° 
~o 


l?o 

g  dz 


_3_  p 

RT  po 


(2.69) 


Density  Change 

Pi  =  Pi 

Force  Change 

f  =  p  V.G  ~  F.V  p 
~1  Ko  m  ° 

where  F  =  Frictional  stress  tensor 

Gm=  Tidal  force  potential  of 

moon  less  that  of  the  earth 

=  4  Kg  M  ^  (i  -  cos2e  ) 
z  D 

M  =  Moon  mass 

Kcp=  Gravitational  constant 

B  =  Moon-zenith  distance 
m 

t  ~  Distance  from  earth  center 

irth 
(2.70) 


D  =  Moon  distance  from  earth 
m 


Entropy 


so  =  Po(<;ptnTo  '  Rtn!‘o  + 


constant)  =  p0(cv-tnpQ  - 


c  -tno  +  constant) 
p  o 


center 


Entropy  Change 


S-,  =  p 


[  fgk  dP  DT  _  f  cLP  1  + 

\J  pdT  T  J  p2T  PJ 

p  r  A!s 1 1 6)dt 

°J  To  (2.71) 

q1  =  Energy  density  flux 

5  =  Frictional  dissipation  = 


—  V.v  as  0 


Latent  heat  of  sublimation 
entropy  change  indicated 
within  square  brackets 

L  -  latent  heat  of  sublimation 
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ct  =  Relative  density  of  water 
vapor  with  respect  co  dry 
air 

ct  =  0.622  (2.72) 

P  =  Saturation  pressure  of 
water  vapor 

P  =  P(T)  (2.73) 


Total  heat  accession 


Heat  accession  change 

^1  ~  ^c  +  %  +  ^conv  +  qcorp  + 


S?ad  ”  ^R 


(2.74) 


Heat  by  conduction 


%(z)  =  Po  v-<KvV  (2'75) 

K  =  Thermal  conductivity 


Heat  by  viscous  dissipation 


q^z)  =  up  % 
If 


|jl  =  Coefficient  of  viscosity 

w  =  Vertical  velocity  dz/dt 

L  =  Characteristic  length 

(2.76) 


Heat  by  convection 


%onv(z)  =  -Pcp  w  |f 


Heat  by  absorption  of  radiation 


-T 

<Wz)  =  Pq]C]C  €x.  nt(z,t)F(Uae  ^ 


X  i 


(2.77) 


rr~* 


I 


=  Optical  thickness  = 


X 


“  n  (z,t) 


z  COBa 


—  dz 


E,  =  Flux  at  top  of  atmosphere 


'X 

t 

cr  = 


a  = 


€  = 


Fhotoionization  cross 
section 

Total  absorption  cross 
section 

Solar  zenith  angle 

Efficiency  factor  for 
conversion  to  heat  (ergs 
photon-1)  (2.78) 


Heat  by  absorption  of  corpuscular  radiation 

qcorp(z)  as  1x1  qradCz)  where  €t 
ergs  per  particle 


X 

(2.79) 


Heat  loss  ly  radiative  emission 

qR(z)  =  Pc 


A 


B^(T)dv  + 


1.68(10)-18e 

o 

1  +  0.6e  —§■  +  0.2e 
o 


-325. 

T 


3 


Mass  absorption  coefficient 
Planck  blackbody  function  (2.80) 
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Other  Approximations  and  Special  Fcrms 

A  number  of  approximations  to  the  complete  equations,  other 
tlian  linearization,  are  commonly  used.  The  other  approximations 
usually  take  (v.r)v  into  account,  but  simplify  equations  to  eliminate 
modes  of  motion  usually  associated  with  gravity  waves. 

None  of  these  are  applicable  at  all  altitudes  and  latitudes 
except  as  rules  of  thumb  for  which  the  crudest  ones  -  the  geostrophic 
equation  for  wind  and  the  hydrostatic  equation  serve  at  mid- latitudes 
and  low  altitudes  with  about  20%  uncertainty. 

1.  Quasi-static  Approximation 

The  quasi-static  approximation  is  one  in  which 


a.  Vertical  acceleration  is  ignored  in  Eq.  2.40 


b.  Horizontal  component  of  Coriolis  force  is  ignored 
in  Eqs.  2.38  and  2.39 

Horizontal  f  =  2Qcos@  =  0 

Vertical  X  =  2osine 


becomes 


Dt 


d.  div^pU)  becomes 


divhc*V 


(2.81) 


If  the  approximation  =  0  is  applied  consistently,  it 
is  found  that  the  eigenfrequencies  (linear  response  to 
small  perturbations)  are  not  all  real.  Such  an  approxi¬ 
mation  completely  changes  the  character  of  the  solution. 
The  application  of  setting  Q  COS  e  =  0  is  also  in  many 
cases  done  carelessly.  If  done  at  all,  O  COS  6  must  be 
set  to  zero  in  all  equations,  not  only  in  Eq.  2.59,  where 
it  makes  the  most  difficulty  for  easy  solutions.  If  n 
COS  0  Ur  is  emitted  from  Eq.  2.59,  but  the  term  Q  COS  8  U, 
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retained  in  Eq.  2.58,  the  equations  cease  to  be  self-adjoint, 
that  is  eigenfrequencies  become  complex,  and  the  nature  of 
solutions  is  completely  altered. 

Although  in  many  circumstances  it  is  possible  to  justify 

jji» 

the  two  approximations  — ^  =  0  and  Q  COS  G  =  0,  the  former 
seems  likely  to  be  untrue  at  high  altitudes ,  and  the  latter 
is  certainly  untrue  in  the  tropics. 

2 .  Geostrophic  Approximations 

The  geostrophic  approximation  includes  the  quasi-static 
approximation  discussed  above,  and,  as  well,  that  the  Coriolis 
parameter  is  constant 

n  cos  0=0 

f  =  2n  SIN  9  =  constant, 

that  the  earth’s  curvature  is  neglected,  that  there  's  no  friction 
F  =  0  in  Eqs.  2.58,  2.59  and  2.60  and  that  all  processes  are 
adiabatic 

q  =  g  =  H  =  0  (2-82) 

in  Eqs.  2.58,  2.61  and  2.62.  The  Eqs.  2.58  to  2.62,  therefore, 
reduce  to 


=  1  =  o 

3x  r  39 


*£  =  .  I—  S£  = 
3y  rcose  3\ 


“20  sin  e  p  Cp0c0) 


=  -2Q  sin  6  p  U 


3z  dr 


g  p 


(2.83) 


(2.84) 

(2.85) 
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This  approximation  yields  handy,  easily  remembered  rules 
of  thumb;  the  hydrostatic  equation 


p  =  _  i|2 
v  g  az 

and  the  geostrophic  equation 


'i?  k  x  p(2  Q  sin  9)  vhp  k  X  p(2  n  sine)  VhT 


(2.86) 


(2.87) 


with  solutions  true  within  10  to  20%  error  in  mid-latitudes , 
and  in  high  latitudes. 

It  will  not  be  valid  at  high  altitudes  and  in  the  tropics. 

The  geostrophic  wind  blows  in  a  direction  normal  to  the 
horizontal  pressure  gradient  with  lower  pressure  to  the  left 
in  the  northern  hemisphere  (to  the  right  in  the  southern  hemi¬ 
sphere).  The  geostrophic  wind  vector  is  directed  along  the 
isobars  in  a  horizontal  surface,  or  what  is  equivalent  along 
the  contour  lines  of  an  isobario  surface. 

A  measure  of  the  error  of  the  geostrophic  wind  equation 
is  given  by  the  geostrophic  departure 


(2.88) 


In  the  absence  of  friction,  the  geostrophic  departure 
is  an  expression  for  the  horizontal  acceleration.  The  geo¬ 
strophic  departures  are  generally  of  the  order  1  to  10  m  sec  1 — 
hence,  on  the  average  are  10  to  20%  of  'the  wind. 

Differentiation  of  the  geostrophic  wind  equation  with 
respect  to  altitude  z  yields  the  Thermal  Wind  Equation 


k  x  2- 


=  _  2  x  2- 
f  lKXaz 


(ijxadp))  =  -  f  [i  *  fl  (p 


(2.89) 


'xxm 


- 


3-  Quasi-geostrophic  Approximation 

The  quasi-geostrophic  approximation  includes  the  geostrophic 
approximation 


dUr 


=  0 


dt 

QcosQ  =  0 
flsine  =  0 
£  =  o 

Q  =  0 

G  =  H  =  0 


in  Eq.  2.58 

in  Eqs.  2.58  and  2.59 

in  Eqs.  2.59  and  2.60 

in  Eqs.  2.58,  2.59  and  2.60 

in  Eqs.  2.58  and  2.62 

in  Eqs.  2.61  and  2.62 


(2.90) 


In  addition  to  the  geostrophic  approximations,  since  the 
horizontal  divergence  is  much  smaller  than  relative  vorticity, 
the  quasi-geostrophic  approximation  assumes  that  the  motion  is 
rotational,  not  solenoidal,  and  hence,  that  the  horizontal 
pressure  gradient  is  quasi-balanced  by  the  geostrophic  wind 
field.  The  equations  also  neglect  certain  terms  required  by 
conservation  of  vorticity. 

The  quasi-geostrophic  equations  are  useful  in  high  and 
mid- latitudes .  The  solution  of  the  equations  requires  both 
upper  and  lower  boundaries  as  well  as  lateral  boundaries. 

The  approximation  is  not  applicable  to  high  altitudes 
and  the  tropics. 

4.  Quasi-solenoidal  Approximation 

The  quasi-solenoidal  approximation  uses  all  the  quasi¬ 
static  approximations 
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dUr 

dt 


0 


Qcos9  =  0 
F  =0 

Q  =0 

G  =  H  =  0 

However,  it  partitions 


(2.51) 


V  =  V  +  V ■ 

'Hi 


so  that  v  is  the  solenoidal  part;  for  which 


V.vt  =  0 

p 


and  is  the  nonrotational  part;  for  which 


Vx  V,  =0 
P  ~€ 


(2.92) 


Since  Vp.v^  «  (Q  =  k.V  x  v)  for  large  motion  systems 
we  use 


v  =  v . 
~h 


(2.93) 


and  with  the  divergence  equation  (by  talcing  divergence  operation 
on  Eqs.  2.59  and  2.60)  so  that  divergence  is  not  produced,  to 
reach  a  double-Poisson  equation  equating  the  fields  of  the 
potentials  of  velocity  and  vorticity. 

5.  Barotropic  Approximation 

The  barotropic  approximation  is  a  two-dimensional  model 
represented  by  horizontal  (isobaric)  maps.  It  may  be  applied 
to  one  or  more  layers,  each  of  which  is  a  boundary  condition 
for  the  other. 
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half  the  grid  spacing.  The  sampling  at  high  frequency  and  small  grid 
spacing  must  be  done  before  filtering  in  order  to  avoid  deleterious 
"aliasing."  In  short,  the  measurements  should  be  instrumented  to 
sample  with  frequency  greater  than  one  per  minute  and  by  spacings 
closer  than  one  per  kilometer,  before  reducing  the  bandwidth  of  the 
data  contained  to  one  sample  per  period  of  computer  updating  and  one 
sample  per  100  to  300  km. 

As  Eliasson  and  Kleinschmidt  (1957)  have  pointed  out,  the  general 
solution  of  the  nonlinear  hydrodynamic  equations  is  handicapped  by  the 
following  considerations : 

1.  The  processes  are  nonlinear  in  their  "mature"  stage 

2.  The  solutions  of  nonlinear  equations  cannot  be  expressed  as 
a  superposition  of  eigensolution,  hence  they  cannot  be  dis¬ 
cussed  analytically,  but  only  by  numerical  or  geophysical 
methods . 

3.  The  solutions,  called  "numerical  forecasting"  or  "physical 
forecasting,"  was  first  attempted  by  L.  F.  Richardson  (1922), 
are  susceptible  to  over sensitivity  of  the  equation  of  motion. 
Margules  pointed  out  that  the  pressure  continuity  equation 
(Eq.  2.62)  would  require  an  accuracy  of  10  ^  sec  1  in  deter¬ 
mination  of  div  v,  which  is  not  attainable  from  wind  observa- 

~  -1  -2 
tions  (i.e.,  10  m  sec  per  kilometer  implies  frequency  10 

sec-1) 

4.  Internal  and  external  gravitational  sound,  inertial  waves, 
and  oscillations  do  not  correspond  to  long  period  waves  seen 
on  weather  maps,  and  introduce  "meteorological  noise"  (a  term 
by  Charney  (1948)).  The  "meteorological  noise"  is  evidence 
of  "aliasing"  or  folding  back  of  higher  frequencies  to  appear 
as  spurious  lower  frequencies. 

Four  procedures  may  be  considered  as  alternatives  in  seeking 
solutions  to  the  fundamental  Eqs.  2.59  to  2.63. 

1.  Calculate  the  motion  with  fictitious  noise  waves  included, 
and  smooth  the  solution.  This  method  is  laborious,  as 
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pointed  out  by  Eliasson  and  Kleinschmidt  (1957),  but  also 
will  yield  results  which  are  fallacious  because  of  aliasing. - 

2.  Change  the  form  of  prognostic  Eqs.  2.59  to  2.63  so  that  they 
are  no  longer  satisfied  by  the  noise  waves.  This  procedure 
is  the  one  followed  by  each  of  the  major  efforts  at  physical 
forecasting.  The  approximations  are  called  "filtering 
approximations,"  since  their  effect  is  to  filter  noise  f r(  m 
the  solution.  Each  requires  a  compromise  with  accuracy  or 
generality. 

3.  Before  time  integration  of  Eqs-  2.59  to  2.62,  the  initial 
values  of  pressure  and  wind  may  be  adjusted  to  each  other  so 
that  no  noise  waves  of  appreciable  magnitude  will  occur  in 
the  solutions.  This  is  tantamount  to  filtering  values  of  U 
and  P  so  that  they  retain  no  frequencies  or  wave  numbers  to 
produce  noise  in  the  solution.  This  step  is  a  desirable  one, 
but  in  itself  may  be  inadequate. 

4.  Compute  the  motions,  by  Laplace  or  Fourier  transformation,  in 
the  (id,  K)  domain  wherein  the  band  limits  of  the  data  and  the 
computed  variables  may  be  rigorously  controlled.  In  order 
that  the  results  of  such  a  laborious  process  may  be  appro¬ 
priately  realized,  it  will  be  worthwhile  to  design  the  (id,  K.) 
characteristics  of  the  input  data  to  avoid  the  errors  of 
aliasing  and  to  exclude  the  unwanted  frequencies  and  wave 
numbers.  Accordingly,  it  may  be  necessary  to  extend  the 
effort  to  the  design  of  instruments  for  meteorological 
observation. 

The  fourth  procedure  listed  above  has  received  small  attention. 

In  the  advent  of  more  capacious  and  faster  computers,  and  of  efficient 
algorithms  for  Fourier  transformations  (Cooley  and  Tukey,  1965),  this 
procedure  should  be  given  appropriate  consideration. 

Solving  the  Complete  Equations 

MacDonald  (1962)  has  described  the  form  of  the  eigen  solution 
equations  for  wave  solutions  of  Eqs.  2.59  to  2.63.  These  eigen 
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solutions  to  the  perturbation  relations  may  perhaps  be  numerically 
calculable.  MacDonald  assumed  solutions  of  the  form 


U  =  U(r, 
P  =  P(r, 
Q  =  Q(r, 


6)  e2CkX'u,t) 

6)  e2(kX-U)t) 
0)  e2(kX-aJt) 


(2.97) 


where  k  is  the  angular  wave  number  associated  with  longitude,  and  it 
is  recognized  that  r  and  9  coordinates  are  not  separable. 

To  obtain  a  partial  differential  equation  in  the  spatial  coordi¬ 
nates  (r,  0,  X)  of  form 


2  2  2 

A  +  2B  +  C  — S  +  lower  order  terms  =  0 

ar2  drde  ae2 


(2.98) 


The  character  of  the  equation  is  given  by  the  sign  of  AC-B  ,  which  is 
given  by  the  sign  of: 

+  N2Q2sin20]  =  u^-u^CN2-*}2)  +  N2fl2sin20  (2.99) 

For  frequencies  larger  than  n  (one  per  12  hr  =2.3  (10)  5  sec  X) 


o 

u)  >  N;  AC  -  B  >0  Elliptic  (i*e,,  2-D  potential  equation) 
2 

u)  <  N;  AC  -  B  <0  Hyperbolic  (i-e.,  equation  of  vibrating 

string) 

2 

u)  =  N;  AC  =  B  Parabolic  (i.e.,  equation  for  linear  heat 

conduction) 


(2.100) 


For  frequencies  less  than  n 

Elliptic,  if  v? (N2  +  n2)  <  u)4  +  N^sin2© 
or  if  N  >  n  and  <  f2sin0 


(2.101) 
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In  short,  for  case 

N  <  id  Differential  equation  is  Elliptic 

flsin©  <  u)  <  N  Differencial  equation  is  Hyperbolic 
0  <  u)  <  Qsin0  Differential  equation  is  Elliptic 

For  the  case  of  solutions  involving  characteristic  lengths  L  >  100  km, 
and  characteristic  periods  greater  than  10  min,  the  equation  is  el¬ 
liptic  for  id  <  flsinfl  and  hyperbolic  for  id  >  QsinB  (as  seen  in  Fig.  26). 


CO 

0.6-1. 5  rad  min  1  N 

(10  min)  1 1.6(10)  ^sec  1 

(12  hr)'1  2.3(  10)~5sec_1  £1 
Osin  8 

E  -  IMPLIES  ELLIPTIC 
p  _  •<  PARABOLIC 

H  -  -  HYPERBOLIC 

FIGURE  26.  Form  of  Differential  Equations  for  Atmospheric  Simulation 


3 

t 

i 

(2.102)  1 
)  . 
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The  solutions  of  the  hyperbolic  differential  equation  are 
analogous  to  those  of  the  equation  of  the  vibrating  string;  those  of 
the  elliptic  differential  equation  are  analogous  to  those  of  the  two- 
dimensional  potential  equation,  and  those  of  the  parabolic  form  to 
those  of  the  linear  heat  conduction  equation. 

In  solving  the  complete  equations  in  the  case  of  global  circula¬ 
tion,  one  may  take  advantage  of  the  notion  that  the  frequencies  of 
interest  are  in  every  case  smaller  chan  the  Vaisala  frequency,  N 

N2  =  (H  +  10  de9/kln)  (2.103) 

which  varies  between  0.6  <  N  <  1.4  for  0  <  z  <  140  km,  as  discussed 
subsequently  in  this  report. 

On  this  basis,  the  equation  may  be  solved  as  a  hyperbolic 
equation  for  cases  in  which  id  >  flsinS,  and  as  a  parabolic  equation 
for  m  =  sin6,  and  as  an  elliptic  equation  for  u)  <  Clsin9. 

The  hyperbolic  equation  may  be  solved  by  integration;  the 
elliptic  differential  equation  may  be  solved  differently  as  a  boundary 
value  problem  in  two  dimensions  and  the  parabolic  differential  equa¬ 
tions  as  -a  boundary  value  problem  in  one  dimension. 

Typical  Wave  Motions 


As  indicated  by  the  equations  for  the  hydrodynamics  and  thermo¬ 
dynamics  of  the  atmosphere,  heat  sources  and  sinks  are  necessary 
in  order  to  maintain  motions  in  the  atmosphere  against  frictional 
dissipation.  However,  heat  sources  and  friction  are  comparatively 
slow  effects  which  do  nc.  directly  lead  to  the  very  rapid  changes 
seen  on  weather  maps  in  disturbed  regions.  Hence,  although  heat 
sources  and  friction  are  of  vital  importance  in  creating  the  supply 
of  potential  energy,  they  are  usually  negligible  in  effect  during 
a  limited  time  interval  of  conversion  from  potential  to  kinetic 
energy.  In  the  limited  time  interval,  the  motion  may  be  treated  as 
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if  adiabatic  and  frictionless.  Most  of  the  work  in  dynamic  meteor¬ 
ology  has  been  done  on  the  basis  of  the  theory  of  small-scale  linear, 
frictionless  and  adiabatic  oscillations. 

Perturbation  theory,  applied  in  the  form  of  the  linearized 
equations  described  above,  may  be  used  to  study  the  small  scale 
oscillation,  and,  in  particular,  the  stability  properties  of  the 
atmosphere.  The  existence  of  any  one  discrete  eigenfrequency  which 
is  imaginary  or  complex  so  that  the  corresponding  eigensolution  will 
grow  exponentially  with  time  suffices  to  prove  that  the  state  is 
unstable. 

In  any  kind  of  oscillation  or  wave  motion,  the  energy  will 
fluctuate  between  kinetic  energy  and  some  other  energy  form.  It 
is  natural  to  classify  the  various  types  of  oscillations  and  waves 
according  to  the  energy  forms  involved.  In  this  respect  there  are 
three  fundamental  types  of  hydrodynamic  oscillations  and  waves: 

a)  Gravity  waves:  In  pure  gravitational  oscillations,  the 
energy  is  stored  as  gravity  potential  when  it  is  not  in  kinetic  form. 

b)  Acoustic  waves:  In  pure  compressibility  or  acoustic 
oscillations,  the  energy  is  stored  as  internal  energy  when  it  is 
not  in  kinetic  form. 

c)  Inertial  oscillations:  In  the  case  of  perturb  t ions  of 
steady  currents,  the  kinetic  energy  of  the  current  represents  a 
supply  of  energy  from  which  the  kinetic  energy  of  the  perturbation 
may  be  fed.  There  exist,  therefore,  pure  inertial  oscillations  in 
which  only  kinetic  energy  is  involved. 

Besides  these  three  fundamental  types  there  are  mixed  types, 
where  apart  from  the  kinetic  perturbation  energy,  two  or  three  forms 
of  energy  are  involved.  Actually  nearly  all  wave  motions  in  the 
atmosphere  are  of  such  mixed  character;  however,  when  one  energy 
storage  form  predominates  the  wave  is  described  as  being  predomi¬ 
nantly  gravitational,  acoustic  or  intertial  in  character. 

The  following  discussion  of  atmospheric  waves  is  along  the  lines 
developed  by  Eliassen  and  Kleinschmidt  in  1957. 


69 


Wave  motion  under  the  influence  of  Gravity  on  a  resting  earth 


Two  topics  may  be  described  without  reference  to  the  rotation 
of  the  earth:  static  stability,  and  combined  gravitational-compress¬ 
ibility  oscillations. 

Static  stability.  As  shown  by  Fjortoft  (1946),  the  criterion 
of  static  stability  can  be  derived  rigorously  from  the  energy  equation. 
The  second  time  derivative  of  the  energy  store  cannot  be  negative  if 
grad  cp  and  grad  sQ  have  everywhere  the  same  direction.  That  is 

ds 

£_  °  >  o  (2.104) 


at  all  levels  implies  static  stability. 

In  a  simpler,  but  less  exact,  concept  consider  an  air  mass  in 
equilibrium,  and  let  a  small  air  particle  at  the  level  zq  where  the 
pressure  is  pQ,  the  density  pQ  and  the  entropy  sQ  be  displaced  to  a 
level  z^  where  the  equilibrium  state  is  characterized  by  p^,  p^ ,  and 
s^.  It  is  assumed  that  the  displacement  of  the  particle  does  not 
disturb  the  equilibrium  pressure  field  so  that  the  particle  arrives  at 
the  level  in  a  state  characterized  by  p1>  p',  and  sQ.  As  a  result 
the  particle  is  then  exposed  to  a  vertical  buoyancy  force  per  unit 


F  =  -g 


(2.105) 


Since  the  density  may  be  considered  as  a  function  of  pressure  and 
entropy , 


dp  =  y  dp  -  p  ds 


(2.106) 


wherein 


Y  =  t -  ,  and  u 

“P 

c  KT 
°v 


(2.107) 
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The  particle  and  the  environment  have  the  same  pressure,  but 
different  entropy  this 

p  '  -  Pl  =  TX (so  sx)  =  p  g§  (zx  -  ZQ)  (2.108) 


and  hence 


F  = 


_  2k  <-§.  (z 

p  ds  ^Z1 


(2.109) 


Since  F  is  reckoned  positive  upward,  the  buoyancy  force  will  accelerate 

ds 

the  particle  towards  the  equilibrium  position  if  p  -r~  >  0,  and  from 

QZ 

this  position  if  this  configuration  is  not  met,  thus  indicating 
stability  in  the  first  case,  and  instability  in  the  second. 


As  a  measure  of  the  static  stability,  we  consider  the  restoring 
buoyancy  force  per  unit  length  of  vertical  displacement.  This 
quantity  has  the  dimension  of  (frequency)  squared.  Its  square  root 
is  termed  the  Brunt-Vaisala  frequency  N 


P  _  qp  ds  _  q  ds 
z,  -  z„  p  dz  c  dz 

i  o  p 


(2.110) 


Also  by  the  gas  equation, 

n2  -  a£fi.  /ot  _  8T/a  z\ 

T  l Dp  Bp73 zl 
For  dry  or  nonsaturated  air. 


1 *2=236=2 /si  +  2_\ 
0  8z  T  laz  +  c  l 

and  for  saturated  air. 


(2.111) 


(2.112) 


(2.113) 


Dry  or  nonsaturated  air  is  statically  stable  when  the  potential 
temperature  increases  with  height,  or  when  the  lapse  rate  of  temper¬ 
atures  does  not  exceed  g/cp  =  0.01  °C  nf1.  On  the  other  hand,  if 
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the  air  is  saturated,  then  tne  lapse  rate  for  moist  air  6  must  be 

m 

used.  The  moist-adiabatic  lapse  rate  is  a  function  of  p  and  T.  At 
low  temperatures, 

6m  ta  (at  low  temperatures  ) 

P 

6  =  0.0065  °C  nf1  (at  0  °C,  1000  mb) 

m  3 

6m  =  0.0044  °C  m-1  (at  20°  C,  1000  mb)  (2.114) 

The  static  stability  of  cloudy  air  is  seen  to  be  considerably  less 

than  that  of  clear  air  with  the  same  lapse  rate.  If  the  air  mass 

is  statically  stable  as  long  as  it  is  nonsaturated ,  the  mass 

will  be  unstable  if  saturation-  occurs .  Such  cases  are  usually  raf erred 

to  as  "conditional  instability." 

The  principal  difficulty  in  dealing  with  systems  where  conden¬ 
sation  occurs  lies  in  the  circumstance  that  the  location  of  the 
regions  of  saturated  air  will  depend  upon  the  motion  itself.  Con¬ 
densation  will  normally  take  place  within  the  ascending  currents, 
whereas  clear  air  prevails  in  regions  of  subsidence.  Combined 
gravitational-compressibility  oscillations . 

Waves  which  combine  the  effects  of  gravity  and  compressibility 
are  important  in  the  troposphere  in  that  they  account  in  general 
for  the  short  periodic  fluctuations  of  winds — the  gusts  and  other 
wind  phenomena  that  seem  to  have  phase  coherence  as  opposed  to 
turbulence  that  lacks  it.  Eckart,  1960,  has  given  a  careful  discussion 
of  the  forms  of  these  waves.  The  following  description  of  the 
fundamental  solutions  of  the  wave  equations  for  an  isothermal  atmo¬ 
sphere,  taken  after  Eliassen  and  Kleinschmidt  (1957),  greatly  over¬ 
simplify  the  problem,  but  illustrate  the  general  nature  of  the  waves. 

For  an  isothermal  atmosphere,  the  coefficients  N2,  r  and  cQ2 
are  taken  as  constant: 
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(Brunt- Vaisala  frequency) 


i 

! 

i 

i 

\ 


(Large  scale  stability  parameter) 


(Speed  of  sound)  (2.115) 


In  the  atmosphere,  the  above  coefficients  (Eq.  2.115)  are  not  constant, 
but  have  values  similar  to  those  illustrated  in  Fig.  27.  The  perturba¬ 
tion  equations  of  motion  are  satisfied  by  solutions  of  the  form, 
wherein  |  and  £  are  displacements  in  the  x,  and  z  directions. 
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FIGURE  27.  Dynamic  Characteristics  of  Field  Equation  Variables 
(After  Eckart,  I960) 
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and  where  A  is  an  arbitrary  amplification  factor,  and  k,  m  and  v 
are  three  constants  satisfying  the  dispersion  relation 


(2.119) 


The  constants  k  and  m  are  the  horizontal  and  vertical  wave  numbers, 
respectively  and  v  is  the  frequency. 

The  different  types  of  fundamental  solutions  may  be  conveniently 
represented  in  the  frequency-wave  number  (v,  k)  diagram  for  the 
isothermal  planar  atmosphere  of  Fig.  28. 


The  horizontal  line  v  =  N  represents  purely  vertical  oscillations 

.  .  2  2  2 

in  an  undisturbed  pressure  field.  The  line  k  =  cQ  v  represents 

Lamb  waves  whose  horizontal  phase  velocity  (-r-)  equals  the  speed  of 

2  2  K 

sound  c  .  For  oiven  values  of  k  and  v  ,  the  dispersion  relation 

°  2 
(Eq.  2.119)  gives  a  corresponding  value  of  m  (square  of  the  vertical 

2 

wave  number).  The  curves  m  =  constant  are  a  family  of  hyperbolas. 

2  2 
The  hyperbola  m  =0  separates  the  domain  of  external  waves  (m  <  0) 

which  have  finite  amplitudes  at  the  boundaries  (such  as  near  the 

2 

earth’s  surface)  from  the  domain  of  internal  waves  (m  >0),  for 
which  the  amplitudes  at  the  boundaries  are  zero.  Tie  internal  waves 
may,  therefore,  be  subdivided  into  low  frequency  internal  waves 
(lower  region)  which  are  dominated  by  the  influence  of  gravity  and 
hence  are  called  "internal  gravity  waves,"  and  high  frequency 
internal  waves  (upper  region)  which  are  dominated  by  compressibility 
terms,  hence  called  "internal  acoustic  waves." 

Any  of  the  eight  combinations  ±k,  ±m,  ±v  will  satisfy  the 
dispersion  relation  (Eq.  2.119).  Linear  combinations  of  the  eight 
solutions  give  the  various  types  of  possible  real  solutions.  These 
include  the  following: 
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FIGURE  28.  Diagnostic  Diagrams  (Wr  k) 
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2  2 

A.  External  waves  for  which  -  m  =  >0 

1.  Standing  oscillations  may  be  represented  by 


pQ^£  =  A  sin  lex  e±mtz  sin  vt 


(2.120) 


2.  Waves  with  horizontal  propagation  may  be  represented  by 


p  £  =  A  e±miz  sin  (kx  +  vt) 


(2.121) 


By  inspection  of  the  term  e±mtZ,  it  may  be  inferred  that  external 
waves  cannot  occur  in  a  layer  of  air  confined  between  two  rigid 
horizontal  planes,  nor  as  free  waves  over  a  level  ground  in  an 
infinite  atmosphere  of  uniform  stability.  External  waves  may  exist 
as  forced  waves ,  produced  by  motions  of  the  boundaries ,  as  in  the 
phenomenon  of  mountain  waves,  or  as  free  waves  in  an  atmospheric 
layer  which  is  bounded  by  other  layers  with  different  properties. 
Propagation  of  external  wave  energy  is  transported  horizontally 
in  the  direction  of  propagation,  but  this  energy  cannot  be  transported 
in  a  vertical  direction. 

o 

B.  Internal  waves,  for  which  m  >  0 

3.  Standing  oscillations,  represented  by 


=  A  sin  kx  sin  m  z  sin  vt 


(2.122) 


The  wave  energy  flux  for  this  foim  vanished  when  averaged  over  one 
period.  The  low  frequency  type,  called  gravity  waves  and  modified 
by  compressibility,  represents  a  circulatory  motion.  The  high  fre¬ 
quency  type,  called  compressibility  or  acoustic  oscillations  and 
modified  by  gravity,  represents  a  more  complicated  streamline  pattern 
of  different  kinematics. 

4.  Waves  with  horizontal  propagation  may  be  represented  by 

p  =  A  sin  m  z  sin  (kx  +  vt)  (2.123) 

o 
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The  wave  energy  has  a  mean  value,  averaged  with  respect  to  z  and  x 
(or  t),  which  is 


Q  = 


A2.  N^k2  -  c02y4 

4  ^  "  c2v2 

o 


(2.124) 


The  mean  wave  energy  flux  is 


p  _  £  k  (N2  -  v2) 

x  4  v  ' 


2  2  ’ 
t  +  m 


(2.125) 


That  is,  the  wave  transports  energy  in  the  direction  of  propagation 
5.  Waves  with  vertical  propagation  may  be  represented  by 


p  2£  =  A  sin  kx  sin  (mz  +  vt) 


(2.126) 


The  mean  wave  energy  (averaged  with  respect  to  x  and  z,  or  t)  equals 
that  of  horizontally  propagated  waves  (Eq.  2.124).  The  mean  wave 
energy  flux  is: 


Fx  =  0, 


.  £  m  v2(N2  - 
4 


v 


2  2 
-L  +  m 


(2.127) 


Comparison  of  the  sign  of  the  vertical  energy  flux  with  that  of  the 

vertical  phase  velocity  (-  v/m)  indicates  that  the  signs  are  equal 

2  2 

for  high  frequency  internal  waves  (v  >  N  )  and  opposite  for  low 
frequency  internal  waves  (v2  <  N2).  Thus,  low  frequency  internal 
waves  (so  called  "gravity"  waves)  with  vertical  propagation  possess 
the  unusual  property  that  they  transport  energy  in  a  direction 
opposite  to  the  direction  of  propagation. 

6.  Waves  with  tilting  phase-lines,  represented  by 


p  z£  =  A  sin  (kx  +  mz  +  vt) 


(2.128) 


77 


These  waves  have  a  horizontal  as  well  as  a  vertical  propagation. 

The  mean  wave  energy  is  twice  the  value  given  by  Eq.  2.124.  The 
mean  horizontal  wave  energy  flux  is  twice  the  value  given  by  Eq.  2.125. 
The  mean  vertical  energy  flux  is  twice  the  value  given  by  Eq.  2.127. 

The  vertical  energy  propagation,  opposite  to  the  vertical  velocity, 
occurs  again  when  v  ‘  <  N"1 . 

Axially  symmetric  wave  motions  in  a  circular  vortex 

In  the  vortex  motions  of  the  atmosphere,  inertial  stability 
considerations  must  be  taken  into  account  in  addition  to  the  static 
stability  considerations  described  in  the  previous  section. 

Pure  inertial  stability.  A  revolving  mass  of  fluid  possesses  a 
particular  kind  of  stability,  usually  referred  to  as  inertial 
stability,  first  found  by  Rayleigh. (1916),  then  treated  by  V.  Bjerknes 
and  Solberg  (1929) . 

The  circulation  c  of  a  circle  R  =  constant,  z  =  constant 
of  a  homogeneous  and  incompressible  fluid  rotating  around  an  axis 
with  no  external  force  is 


c 


=  2ttR2w 


(2.129) 


where  R  is  the  distance  from  the  axis  and  z  is  the  distance  from  a 
plane  normal  to  the  axis.  The  kinetic  energy  of  the  vortex  bounded 
by  a  rigid  surface  of  revolution  symmetric  with  respect  to  the  axis 
and  enclosing  a  volume  t  is 


K. 


J  8u2R2 


dT 


T 


(2.130) 


Suppose  that,  in  addition  to  the  vortex  motion,  there  is  also  an 
axially  symmetric  meridional  motion,  characterized  by  a  meridional 
velocity  field  v^,  The  total  kinetic  energy  of  the  fluid  is 

Km  +  =  constant 
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wherein 


K 


m 


■/ 


2 

V  mpdT 


(2.131) 


The  field  corresponding  to  an  extremum  of  Kv  may  be  obtained  by 
requiring  the  variation  to  vanish  for  a  virtual  meridional  dis¬ 
placement,  so  that 


pcQ  grad  +  grad  pq  --  0 


(2.132) 


The  above  equation  shows  that  the  pressure  increases  with  R  and  is 
independent  of  z.  So  also  does  cQ  depend  only  on  R,  being  inde¬ 
pendent  of  z.  The  balanced  vortex  is  stable  with  respect  to  axially 

p 

symmetric  motions  when  Ky  is  a  minimum;  that  is,  when  cq  increases 
monotonically  with  R.  The  criterion  of  inertial  stability  may, 
therefore,  be  written 


d  CQ  /d  R  >  0 


for  all  values  of  R 


(2.133) 


The  stability  criterion  can  also  be  derived,  in  a  less  rigorous 

manner,  by  means  of  the  particle  method.  The  particle  is  chosen  as 

a  symmetric  circle  which  is  assumed  to  expand,  or  contract,  without 

2 

change  in  circulation,  through  an  undisturbed  environment.  If  cQ 
increases  with  R,  then  a  symmetric  circle  which  has  been  brought  to 
expand  will  arrive ■ in  z  environment  which  rotates  faster  than  the 
circle  itself.  The  pressure  gradient  which  balances  the  centrifugal 
force  of  the  environment  will  exceed  the  centrifugal  force  of  the 
displaced  circle  and  thus  give  a  net  restoring  force ,  which  may  be 
termed  a  "dynamic  buoyancy."  The  vortex  is  thus  inertially  stable, 
and  the  restoring  force  per  unit  mass  and  per  unit  length  displacement 
normal  to  the  axis  is  found  to  be 


2  2 
v  =  grad  c 
t  o 


(-  grad  m  )  = 


2n2R3 


dc 
_ o 

dR 


(2.134) 
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This  quantity  has  the  dimension  of  a  frequency  squared  and  may  be 

o 

taken  as  a  measure  of  the  inertial  stability.  If  cQ  decreases  with 
R,  the  dynamic  buoyancy,  acting  on  an  expanding  or  contracting  circle, 
will  be  directed  away  from  its  equilibrium  position,  thus  indicating 
inertial  instability. 

In  an  inertially  stable  vortex  an  axially  symmetric  disturbance 

of  the  balanced  state  will  result  in  inertial  oscillations.  In  these 

motions,  the  meridional  displacements  satisfy  differential  equations 

analogous  to  the  perturbation  equations  for  two-dimensional  gravi- 

2  2 

tational  waves,  with  v  instead  of  v  .  For  inertial  oscillations, 

*  t  s 

the  displacements  normal  to  the  axis  are  the  ones  which  give  rise  to 
restoring  forces,  while  the  displacements  parallel  to  the  axis  are 
inactive  and  must  be  driven  by  the  pressure  force  alone. 


Combined  effect  of  static  and  ineruial  stabilit 


The  stability  with  respect  to  axially  symmetric  motions  of  a 
compressible  baroclinic  vortex  in  a  gravity  field  is  an  example  of 
the  combined  effects  of  static  and  inertial  stability. 

The  energy  equation  of  the  system  is 


/  Kv  -  $  -  E  =  constant  (2.135) 

wherein  is  the  kinetic  energy  due  to  the  meridional  motions,  Kv 
is  that  due  to  the  vortex  motion,  5  is  the  gravity  potential  energy  of 
mass,  and  E  is  the  internal  energy  (all  of  these  are  energy  per  unit). 

The  motion  is  such  that  a  fluid  line  which  at  one  time  forms  a 
symmetric  circle  around  the  axis  will  always  remain  a  symmetric  circle. 
Its  circulation  c  as  well  as  its  specific  entropy  s  will  remain  un¬ 
changed.  Consequently,  the  change  of  the  quantity  J  =  Kv  +  $  -  E 
will  depend  only  upon  the  meridional  displacements  6rm.  J  may  be 
considered  an  "effective  potential  energy"  or  energy  store  for  the 
meridional  motion.  Differentiation  of  J  with  respect  to  time,  and 
the  continuity  of  potential  cp,  implies  that  J  will  be  an  extremum 
when 
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grad  t p  +  cq2  grad  cpv  +  j  grad  pQ  =  0 


The  second  variation  of  J  is,  given  by 


62J  = 


J  /  [6^  ‘  M  *  62m  +  rff  ^ 

c  RT  grad  p  \2"j 

+  -2- —  sr  .  - — )  p  dr 

cv  ^  P  /  J  ° 


-t —  6r  . 

cv 


(2.136) 


(2.137) 


where  M  is  a  tensor  defined  by 

1  2 
M  =  — 1  grad  sq  (~aQ  grad  pQ)  +  grad  cq  (-grad  cpv)  (2.138) 

CP  ''  °  ° 

2 

If  the  quadratic  form  6^  .  M  .  6r^  is  positive  definite,  then  6  J 
is  positive  for  all  virtual  axially  symmetric  meridional  displace¬ 
ments,  J  is  a  minimum,  K  is  bounded,  and  the  vortex  is  stable  with 
respect  to  axially  symmetric  disturbances. 


The  criterion  for  a  stable  vortex  is  then 


Det(M) 


1  2 

=  —  grad  s  x  grad  c  .  grad  cp  x  (-grad  cn  )  >  0 


Cp  -  o 


(2.139) 


Trace (M)  =  ~  grad  sQ  .  I - 
CP  \ 


grad  p,. 


+  grad  cQ  .  (-grad  cpv)  >  0 


(2.140) 


Stability  criteria  for  the  baroclinic  vortex  applied  to  the  atmosphere. 

In  an  absolute  frame  of  reference,  the  motion  of  the  atmosphere 
is  nearly  a  circular  vortex  motion  around  the  earth’s  axis.  Consider 
the  earth  as  a  sphere  with  radius  a.  In  a  point  with  latitude  $  and 


height  z  above  the  earth’s  surface, 


the  distance  from  the  axis  is 


R  =  (a  cp  z)  cos? 


(2.141) 


and  the  absolute  circulation  of  a  symmetric  circle  is 
c  =  2n  (QR2  +  Ru) 

a 


(2.142) 


where  Q  denotes  the  angular  speed  of  rotation  of  the  earth,  and  u 
the  relative  zonal  wind  speed,  assumed  independent  of  longitude. 

Under  these  circumstances,  the  criteria  for  stability  may  be 
written 


(Stable) 


(2.143) 


where  the  index  s  means  differentiation  at  constant  s,  i.e.  in  an 
isentropic  surface,  also,  it  may  be  written,  if  z  is  neglected  with 
respect  to  a,  as 


sin  ? 


•K«),  ”•-=(«),] 


>  0  (2.144) 


(Stable) 


Since  the  slope  of  the  isentropic  surface  is  very  small,  the  term  in 
cos?  may  also  be  neglected,  so  that  the  criterion  becomes 


1  + 1  *  -  (ly)s  >  0 


(Stable ) 


(2.145) 


If  use  is  not  too  great,  the  second  term  may  be  neglected  so  that 

\  -  >  0  Stable  (2.146) 

In  this  form,  the  criterion  applies  strictly  to  a  straight,  parallel 
current  of  arbitrary  direction  on  a  flat  earth,  rotating  with  angular 
velocity  Q  sin  § . 
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It  is  noteworthy  that  the  direction  of  least  stability  is  nearly 
along  the  isentropic  slope.  The  result  is  commonly  superimposed 
meridianal  circulations  in  flat  cells  oriented  along  the  isentropic 
slope . 

Application  of  the  above  criterion  to  atmospheric  currents  has 
indicated  two  kinds  of  flow  patterns  where  the  stability  with  respect  f 

j 

to  meridianal  or  transverse  circulations  is  normally  small  and  where  \ 

I 

the  current  may  be  unstable  in  pronounced  cases.  The  first  type  is  • 

the  frontal  zone  in  the  lower  troposphere,  in  which  the  unstable 
direction  coincides  with  the  frontal  slope.  The  significance  of 
the  instability  is  that  the  air  makes  little  resistance  against  an 
upgliding  motion  in  the  direction  of  the  frontal  slope.  The  other 
type  of  flow,  where  the  stability  is  small,  is  found  in  the  upper 
troposphere  to  the  south  of  strong  jets.  Here  the  isentropic  surfaces 
are  almost  horizontal,  so  that  there  is  little  resistance  just  south 
of  a  jet  stream  maximum  against  meridianal  or  transverse  circulations 
in  flat  horizontal  cells  superimposed  upon  the  current.  In  all  other 
parts  of  the  atmosphere  a  pronounced  stability  with  respect  to  meridi¬ 
anal  motions  seems  to  prevail. 

Two  dimensional  wave  disturbances  of  linear  currents 

If  a  mean  current  u  is  defined  by  averaging  with  respect  to  x, 

then 


u  =  u  +  U*  w  =  w*  w  =  0  (2.147) 

the  transport  of  momentum  through  a  line  z  =  constant  is,  apart  from 
a  density  factor, 


w„  =  uTwT 


the  change  in  mean  velocity  is 


Bu 

at 


a_ 

az 


u  *  w* 


(2.148) 


(2.149) 
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and  the  rate  of  change  of  kinetic  energy  of  the  mean  flow  is 


-H.fi 


(2.150) 


so  that  integration  through  the  layer  gives  the  total  kinetic  energy 
of  the  mean  motion 


=  -/ 


rz  2 

=  I  ufwl 


dz  >  0  (Stability) 


(2.151) 


For  amplifying  perturbations,  d  Km/d  t  must  be  negative,  that  is, 
the  momentum  transport  must  be  directed  from  layers  with  large  u 
values  to  layers  with  small  u  values. 

Stability  of  zonal  currents  in  the  case  of  horizontal 

non-diverqent  motion 

Fjortoft,  R.  (1950)  has  generalized  the  theory  of  stability  of 
linear  currents  so  that  it  applies  to  two-dimensional  non-divergent 
motion  on  a  sphere.  This  case  is  of  interest  in  meteorology,  since 
at  a  level  near  500  mb,  the  large  scale  motions  may  be  approximated 
as  horizontal  and  non-divergent. 

Consider  the  earth  as  a  sphere  with  radius  a,  and  latitude  §. 
The  area  of  a  ring  between  $  and  $  +  d$  is 


dF  =  2na  cos$d$  =  2TTaRd§ 


(2.152) 


The  absolute  circulation  along  a  latitude  circle  is 


C  (§)  =  2ttRu 

3  '"'a. 


(2.153) 
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where  a  bar  designates  mean  values  over  a  latitude  circle.  The  two- 
dimensional  velocity  field  defines  a  scalar  vorticity,  i.e.  the 
absolute  vorticity  component  £  normal  to  the  sphere,  which  is 
individually  conserved. 

From  theorems  of  Stokes  and  Kelvin  for  the  conservation  of 
circulation 


£^dF  =  constant  (2.154) 

F, 

X 

From  the  principle  of  conservation  of  total  absolute  angular  momentum 

Ff^dF  =  constant  (2.155) 

From  Eqs.  2.154  and  2.155,  a  zonal  current  is  stable  if  £  varies 

monotonically  with  latitude.  If  c  for  such  a  current  is  plotted 

3. 

as  a  function  of  F,  the  curvature  of  the  profile  curve  will  be  of 
one  sign  in  all  latitudes. 

If  the  mean  absolute  angular  speed  of  rotation  along  a  latitude 
circle  is 


fi  =  -X  (2.156) 

a  2ttIT 

then,  as  shown  by  Fjortoft  (1950),  a  zonal  current  is  stable  if  the 

values  of  in  all  latitudes  w^ipre  6£a/d§  <  0  are  larger  than  any 

value  of  Q  in  latitudes  where  >  0.  These  stability  criteria 
~a  09 

are  sufficient,  but  not  necessary. 

Introducing  the  relative  circulation 


c  =  2ttRu 


(2.157) 


ca  =  2ttR0  +  c 


,£a  ~  BF 


o  a_c 

SF  '  „c2  *1? 


(2.158) 


a=  n  + 

'"•'3 


The  time  rate  of  change 


of  kinetic  energy  of  the  mean  flow  is 


,  d_  /"  ^  u  2dF 

“dt  '  dt  J  ~a 


/9a  ^Ca_  dF  >  0  (Stable) 
2rr  Bt 


(2.159) 


Then  if 


a°s .  !3» 

-55-  2”  SR 


Wtosta  «„  is  tb.  northward  eddy  »*  °*  enguler 


W  =  2ttR  u’v’ 
m 


(2.160) 


(2.161) 


then 


(Stable) 


(2.162) 


Measurements  of  the  mean  meridianal  eddy  flux  of  angular  momentum 
indicate  a  pronounced  northward  flux  south  of  the  latitude  of  maximum 
westerlies,  in  agreement  with  the  fact  that  a  northeast-southwest 
tilt  of  troughs  and  ridges,  corresponding  to  a  positive  correlation 
between  u'  and  v’  is  predominant  in  these  latitudes.  North  of  the 
maximum  westerlies,  on  the  other  hand,  only  small  values  of  Wm  are 
observed  which  are  partly  directed  southward.  On  the  average, 
therefore,  d  Km/d  t  is  positive  in  the  atmosphere,  indicating  that 
kinetic  energy  of  the  perturbations  is  continually  being  converted 
into  kinetic  energy  of  the  mean  current.  Perturbations  in  general 
must  take  their  energy  from  sources  other  than  the  kinetic  energy  of 
the  mean  current.  This  does  not  exclude  the  possibility,  however, 
that  some  growing  disturbances  (perhaps  the  fastest  growing  ones) 
develop  their  instability  as  one  of  this  type,  taking  part  of  their 
energy  from  the  kinetic  energy  of  the  mean  zonal  current. 

Rossby  waves .  The  application  of  the  theory  of  horizontal, 
non-divergent  motion  to  large  scale  currents  in  the  atmosphere  is 
due  to  Rossby  and  colleagues  (1939). 


Noting  the  conservation  of  the  vertical  component  of  the  absolute 
vorticity  C_  =  (X.+£)  =  0,  with  the  notation 

a 

O  -  1  51  _  2g  cos$  (2.163) 

p  a  dtp  a 
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the  conservation  may  be  described  by 


H  =  -v  .  ve  -  pv 


(2.164) 


Rossby  derived  a  wave  of  form 


V 


^  .  sin  K(x  -  ct) 


(2.165) 


in  a  non-divergent  flow,  and  constant  g,  to  have  a  velocity  of 
propagation 


C  =  °R  =  Uo  -  h 


(2.166) 


Rossby* s  formula  says  that  the  waves  propagate  westward  relative  to 
the  air  at  a  speed  proportional  to  the  square  of  its  wavelength. 

The  waves  of  wave  number  k 


ks=V?r 


(2.167) 


are  stationary.  Waves  with  wavelengths  greater  than  the  stationary 
wavelength  Lg  =  2u/ks  move  westward  and  shorter  waves  eastward. 

Shearing  instability 

A  linear  current  of  a  homogeneous  and  incompressible  fluid  is 
possibly  unstable  if  the  shear  of  the  velocity  profile  has  a  maximum 
somewhere  within  the  layer.  Instability  of  this  kind  is  particularly 
pronounced  if  the  shear  of  the  current  is  concentrated  in  a  shear 
line,  which  represents  a  discontinuity  in  the  velocity  field. 

Hoiland  (1942)  has  derived  the  shearing  instability  by  adapting  the 
vorticity  equation.  Introducing  a  unit  vector  t  tangential  to  the 
line  of  the  shear  discontinuity  s ,  and  letting  the  sliding  vorticity 
Q  in  a  point  on  s  be  defined  by 


C  =  Cv, 


(2.168) 
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so  that 


ft  =  ^  *  ^rad  £  •  t  (2.169) 

This  may  be  generalized  to  the  case  when  the  two  fluid  layers,  of 
velocity  =  II 2  have  different  densities  p-^,  p2>  and  there 

is  a  uniform  gravity  field  directed  normal  to  the  undisturbed  surface. 
If  the  lower  fluid  is  the  heavier  one,  the  effect  will  act  to  stabi¬ 
lize  the  system,  waves  longer  than  a  certain  critical  wavelength  will 
be  stable,  and  shorter  waves  unstable.  Perturbation  theory  gives, 
in  this  case,  a  formula  for  the  velocity  of  propagation. 


c 


plUl  +  P2U2  , 
p;,  +  P2 


pl  -  p2 
P1  +  p2 


2  tUl  U2 ' 


(pl  +  P2) 


(2.170) 


The  system  described  by  the  above  equation  is  unrealistic  because 
strict  velocity  discontinuities  do  not  exist  in  nature.  If  the 
velocity  discontinuity  is  replaced  by  a  transition  layer  of  continuous 
velocity  variation,  this  will  not  have  noticeable  effect  upon  waves 
of  wavelength  considerably  longer  than  the  depth  of  the  transition 
layer,  but  waves  of  comparable  Or  shorter  dimension  will  not  be 
stable.  In  this  case,  the  amplifying  fcaves  are  confined  to  an  inter¬ 
mediate  band  of  wavelengths,  whose  width  depends  upon  the  differences 
in  velocity  and  density  in  the  upper  and  lower  layer  and  the  depth 
of  the  transition  layer. 

Billow  clouds,  which  are  formed  in  a  layer  of  high  static 
stability  and  strong  vertical  wind  shear,  are  generally  assumed  to 
result  from  instability  of  this  kind. 

Wave  motions  in  zonal  currents 

Considering  a  steady,  linear,  barotropic  current,  Eliassen  and 
Kleinschmidt  (1957)  describe  an  equation  in  third  degree  in  c,  the 
velocity  of  propagation  of  wave  motion  in  which  the  particles  move 
horizontally  under  a  horizontal  pressure  force  which  is  independent 
of  height. 
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t  -  <u°  ■ c)2]  (u°  - c  -  $ = $ 


(2.171) 


wherein  H  =  RT/g,  X  =  2flsin$ .  The  equation  has  two  large  roots, 
approximately  given  by 


C  ±  S gH 


(2.172) 


which  are  predominanrly  sound  waves,  slightly  modified  by  gravity. 

The  third  root  of  2.171,  which  is  much  smaller  is  found  by  neglecting 
2 

(u  -  c)  versus  gH,  yielding 


-  fi_ 


(2.173) 


This  root  agrees  with  the  observed  speed  of  '-raves  in  the  westerlies 
in  order  of  magnitude,  and  is  an  improved,  but  still  approximate 
form  descriptive  of  the  Rossby  waves.  The  difference  is  slight  for 
wavelengths  less  than  6000  km. 

Fjortoft  (1951)  obtained  a  generalization  of  the  Rossby  formula 
to  the  case  of  the  baroclinic  current.  This  simple  but  crude  theory 
introduced  the  advective  assumption  that  entropy  changes  result  from 
horizontal  advection  only.  Using  the  hydrostatic  and  geostrophic 
equations,  Fjortoft  derived  the  formula  for  velocity  of  propagation 


C  =  -  fc6* 


f^r  -  (a2  -  a?) 


(2.174) 


The  radicand  vanishes  for  waves  whose  wavelength  is 


Lx  =  2tt 


fT{* 


(2.175) 


UXtoA JusSt  *•*•**•«  * 


This  quantity  may  be  calculated  if  uQ  is  known  at  all  levels. 

Clearly,  increases  with  increasing  wind  shear.  For  wind  profiles 
observed  in  the  westerlies,  is  the  order  of  several  thousand 
kilometers,  even  for  moderate  wind  shear. 

A  mere  complicated  solution,  for  the  two  layer  baroclinic 
atmosphere  middle,  has  been  worked  out  by  P.  D.  Thompson  (1961). 

According  to  Charney  (1947)  amplif ic~ tion  will  take  place  pro¬ 
vided  the  wavelength  is  shorter  than  the  critical  value  L^.  According 
to  Eady  and  Fjortoft  sufficiently  short  waves  are  again  stable,  so 
that  amplification  is  confined  to  an  intermediate  band  of  wavelengths. 
The  growth  rate  for  a  given  wind  shear  has  a  maximum  for  a  certain 
wavelength  i,  x  of  the  order  of  3000  km,  which  depends  upon  both 
latitude  and  static  stability. 


max  20sin$ 


(2.176) 


When  the  polar  front  had  been  discovered  by  the  Norwegian 
school  about  1920,  V.  Bjerknes  and  his  collaborators  put  forward 
the  hypothesis  that  the  extratropical  cyclones  are  amplifying  wave 
motions  which  develop  spontaneously  on  the  frontal  surface  as  a 
result  of  some  kind  of  instability  connected  with  it.  H.  Solberg 
(1928)  investigated  theoretically  the  possible  wave  types  in  a 
system  consisting  of  two  statically  stable,  barotropic  layers  of 
different  density  moving  zonally  at  different  speeds  and  separated 
by  a  sloping  interface.  Solberg  found  that  two  types  of  amplifying 
waves  are  possible  ir.  this  system.  The  first  type  is  of  very  short 
wavelength,  up  to  a  few  kilometers,  and  has  clearly  nothing  to  do 
with  cyclones.  The  second  type  is  found  in  a  band  of  wavelengths 
of  the  order  100G  to  2000  km  which  is  in  agreement  with  the  dimensions 
of  cyclones.  Since  the  latter  type  have  certain  kinematical  features 
in  common  with  growing  cyclones ,  it  represents  one  possible  explana¬ 
tion  for  the  growth  of  cyclones.  One  explanation  for  the  two  wave 
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types  is  the  following.  The  growth  of  very  short  unstable  waves 
is  due  to  the  shearing  instability  connected  with  the  wind  discon¬ 
tinuity  at  the  frontal  surface.  At  sufficiently  short  wavelengths, 
the  shearing  instability  dominates  over  the  static  stability,  so 
that  the  disturbances  amplify.  However,  as  the  wavelength  increases, 
the  static  stability  becomes  relatively  more  important,  and  for 
wavelengths  greater  than  a  few  kilometers,  the  waves  are  stable. 

For  all  the  shorter  wavelengths,  the  earth’s  rotation  has  no  notice¬ 
able  effect;  the  fluid  particles  move  in  orbits  in  planes  nearly 
vertical.  As  the  wavelength  is  further  increased,  however,  the 
earth’s  rotation  has  the  effect  of  turning  the  orbital  planes  towards 
the  horizontal,  so  that  the  static  stability  becomes  less  effective. 
From  a  certain  wavelength  on,  the  role  of  the  static  stability  has 
been  thus  reduced  so  that  the  shearing  instability  again  dominates, 
and  the  waves  become  instable.  These  longer  instable  wavelengths 
are  the  cyclone-like  long  waves ,  with  quasi-horizontal  particle 
orbits.  For  still  longer  wavelengths,  the  inertial  stability  due  to 
the  earth’ s  rotation  comes  to  dominate,  and  the  waves  are  again 
stable . 

Solberg’s  long  unstable  waves  take  their  energy  from  the  kinetic 
energy  of  the  mean  zonal  motion,  whereas  the  potential  energy  of  the 
system  is  increased  during  their  growth.  Since  it  is  known  that  the 
perturbations  in  the  westerlies  on  the  average  lose  energy  to  the 
kinetic  energy  of  the  mean  zonal  flow,  it  does  not  seem  possible 
that  Solberg’s  type  of  instability  is  responsible  for  the  growth 
of  all  kinds  of  disturbances  in  the  westerlies.  Certain  types  of 
frontal  waves  of  relatively  short  wavelengths  may  develop  in  accord¬ 
ance  with  Solberg’s  theory. 
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GENERAL  CIRCULATION 


The  conception,  "general  circulation",  comprises  the  largest 
wind  systems  in  the  spectrum  of  atmospheric  motions  down  to  systems 
of  the  magnitude  of  the  cyclones  and  anticyclones.  Averaging  over 
sufficiently  long  periods  of  time  does  away  with  the  smaller  among 
these  systems.  The  remaining  ones  are  large  streams  encircling 
the  entire  globe  which  remain  nearly  unchanged  when  annual  averages 
are  taken.  These  are  in  the  first  instance  implied  by  the  term 
general  circulation.  Certain  irregularities  are  caused  by  land  and 
sea  distribution  and  by  mountains.  The  following  treatment  is  largely 
taken  from  reviews  by  Lorenz  (1967),  Starr  (19S8),  and  Mintz  (1961). 

The  motion  in  these  currents,  which  have  a  crude  symmetry  with 
respect  to  the  earth's  axis,  is  predominantly  zonally  directed,  but 
has  also  a  meridional  component  that  transports  energy  poleward  from 
the  equator. 

The  principal  problem  pertaining  to  these  streams  is  the  ques¬ 
tion  of  their  maintenance.  Friction  opposes  their  continuation, 
presumably  seeking  to  reduce  all  velocities  to  that  of  a  solid 
rotation  with  the  earth,  that  is  to  produce  a  general  calm.  To  main¬ 
tain  the  general  circulation  against  friction,  heat  sources  and  sinks 
are  necessary.  The  theory  of  the  general  circulation  thus  differs 
from  much  of  the  other  theoretical  work  that  has  been  performed  in 
that  both  friction  and  heat  sources  and  sinks  are  drawn  into  considera¬ 
tion.  Friction  and  heating  energetics  lead  to  flow,  which  is  taken  as 
an  assumed  parameter  by  most  mathematical  treatments. 

The  eddy  motions  evolving  from  the  predominantly  zonal  currents 
comprise  the  convective  cells  for  the  conversion  of  essentially  heat 
energy  into  kinetic  energy  in  a  rotating  atmosphere.  There  are  two 
modes  of  response  of  a  convective  model:  symmetric  and  cellular. 

The  symmetric  mode,  which  occurs  when  the  rate  of  rotation  is  small, 
is  meridionally  similar  to  the  field  of  motion  induced  in  a  non¬ 
rotating  stratum  of  gas  by  a  constant  and  symmetric  distribution  of 
heating  at  low  and  cooling  at  high  latitudes.  Such  non-rota tional 
motion  is  devoid  of  zonal  components  and  takes  place  in  meridianal 
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planes,  being  identical  for  all  longitudes.  The  symmetric  mode  is  that 
field  of  motion  induced  in  a  curved  stratum  of  gas  having  a  relatively 
slow  rate  of  basic  rotation  about  the  polar  axis  and  heated  by  a  con¬ 
stant  and  symmetric  distribution  of  heating  at  low  latitudes  and  cool¬ 
ing  at  high  latitudes.  It  is  characterized  in  the  meridianal  plane  as 
a  single  cell  with  air  rising  near  the  equator  and  sinking  near  the 
pole.  As  compared  with  the  non-rotating  case,  relative  westerlies 
appear  at  higher  and  relative  easterlies  at  lower  latitudes  with  a 
sloping  surface  of  separation  between  the  two.  Ihe  flow  is  still  the 
same  at  all  longitudes,  thus  being  symmetric  about  the  polar  axis. 

The  cellular  mode,  called  also  a  high-rotation  regime  or  Rossby 
regime,  is  developed  by  passing  from  the  low-rotation  regime  by  keeping 
other  quantities  the  same  but  increasing  the  rotation  a  moderate  amount. 
Such  motions  are  turbulent,  forming  in  the  meridianal  plane  three  cells. 
Of  these,  the  cell  in  middle  latitudes  has  a  mean  circulation,  in  the 
meridianal  plane,  which  is  contrary  to  the  local  thermal  forcing,  in¬ 
volving  a  descent  of  relatively  warm  gas  and  an  ascent  of  colder  gas. 
This  forced  effect  is  a  result  of  the  turbulent  character  of  the 
motions.  The  relative  mean  flow  of  the  cellular  mode  also  shows 
westerlies  at  higher  latitudes  and  easterlies  at  lower  latitudes. 

The  convection  in  the  turbulent  regime  is  difficult  to  picture 
simply  because  of  the  random  components  of  motion  associated  with  it. 
However,  convection  requires  a  net  systematic  statistical  effect,  in 
the  sense  of  cold  gas  sinking  and  warm  gas  rising. 

The  radiative  cooling  at  high  latitudes  and  heating  at  low  lati¬ 
tudes  creates  a  pool  of  dense  and  a  pool  of  light  gas,  creating  a 
zonally  distributed  store  of  available  internal  and  potential  energy 
for  conversion  into  kinetic.  The  effect  of  the  turbulence  is  to  break 
off  smaller  masses  of  cold  gas  at  individual  longitudes  here  and  there, 
which  then  migrate  to  middle  and  lower  latitudes  and  do  their  sinking 
sporadically  in  these  places.  Also,  the  warm  gas  is  first  moved  in 
individual  poleward-moving  streamers  toward  middle  and  higher  latitudes 
over  narrow  longitude  sectors,  to  be  then  displaced  upward  by  the  sink¬ 
ing  and  horizontally  spreading  cold  gas  blobs.  In  the  turbulent 
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process  are  four  steps  of  the  energy  cycle  of  the  earth Ts  atmosphere: 

(1)  the  generation  of  zonal  available  potential  plus  internal  energy 
by  differential  heating  through  radiation;  (2)  the  transformation  of 
the  zonal  available  into  the  eddy  available  form  of  this  energy;  (3)  the 
transformation  of  eddy  available  potential  plus  internal  energy  into 
eddy  kinetic  energy,  evidenced  by  a  pattern  along  a  single  zone  of  warm 
air  in  general  rising  and  cold  air  sinking;  and  (4)  a  transport  of 
momentum  either  in  the  direction  of  or  contrary  to  the  gradient  of  mean 
flow.  Required  for  the  latter  transport,  called  an  effect  of  negative 
viscosity,  is  a  supply  of  eddy  kinetic  energy,  drawn  from  the  eddies  by 
conversion  of  other  forms  of  energy  (potential,  internal,  etc.)  within 
the  eddies  or  transferred  from  another  scale  of  eddies,  including 
smaller  ones  in  which  case  the  smaller  ones  must  receive  a  supply  from 
other  sources.  In  the  case  of  a  system  in  which  the  eddies  are 
thermally  driven,  the  heat  dissipated  by  braking  of  the  mean  flow  by 
friction  or  the  heat  produced  by  hydroma gnetic  action,  by  the  second 
law  of  thermodynamics,  cannot  again  be  used  to  drive  the  eddies. 

Following  the  notation  of  Starr  (1968),  u,  v,  and  w  are  the  zonal, 
meridional,  and  vertical  components  of  wind  velocity,  R,  $,  and  6  are 
the  radius,  latitude,  and  longitude  coordinates  of  the  earth,  and  p 
and  p  are  the  pressure  and  density  of  the  atmosphere.  Also,  a  bar  (— ) 
denotes  the  time  average  of  any  quantity,  a  prime  (’)  the  departure  of 
a  quantity  from  its  time  average,  brackets  ([])  the  average  of  any 
quantity  with  respect  to  longitude,  and  a  star  (*)  the  departure  of  a 
quantity  from  its  longitudinal  average.  The  operators  —  ,  ’,  [] ,  and 
*  are  commutative. 

The  wind  field  vector  components  may  be  resolved  with  such  nota¬ 
tion  according  to  the  formulae  for  U 


u  =  u  +  u’ 

(2.177) 

U  =  tu]  +  u* 

(2.178) 

thus,  in  greater  detail. 

U  =  [U]  +  U*  +  [U]»  +  U*T 

(2.179) 
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Resolutions  of  the  fields  of  temperature  T,  specific  humidity  q,  and 
other  quantities  may  be  similarly  performed.  Since  his  data  consist 
mainly  of  observations  at  standard  pressure  levels,  Starr  lets  U  repre¬ 
sent  the  time  average  for  fixed  X,  cp,  ?  (not  the  same  as  for  fixed  X, 
cp,  z).  Similarly,  [U]  is  the  average  with  respect  to  longitude  along 
a  latitude  circle  on  an  instantaneous  isobaric  surface. 

To  a  considerable  extent  the  features  appearing  respectively  in 
the  fields  of  [U],  11*,  [U] ’  and  U*f  are  those  appearing  in  the  field 
of  [U] ,  in  U  but  not  [U],  in  [U]  but  not  [Uj,  and  in  U  but  neither  U’ 
nor  [U].  Similar  remarks  apply  to  the  features  of  T  and  q.  However, 
some  meteorological  features,  such  as  the  jet  stream,  do  not  clearly 
fall  into  any  one  category. 

The  fields  of  11  and  UT  in  (2.177)  are  called  the  long-term  or 
time-averaged  or  standing  motion  and  the  transient  motion.  The  compo¬ 
nents  [U]  and  TV]  are  called  the  zonal  circulation  and  the  meridional 
circulation,  and  the  components  of  U*  and  V*  are  called  eddies.  The 
field  of  [W]  demanded  by  continuity  is  also  included  as  part  of  the 
meridional  circulation.  Thus  the  terms  in  (2.179)  are  respectively 
the  time-averaged  or  standing  zonal  and  meridional  circulations,  the 
time-averaged  or  standing  eddies,  the  transient  zonal  and  meridional 
circulations,  and  the  transient  eddies. 

To  avoid  ambiquity  in  use  of  the  terms  "zonal"  and  "meridional” 
in  this  report,  a  "zone"  means  a  latitude  circule  or  a  region  extend¬ 
ing  along  a  latitude  circle,  "zonal  motion"  means  motion  parallel  to 
the  zones,  and  is  synonymous  with  U,  "meridional  motion"  means  motions 
parallel  to  the  meridians  or  meridional  planes.  A  "zonal  average" 
means  an  average  within  zones,  or  with  respect  to  longitude,  and 
"zonal  symmetry"  denotes  invariability  within  zones. 

Figure  29  shows  the  cellular  meridional  flow  Suggested  by  Rossby 
(1945).  Although  some  features  of  his  model  are  no  longer  acceptable, 
Rossby *s  picture  relates  the  meridional  cells  to  the  wavy  bounds  of 
the  polar  front  zone  and  the  horse  latitudes,  which  separate  the  polar 
easterlies  and  the  mid- latitude  westerlies,  in  one  case,  and  the  mid¬ 
latitude  westerlies  and  the  equatorial  easterlies  in  the  other. 
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FIGURE  29. 


Cellular  Meridional  Circulation  on  a  Rotating  Earth 
(After  G.C.  Ross  by,  1945) 


The  vertical  distribution  of  annual  mean  zonal  velocities  in  the 
atmosphere,  averaged  for  all  longitudes  [U],  adapted  from  observational 
material  compiled  by  Buch  (1954)  for  the  Northern  Hemisphere  and  by 
Obasi  (1965)  for  the  Southern  Hemisphere  is  shown  in  Fig.  30.  More 
recent  compilations,  such  as  that  by  Starr  (1969)  given  in  Figs.  3 land 
32,  suggest  that  Buch’s  maximum  westerlies  are  slightly  too  low. 

The  vertical  distribution  of  annual  mean  meridional  velocities  in 
the  atmosphere  averaged  for  all  longitudes  [V]  is  shown  in  Fig.  33. 
These  data  given  by  Starr  (1968)  were  adapted  from  Lorenz  (1967)  who 
in  turn  utilized  observational  material  from  Buch  (1954)  for  the 
Northern  Hemisphere  and  by  Obasi  (1963)  for  the  Southern  Hemisphere. 

The  horizontal  arrows  in  Fig.  33  were  obtained  by  calculation  from 
measurements  of  momentum  transports  and  other  observed  data.  The 
vertical  arrows  are  only  schematic  and  are  not  to  scale.  A  more  recent 
compilation  of  Northern  Hemisphere  data  is  given  in  Figs.  31  and  32. 

The  vertical  distribution  of  the  mean  annual  eddy  transport  of 
o  2  _  _ , 

angular  momentum  2nR  cos  $  {[UTVT]  +  in  the  actual  atmosphere 

is  given  in  Fig. -34.  The  data  of  Fig.  34  were  adapted  by  Starr  (1968) 
after  Lorenz  (1967),  who  used  data  compilations  of  Buch  (1954)  for  the 
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Northern  and  of  Obasi  (1963)  for  the  Southern  Hemisphere.  Both  the 
effect  of  transient  eddies  and  of  the  standing  eddies  are  included. 


The  annual  average  curves  derived  from  observations  in  the  actual 
atmosphere  giving  the  relative  mean  angular  velocity  about  the  polar 
axis  [U]/R  cos$  (dotted  line),  the  northward  angular  momentum  transport 

O  r\ 

by  turbulent -eddies  2ttR  cos  |/  pCU^V’  ]dz  (smooth  line)  and  the  angular 

o  o  _  __ 

momentum  transport  by  mean  meridianal  circulations  2ttR  cos~§f  p[U][V]dz 
(dashed  line)  are  shown  in  Fig.  35.  •  The  two  transports  are  integrated 
vertically  through  almost  the  entire  mass  of  the  atmosphere  (up  to  20 
Km).  The  angular  velocity  represents  a  vertical  average  through  the 
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FIGURE  31 .  Long-Term  Averages  of  Wind  Velocities,  Mass  Flow  and  Momentum 
Parameters  (After  V.P.  Starr  et  a!.,  to  be  Published,  1969) 
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same  mass.  All  curves  are  by  Starr  (1968)  who  adapted  them  from  the 
review  by  Lorenz  (1967)  who  in  turn  used  data  analyses  by  Buch  (1954) 
for  che  calendar  year  1950  for  the  Northern  Hemisphere  and  by  Obasi 
(1965)  for  the  calendar  year  1958  for  the  Southern  Hemisphere.  More 
recent  data  by  others,  including  the  data  by  Starr  (1969)  given  in 
Figs.  31  and  32,  corroborate  in  a  general  way  these  data.  The  dif¬ 
ference  between  the  hemispheres,  shown  in  Fig.  35,  is  real.  Also,  a 
negative  eddy  viscosity  effect  is  illustrated  by  the  data  that  have 
largest  numerical  values  of  eddy  transport  located  in  the  zones  of 
greatest  shear  in  angular  velocity. 


Equator 


FIGURE  35.  Angular  Velocity  Averaged  over  Vertical  and  Angular  Momentum 
Transports  Integrated  over  Vertical  (After  Starr,  1968) 
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Figure  31  shows  data,  recently  compiled  by  Starr  (1969)  from  re¬ 
ports  of  about  800  stations  over  the  five-year  period  starting  Ffey  1, 
1958,  which  describe  the  long-term  (60-month)  time-  and  zonally- 
averages  in  the  Northern  Hemisphere  of  zonal  wind  [U],  meridional  wind 
[V]  and  stream  function  The  stream  function  measures  mass  flow. 
Between  adjacent  lines  of  the  stream  function,  mass  transport  is  con¬ 
stant.  The  mass  flow  is  dependent  on  cos§,  since 

1  =  2riR  cosf  [V] 

up 

q  ||  =  -2t®2  cosf  [W] 

In  Fig.  32  are  shown  the  seasonal  averages  of  the  same  data  used  for 
Fig.  31. 

Also  shown  in  Figs.  31  and  32  are  the  momentum  parameters: 

[U * V r ]  cos^$  wiich  is  proportional  to  the  northward  angular  momentum 

transport  by  transient  eddies,  [U]/cosi  which  is  proportional  to  the 

~  _  2 

relative  angular  velocity  and  [U*V*]cos  §  which  is  proportional  to  the 
northward  angular  momentum  transport  by  standing  eddies. 

Figure  36  displays  vertical  meridional  cross  sections  through  the 
atmosphere  showing  (a)  the  distribution  of  the  mean  temperature  [T], 

(b)  the  distribution  of  northward  heat  transport  by  transient  eddies 
to  be  obtained  from  the  values  of  [V’T* ]cos$  by  multiplying  the  values 
by  2nCpR,  and  (c)  the  distribution  of  northward  heat  transport  by 
standing  eddies  to  be  obtained  from  the  values  of  [V*T*]cos§  by  multi¬ 
plying  one  values  by  2nCpR.  These  data  were  obtained  by  Starr  (1968) 
by  automatic  computing  from  five  years  of  daily  data  at  ten  levels  for 
about  700  stations  in  the  Northern  Hemisphere. 

Figure  37  shows  the  transport  of  water  vapor  across  latitude 
circles  by  the  motions  of  the  atmosphere  on  a  mean  annual  basis  for 
the  calendar  year  1958,  adapted  by  Starr  (1968)  from  transport  calcula¬ 
tions,  by  Peixoto  and  Crisi  (1965)  over  the  Northern  Hemisphere,  and 
by  Starr,  Peixoto,  and  McKean  (1968)  over  the  Southern  Hemisphere. 
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FiGURE 36.  Temperature  and  Transpor'  of  Heat  (After  Starr,  1968) 


104 


r 


■JTOL 


FIGURE  37.  Transport  of  Water  Vapor  (After  Starr,  1968) 


Just  as  there  exists  a  transfer  of  kinetic  energy  from  the  eddies 
to  the  0  harmonic  (i.e. ,  the  mean  zonal  kinetic  energy),  a  redistribu¬ 
tion  of  kinetic  energy  takes  place  among  the  wave  numbers  from  1  to  15 
and  probably  beyond.  Figure  38  shows  the  results  of  a  spectral  study 
based  upon  daily  hemispheric  wind  data  for  a  period  of  nine  years  at 
the  5  km  level  (near  500  mb)..  The  main  loss  to  other  wave  numbers  is 
from  the  central  band,  from  wave  numbers  5  to  11  inclusive,  correspond¬ 
ing  to  wavelengths  at  45  deg  latitude  of  about  2000  to  about  5000  km. 

This  energy  is  given  up  to  both  lower  and  higher  wave  numbers,  with  the 
loss  being  made  good  again  by  a  transformation  of  potential  and  heat 
energy  into  kinetic  in  this  central  band.  At  the  right  in  the  diagram 
is  shown  the  loss  of  kinetic  energy  from  each  wave  number  to  the  zonal 
flow.  All  15  wave  numbers  feed  energy  to  the  zonal  flow.  The 
especially  large  loss  from  the  box  for  wave  number  2  probably  reflects 
a  monsoonal  action. 

* 
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UNITS  (  10  ^ergs  cm  ^sec  'mb  '  ) 

FIGURE 38.  Spectrum  of  Kinetic  Energy 
Exchange  at  5-km  Altitude 
(After  Starr,  1968, 

Adapted  from  Saltzman 
and  Teweles,  1964) 


Figure  39  is  after  Oort 
(1964)  and  Starr  (1968)  and  is  a 
schematic  diagram  illustrating 
the  balance  of  four  types  of 
energy  for  the  annual  average 
condition  of  the  entire  atmos¬ 
phere.  The  parentheses  enclose 
Starr's  estimates. 

The  arrow  C  indicates  a  flow 
of  energy  from  the  eddies  to  the 
mean  flow.  The  heating  of  the 
atmosphere  indicated  at  the  top 
of  the  diagram  is  due  to  such 
sources  as  heating  by  condensa¬ 
tion  of  water  vapor  and  by 
radiation. 

Table  8  shows  five-year  and 
seasonal  averages  for  the  trans¬ 
formations  of  energy  in  the  at¬ 
mosphere,  derived  by  Starr  (1968). 
His  estimates  are  derived  from 
available  data  for  the  Northern 


Pe  r  m  *  Hemisphere  atmosphere  for  the 

and  Teweles,  1964) 

period  May  1,  1958,  to  April  30, 

1963.  Because  of  nonlinear  effects,  the  figures  for  the  total  five- 
year  period  are  not  equal  to  the  means  for  the  four  individual  seasons. 

Wave  Regime  of  Thermal  Circulation 

Mintz  (1961)  has  noted  that  there  are  two  regimes  of  large-scale 
thermal  circulation  that  can  develop  in  a  thin  layer  of  fluid  which  is 
gravitationally  held  to  a  rotating  sphere  and  is  heated  at  the  equator 
of  rotation  and  cooled  at  the  poles .  In  each  of  these  general  circula¬ 
tion  regimes  the  fluid  performs  two  basic  functions  wnich  can  hold  for 
a  rotating  sphere :  it  maintains  thermal  equilibrium  by  transporting 
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FIGURE  39.  Energy  Balance  and  Transport  (After  Oort,  1964, 
and  Starr,  1968) 
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TABLE  8.  TRANSFORMATIONS  OF  ENERCY  IN  THE  ATMOSPHERE, 
on 

UNITS  OF  10  ergs  sec  (After  Starr,  1968) 


KE  signifies  kinetic  energy  and  P+I  is  potential 
plus  internal  energy. 


heat  from  the  equatorial  heat  source  to  the  polar  cold  source,  and  it 
produces  the  kinetic  energy  necessary  to  maintain  the  circulation 
against  frictional,  dissipation.  The  earth's  atmosphere  is  characterized 
by  the  v,ave  region.  Discussion  of  the  symmetric  region  is  given  to 
illustrate  why  the  earth  is  in  the  wave  regime. 


In  the  symmetric  regime  of  general  circulation,  the  fluid  in 
all  longitudes,  ascends  near  the  equator  of  the  rotating  sphere, 
flows  poleward  in  the  upper  levels ,  descends  near  the  pole  and 
flows  equatorward,  increasing  in  intensity  with  height  in  the 
direction  of  rotation  of  the  sphere.  In  this  regime  the  meridional 
and  zonal  components  of  the  motion  are  the  same  in  all  longitudes. 

The  circulation  is  symmetric  about  the  pole,  as  shown  in  Fig.  40. 

In  the  other  form  of  general  circulation,  called  the  "wave 
regime,"  the  flow  is  characterized  by  large  amplitude  horizontal  waves 
(baroclinic  Rossby  waves  described  earlier  by  Eq.  2.174)  in  the  middle 
and  upper  levels,  and  by  large  horizontal  eddies  (cyclones  and  anti¬ 
cyclones)  in  the  lower  levels  of  the  fluid,  as  shown  schematically  in 
Fig.  41.  In  each  hemisphere,  the  centers  of  the  cyclones,  in  the  mean, 
lie  c3oser  to  the  pole  than  the  centers  of  tne  anticyclones  and,  as  a 
result,  the  relative  zonal  flow  at  low  levels,  when  averaged  over  all 
longitudes,  is  easterly  near  the  equator  and  easterly  near  the  pole, 
but  westerly  in  middle  latitudes.  In  addition,  in  the  wave  t'egime  the 
zonaily  averaged  meridional  circulation  is  reversed  in  middle  latitudes, 
as  shown  in  the  figure. 

In  the  wave  regime  of  circulation,  there  are  the  tongues  of  warm 
fluid  moving  poleward  in  some  longitudes ,  while  at  the  same 
elevation  tongues  of  colder  fluid  move  equatorward  in  other  longi¬ 
tudes,  which  produce  the  poleward  heat  transport  that  balances 
the  differential  latitudinal  heating.  Unlike  the  wave  regime, 
in  the  symmetric  regime  a  net  poleward  heat  transport  occurs  only 
when  the  vertical  lapse  rate  of  temperature  is  less  than  the 
adiabatic  rate,  so  that  the  northward  moving  fluid  of  Fig,  41. 
has  a  higher  potential  temperature  than  the  southward  moving  fluid. 

Critical  poleward  temperature  gradient 

Mintz  (1961)  has  described  the  critical  poleward  temperature 
grad '  nt  for  dynamic  stability  of  each  wave  number  in  the  symmetric 
ren-n-i.  Change  in  poleward  temperature  gradient  may  affect  the 
d  -tv.,  ic  stability.  For  temperature  gradient  smaller  than  the  critical 
temperature  gradient,  the  symmetric  regime  is  stable;  for  larger 
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STREAMLINES  OF  THE  FLOW  AT  THE  UPPER  LEVELS. 

:ross  section  showing  the  meridional  projection  of  the  circulation  and 
ISOTHERMS  OF  POTENTIAL  TEMPERATURE  6 . 


FIGURE  40.  Symmetric  Regime  of  General  Circulation 
(After  Mintz,  1961) 


PRESSURE 


POLE  — - LATITUDE  EQ. 


A)  STREAMLINES  OF  THE  FLOW  AT  THE  MIDDLE  AND  UPPER  LEVELS  (HEAVY 
LINE)  AND  NEAR  THE  GROUND  (THIN  LINE).  L  =  LOW  PRESSURE, 

H  =  HIGH  PRESSURE. 

B)  CROSS  SECTION  SHOWING  THE  ZON ALLY-AVERAGED  MERIDIONAL 
CIRCULATION;  AND  THE  ZON  ALLY-AVERAGED  ZONAL  WIND,  WHERE 
W  =  WESTERLY  WIND,  E  =  EASTERLY  WIND. 


FIGURE  41.  A  Wave  Regime  of  General  Circulation  (After  Mi ntz,  1961) 
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gradients,  the  symmetric  regime  is  unstable  and  waves  develop.  In  his 
development  Mintz  uses  the  two-level  baroclinic  model,  applying  the 
quasi-geostrophic  vorticity  equation 


__2  =  _v  •  v(C  +f)-fv*^+k*yxF  (2.180) 
3t  ~g  v  g  ' 

to  levels  1  and  1,  taken  as  the  250  and  750  mb  levels,  and  the 
thermodynamic  equation 


ae  = 
at 


-v 


ve 


9 

CPT 


q 


(2.181) 


to  level  2,  taken  at  the  500  mb  level.  Where  £  =  k  •  V  x  v 

is  the  vertical  component  of  the  curl  of  the  geostrophic  wind, 
w  =  ap/dt  is  the  vertical  velocity  in  x,  y,  p- space,  e  is  the 

o 

potential  temperature  and  q  is  the  rate  of  non- adiabatic  heating 
per  unit  mass.  Mintz  makes  further  use  of  the  geostrophic  wind 
equation 

Xg  =  |kxvz1  (2.182) 

the  equation  of  continuity  of  mass 


7  *  y  =  -dw/ap 


(2.183) 


Poisson's  equation  for  potential  temperature 


R 


the  equation  of  state 


and  the  hydrostatic  equation 


pg 


= 

szi 


(2.184) 


(2.185) 


(2.186) 
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In  the  above  expressions,  f  =  2  Q  sin  §  is  the  Coriolis  force, 

Q  is  the  rate  of  rotation  of  the  planet,  §  is  the  latitude,  F  is 
the  horizontal  friction  force ,  V  is  the  horizontal  del  operator  on 

A 

the  isobaric  surface,  k  is  the  unit  vertical  vector,  t  is  the  time, 
z  is  the  height  of  the  isobaric  surface,  v  is  the  horizontal 


velocity,  g  is  the  acceleration  of  gravity,  p  is  the  vertical  eddy 
viscosity,  6  is  the  density,  s  =  (1  -  is  the  static  stability. 


d  a 

emperature  gradient ,  =  -  •£—  is  the 


Y  =  dT/3z  is  the  vertical  temperature  gradient,  =  ~  is  the 
adiabatic  vertical  temperature  gradient,  R  is  the  universal  gas 


constant,  m  is  the  mean  molecular  weight,  is  the  coefficient  of 
specific  heat  at  constant  pressure,  a  is  the  radius  of  the  earth, 


l ft  is  the  net  heating  of  the  atmosphere  on  the  equatorward  side  of 


the  central  latitude . 


Using  the  relations  developed  by  P.  D.  Thompson  (1961)  for 
the  two  level  baroclinic  atmosphere,  and  linearizing  by  neglecting 
the  heating  arid  friction  terms  and  by  requiring  that  the  zonal 
component  of  the  geostrophic  wind  be  horizontally  uniform  and 
independent  of  time,  Mintz  describes  the  behavior  of  any  sinusoidal 
perturbation  of  zonal  flow  as  being  determined  by  the  parameter 


where 


(~^2\  _  a2n2ml2 
\  95  /  Rn2(X4-n4)^ 


(2.187) 


(2.188) 


and  n  is  the  planetary  wave  nu Tiber ,  i.e.,  the  number  of  waves  around 
the  hemisphere.  When  5  <  0,  the  perturbation  at  both  levels  oscillates 
about  its  initial  amplitude.  When  6  >  0,  the  perturbation  grows  in 
amplitude  without  limit- -that  is,  the  waves  amplify  and  the  wave 
perturbation  acts  to  reduce  the  poleward  temperature  gradient. 
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-3T ’^b(p),  IN  UNITS  OF  KT 


Critical  temperature  gradient  and  doming r>t  wave  number 

A  plot  of  the  neutral  curve  6  =  0  as  a  function  of  poleward 
temperature  gradient  (-  dTg/d  $)  and  planetary  wave  number  n  is  shown 
in  Fig.  42.  As  may  be  seen  in  the  figure  there  is  a  critical  value 
of  the  poleward  temperature  gradient. 


FIGURE  42.  Stability  of  the  Symmetric  Circulation  (After  Mintz,  1961) 


Mi 


(2.189) 


below  which  6  <  0  for  perturbations  of  all  wavelengths.  When 

(-  aTo/d0)  barely  exceeds  (R/c  )T0s,  perturbations  of  planetary 
2  u  p  2 

wave  number  n  =  X/2  have  6  >  0  while  longer  and  shorter  waves 
still  have  6  <  0.  This  particular  value  of  wave  number  is  a 
critical  one,  and  is  termed  the  "dominant  wave  number": 


Should  the  poleward  temperature  gradient  be  larger  than  the  critical 
value  (R/cs)T2s,  then  6  is  positive  over  a  broader  range  of  wave¬ 
lengths.  However,  it  is  never  positive  for  wave  numbers  larger 
than  X. 

Velocity  of  wave  propagation 

The  dominant  wave  number  nd  is  important  because  the  critical 
value  of  ( -QT2/a0)  never  can  be  very  much  exceeded  except  when 
X  £  1.  We  may  imagine  a  differential  heating  AQ  which  exceeds  the 
heat  transport  characteristic  of  the  symmetric  regime  and  thereby 
makes  the  poleward  temperature  gradient  ( -t^T2/A(?)  gradually  increase 
to  approach  its  critical  value  from  below.  Because  of  irregularities 
of  the  surface  of  the  atmosphere ,  small  amplitude  perturbations  in 
all  wavelengths  will  always  be  present.  As(-$T2/B0)  passes  its 
critical  value,  perturbations  whose  wavelengths  lie  near  the  dominant 
wavelength  will  start  amplifying,  and  transport  heat  polewards. 
Thermal  equilibrium  will  first  be  reached  when  the  waves  have  grown 
large  enough  in  phase  lag  and  amplitude  for  the  poleward  transport 
to  balance  the  heating.  As  the  waves  continue  to  grow,  the  heat 
transport  becomes  larger  than  the  AQ  heat  input,  (-8T2/dQ)  in  time 
is  reduced  to  less  than  its  critical  value,  giving  rise  to  an 
oscillation  in  which  thermal  equilibrium  is  maintained  by  the 
waves  only  in  the  long  term  average. 
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wave 


The  criterion  for  the  establishment  and  maintenance  of  the 
regime,  is  whether  (-  AT ^ /d0)H  _aq  is  larger  or  smaller 

than  (-  dT2^a^crit  wherein  (-  ST2/d0)H  is  the  poleward 
temperature  gradient  sufficient  to  tranlport  the  net  heating 
of  the  atmosphere  on  the  equatorward  side  of  the  central  latitude 
under  steady  state  symmetric  circulation 


and  the  critical  value  of  net  heating  is 


That  is,  when 


ACWt 


2~  2 
ngR  T^s 

VTmc  n2 

P 


(2.191) 


(2.192) 


(2.193) 


and 


AQ  > 


A!W 


the  wave  regime  of  general  circulation  will  form. 

In  the  long  time  average,  the  heat  transport  equals  the  net 
heating  equatorward  of  the  central  latitude.  The  total  poleward 
heat  transport  across  the  central  latitude  by  the  waves  when  the 
surface  winds  are  weak  compared  to  the  upper  levels  is  compared  to 
the  upper  levels  is 

rra2c  HP  V?  sin  a  , ,  - 

7j  P  S  2 _  \2.l 

%  ~  Rgn 


where  a.  is  the  phase  lag  between  the  waves  at  levels  1  and  3.  Hw 

depends  not  only  on  the  magnitudes  of  and  a  but  also  on  the 
correlation  between  them.  An  order  of  magnitude  estimare  of  the 


wind  speed  average  of  levels  1  and  3  required  for  wave  regime  trans¬ 
port  of  net  heating  AQ  is 


(2.195) 


Wave  dynamic  parameters  of  the  earth 


The  theoretical  determination  cf  the  regime  of  general 
circulation  of  the  earth’s  atmosphere  is  obtained  by  comparing 
AQ  with  AQcrit.  To  a  first  approximation,  Q  for  the  earth  is 
independent  of  the  temperature  distribution.  The  net  heating  per 
unit  horizontal  area,  as  a  function  of  latitude  is  (S(l  -  A)  -  W), 
where  3  is  the  incident  solar  radiation  derived  from  astronomical 
relations,  A  is  the  planetary  albedo,  taken  as  a  constant  A  =  0.34, 
and  W  is  the  long- wave  radiation  emitted  by  the  planet.  This  yields 
the  mean  annual  distribution  of  S(1  -  A)  shown  in  Fig.  43. 


Because  the  earth  has  extensive  oceans  that  are  heated  by  the 
sun,  there  is  widespread  and  persistent  upward  diffusion  of  water 
vapor  from  the  ocean  source,  and  the  air  in  all  latitudes  is  in  a 
state  of  near  saturation.  Mints  assumes  for  his  model  of  the  earth 
that  the  relative  humidity  is  uniformly  high  in  all  latitudes  and 
that  there  is  sufficient  water  vapor  in  the  atmosphere  so  that  all 
the  outward  radiation  may  come  from  the  water  vapor  in  the  air. 

Under  these  conditions,  where  the  temperature  is  high,  as  it  is  near 
the  equator,  the  water  vapor  which  emits  the  outgoing  radiation  will 
be  at  relatively  high  levels  in  the  atmosphere;  where  the  temperature 
is  lew,  as  it  is  in  the  troposphere  of  the  poles,  the  emitting  water 
vapor  will  be  at  lower  levels.  But  the  intensity  of  the  outgoing 
radiation,  which  depends  essentially  only  on  the  temperature  of  the 
water  vapor  and  not  on  its  elevation,  will  be  of  the  same  magnitude 
everywhere.  Hence,  in  Mintz’ model,  the  outgoing  radiation  W  is  taken 
as  constant  with  latitude.  For  thermal  equilibrium,  the  magnitude 
of  VJ  must  equal  the  latitude  mean  of  S(1  -  A)  which  is  known.  Hence, 
there  results  a  net  heating  between  the  equator  and  latitude  37°, 
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(2.196) 


IQ  -  T 

and  the  pole  of  10.6  (10)  cal  day 

AQ  ss  10.6(10)19  cal  day”1 
Further,  Mintz  takes  the  following  values  of  the  earth  parameters: 


a  =  6. 37(10)  ra 

-2 

g  =  9.81  msec 
ft  =  7.29  (10)" 5  sec”1 
T2  =  260°C 

cp  =  1004  KJ  ton-1  deg”1 


R/m  =  287  KJ  ton"1  deg”1 
H  =  2.2(10)~2ton  m  1  si 
Yd  =  -g/c  =  -10°C/Km 


Y  =  Yd/2 


S  =  (1 


-  JU  = 


■)  =  0.5 


to  derive  the  computer  parameters 

AQ  =  5. 1(10) 12  KJ  sec”1 
AQcrit  =  °-6(10)12  ^  sec”1 


=  310°C/rad 


(— 2) 

V  W  4s=aq 

(-S?2)  .  '  38°c'r*d 

\  r  '  en  t 


crit 

Y  =  6.36 


(v2) 

z  H 


nd  =  5.4 
L.  =  5300  Km 


=  13  m  sec 


-1 


w=AQ 


SMS  (V2) 


=  9  m  sec 


-1 


Hw=° 


(2.197) 


(2.198) 
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For  the  Mintz  model  of  the  earth's  atmosphere  the  differential  heating 
AQ  is  many  times  larger  than  AC^  ,  the  maximum  for  a  dynamically 
suable  symmetric  circulation.  Since  with  only  symmetric  circu¬ 
lation  the  poleward  temperature  gradient  required  to  transport  the 
heat  from  the  equator  to  the  poles  would  require  a  gradient  of 
310°C  per  radian,  but  the  critical  dynamical  limit  of  the  poleward 
temperature  gradient  is  only  38°C  per  radian,  planetary  waves  must 
develop  in  the  atmosphere  whose  greater  efficiency  can  transport  the 
net  heating.  The  number  of  waves  around  the  earth  will  be  the 
nearest  whole  number  to  the  dominant  wave  number  n^  =  5.4, 
corresponding  to  the  dominant  wavelength  at  latitude  45  deg  of 
5300  km.  The  root  means  spare  meridional  wind  is  9  m  sec  \ 

The  computations  described  above  are  for  an  atmosphere  in  which 
there  is  no  release  of  latent  heat  of  condensation.  If  condensa¬ 
tion  is  taken  into  account  then  the  saturated  adiabatic  process 
rate  ys  must  replace  the  unsaturated  adiabatic  process  rate  y^ 
over  about  one- half  of  the  wave  domain,  so  that 


(2.199) 


The  effect  of  replacing  s  in  the  computations  by  s*,  after  the 
waves  have  developed  to  the  extent  of  producing  condensation,  has 
two  important  consequences:  it  reduaes  the  critical  poleward  tempera¬ 
ture  gradient  with  only  a  small  change  in  the  dominant  wave  number, 
and  it  produces  a  condition  of  hypercritical  instability  leading  to 
a  rapid  rate  of  wave  amplification  in  the  one  dimensional  case  and 
a  rapid  rate  of  cyclogenesis  in  the  two  dimensional  system. 

fls  illustrated  by  Fig.  44,  this  hypercritical  instability  ensues 
when  the  vertical  motion  of  the  amplifying  perturbations  brings 
about  condensation  in  the  ascending  branches  of  the  waves,  causing 
the  stability  criterion  to  change  abruptly  from  — j2-—  to  -^2-S—  . 

P  P 
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FIGURE  44.  Illustration  of  Hypercritical  Instability  of  Atmospheric  Waves  Resulting 
from  Wave-Produced  Condensation  (After  Mintz.  1961) 


Because  s*/s  -  0.6,  about  two- thirds  of  all  the  potential  energy 
that  was  slowly  stored  in  the  buildup  of  the  temperature  gradient 
to  its  magnitude  —  TjS  now  is  available  for  rapid  conversion  to 
kinetic  energy.  Therefore,  in  this  model  atmosphere,  waves  which 
would  grow  only  slowly  and  to  moderate  amplitudes  (rms  V2  =  13  m  sec 
when  the  temperature  gradient  is  brought  past  its  critical  value 
from  below  and  no  condensation  is  permitted,  will  with  the  conden¬ 
sation  grow  rapidly  and  for  a  shorter  time  to  larger  amplitudes. 

With  the  larger  amplitudes,  H  »  AQ  and  the  poleward  temperature 

”  7? 

gradient  falls  rapidly  to  the  saturation- critical  value  ~  TgS*. 

Only  after  (  -|^)  falls  below  do  the  waves  die  oBt ,  so 

that  the  cycle  may  begin  again  with  —  TgS  again  as  the  determining 
value  for  the  next  releave  of  hypercritical  instability. 


-1, 
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The  following  Table  9  shews  Mintz'  theoretically  predicted 
values  of  the  critical  'and  hence  the  limiting)  poleward  temperature 
gradients,  the  dominant  wave  numbers,  and  the  root  mean  square 
meridional  wind  for  his  model  of  the  earth’s  atmosphere,  compared 
with  the  observed  values  for  the  real  earth  atmosphere.  Both 
sets  are  yearly  averages. 


TABLE  9.  THEORETICAL  AND  OBSERVED  VALUES  OF  CRITICAL  POLEWARD 
TEMPERATURE  GRADIENT  DOMINANT  WAVE  NUMBER 
AND  MERIDIONAL  WIND  (After  Mintz,  1961) 


PREDICTED 

OBSERVED 

(»*  /cn,,S  =  0.5'38C°/">« 

(a,2).*.,  S-04=23C°/rad 

(^observed  =  29C°/rad 

nd(s  —  0.5)  =(5-4>"5 
nd(s*  =  0.3)  =  7 

nd,  observed 

rms(V2)Hw=AQ~9mS"1 

“lVdra!jslu-i'1 

The  agreement  between  the  theoretically  derived  quantities  for 
the  Mintz  model  atmosphere  and  the  observed  quantities  for  the  real 
atmosphere  of  the  earth  is  good.  The  observed  poleward  temperature 
gradient  falls  almost  halfway  between  the  two- limiting  values,  the 
observed  dominant  wave  number  is  in  close  agreement  with  the  theo¬ 
retical  numbers,  and  the  observed  rms  meridianal  wind  is  approximately 
the  same  as  the  theoretical  estimate. 


With  respect  to  seasonal  changes  of  the  general  circulation, 
Mintz  notes  first  that  the  intensity  of  the  atmospheric  waves 
depends  direg^ly  upon  AQ.  Neither  the  critical  poleward  temperature 
gradient  (  and  hence  the  actual  poleward  temperature 

gradient,  nor  the  dominant  wave  number  n^  depend  directly  upon  AQ. 

Using  the  measured  poleward  transport  of  sensible  heat  at  the 
central  latitude,  which  is  about  four  times  greater  in  summer  than 
in  winter,  to  arrive  at  a  rough  estimate  of  the  seasonal  values 
of  the  differential  heating  of  the  atmosphere,  Mintz  obtains 

AQ  ~  2(10)12  XJ  sec-1 

summer 

AQwinter  ~  8(10)12  sec--"  (2.200) 

From  this,  Mintz  concludes  that  the  rms  meridianal  wind  must  be  twice 
as  large  in  winter  as  in  summer;  this  is  approximately  what  is 
observed. 

On  the  other  hand,  the  static  stabilities,  s  and  s*,  change 
relatively  little  with  season  in  the  earth's  atmosphere,  since 
the  greater  winter  upward  heat  transport  is  compensated  by  the 
greater  heat  of  condensation  produced  by  the  waves  and  released  in 
the  lower  half  of  the  atmosphere.  Observations  show  that  at  the 
central  latitude, 

(~5r)  =  22°c/rad 

\  /summer 

(~4?r)  =  36°C/rad  (2.201) 

\  ^  /winter 

This  seasonal  change  of  poleward  temperature  gradient  is  relatively 
small  compared  to  the  fourfold  increase  in  AQ.  The  dominant  wave 
number  depends  inversely  on  the  square  root  of  the  static  stability. 
Therefore,  in  spite  of  the  large  seasonal  change  in  AQ,  the  wave 
number  can  decrease  only  slightly  from  summer  to  winter  (from  7  or  8 
to  5  or  6) ,  as  is  observed. 
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Effects  of  the  oceans  and  continents 


Thei'e  are  several  important  ways  in  which  the  oceans  influence 
the  general  circulation  of  the  atmosphere  of  the  earth.  One  of  these 
is  that  evaporation  from  the  oceans  keeps  the  atmosphere  in  a  state 
of  near  saturation.  With  its  characteristic  temperatures,  the 
order  of  magnitude  of  which  is  given  by  the  effective  insolation 
S(1  -  A),  near  saturation  means  sufficient  water  vapor  in  the 
earth's  atmosphere  to  make  the  infrared  emission  to  space  almost 
constant  with  altitude.  As  a  result,  the  differential  heating 
AQ.  exceeds  the  critical  value  for  a  stable  symmetric  circulation  and 
the  wave  regime  results. 

A  second  important  effect  of  water  in  the  earth's  atmosphere 
is  that  once  the  static  stability  s  is  reduced  to  its  new  value  s* 
by  condensation  in  the  waves,  the  heat  of  condensation  thereafter 
keeps  s*  constant.  In  the  wave  regime  there  is  a  positive  correlation 
between  the  ascending  motion  and  the  temperature  to  produce  a  net 
upward  transport  of  heat  which,  by  itself,  would  make  the  space- 
averaged  vertical  gradient  of  temperature  more  stable.  The  heat  of 
condensation,  however,  is  almost  entirely  released  in  the  lower  half 
of  the  atmosphere,  in  contrast-  with  the  principal  cooling  term,  the 
divergence  of  the  net  vertical  infrared  radiative  flux,  which  is 
relatively  constant  with  elevation.  Being  released  in  the  lower 
half  of  the  atmosphere,  the  heat  of  condensation  therefore  opposes 
the  reduction  of  the  space  averaged  vertical  temperature  gradient  by 
the  upward  heat  transport  of  the  waves.  In  this  way,  s*  remains 
relatively  constant  in  value. 

In  addition  to  the  effect  of  evaporation  from  the  ocean  producing 
the  water  vapor  of  the  atmosphere ,  the  oceans  also  importantly 
influence  the  atmosphere  through  the  seasonal  storage  of  insolation. 
Since  little  insolation  is  absorbed  by  the  earth's  atmosphere 
directly,  by  far  the  larger  part  of  the  effective  insolation  S(1  -  A) 
is  absorbed  by  the  underlying  surface.  Because  land  has  small  thermal 
conductivity,  especially  when  dry,  little  heat  is  stored  in  the  land 
as  the  season  changes.  In  contrast,  the  upper  layer  of  the  oceans 
has  a  high  chermal  eddy- conductivity  and  produces  a  large  seasonal 
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storage  of  insolation.  Since  the  temperature  of  the  ocean  surface 
changes  little  with  the  season,  the  heat  transfer  from  ocean  to  over- 
lying  atmosphere  has  relatively  seasonal  variation.  If  a  planet  were 
completely  ocean  covered,  the  would  have  only  a  small  seasonal 
variation.  On  earth  there  is  sufficient  ocean  area,  even  the  northern 
hemisphere,  for  tQ  to  maintain  a  preponderance  over  throughout 

the  year,  and  thereby  maintain  the  normal  wave  regime  year  around. 

As  a  consequence  of  the  storage  process  and  the  limited 
longitudinal  extent  of  the  oceans,  there  is  the  important  effect 
of  the  seasonal  storage  of  insolation  by  the  oceans  in  producing 
large  zonal  variations  in  the  heating  of  the  atmosphere.  In 
the  winter  of  the  northern  hemisphere,  for  example,  in  middle  and 
high  latitudes  there  is  strong  relative  cooling  of  the  air  over 
the  continents  and  heating  over  the  oceans  as  the  air  of  the  tropo¬ 
sphere  flows  from  west  to  east.  Over  the  continents  the  air  cools 
because  S(1  -  A)  is  smaller  than  W,  while  the  temperature  of  the 
ground  adjusts  itself  to  the  temperature  of  the  air  because  the 
thermal  conductivity  of  the  ground  is  small.  When  the  air  that  has 
been  cooled  over  a  continent  moves  over  the  ocean  again,  it  is  heated 
because  it  is  colder  than  the  warm  ocean  surface,  while  the  temperature 
of  the  ocean  surface  remains  almost  unchanged  because  the  thermal 
eddy- conductivity  of  the  oceans  is  large.  There  are  two  important 
consequences  of  this  zonally  varying  heating  of  the  atmosphere. 

One  is  that  low-level  convergence  is  produced  in  the  longitudes 
where  there  is  heating  of  the  air  and  low-level  divergence  when  there 
is  cooling.  As  a  result,  in  the  winter  season  there  will  be  general 
low  level  horizontal  convergence  over  the  oceans  and  divergence 
over  the  continents,  and  therefore,  lower  mean  surface  pressure  over 
the  oceans  than  over  the  continents.  In  the  summer  season,  the 
reverse  pattern  of  heating  and  cooling  and  reverse  distribution  of 
high  and  low  surface  pressures  over  oceans  and  continents  will  occur. 

At  the  same  time,  the  heating  established  in  the  upper  levels  an 
ultra- long- wave  pattern  of  opposite  phase  to  the  surface  pressure 
distribution. 
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The  second  important  consequence  of  the  zonally  varying  heating 
due  to  the  ocean- land  contrast  is  that  it  affects  the  horizontal 
temperature  gradients,  and  hence  the  baroclinic  stability  at  the 
coast  lines.  Especially  near  the  eastern  coastlines  of  the  two 
northern  hemisphere  continents  in  the  winter  season,  the  contrast 
Detween  cold  continental  air  and  warm  maritime  air  produces  a 
zone  of  large  horizontal  temperature  gradient  and  hence  a  zone  of 
maximum  baroclinic  instability.  Thus,  the  eastern  coastlines 
become  preferred  regions  of  baroclinic  wave  development.  As  a 
baroclinic  disturbance  moves  eastward  from  its  source  region  of 
maximum  horizontal  temperature  gradient  and  moves  through  the 
region  of  maximum  heating,  and  consequent  surface  convergence, 
of  the  ultra- long  atmospheric  wave,  the  low-level  cyclonic  part  of 
the  moving  disturbance  is  amplified.  In  this  way  the  baroclinic 
wave- cyclones ,  instead  of  being  uniformly  distributed  around  the 
globe,  are  concentrated  in  the  oceanic  longitudes,  especially  in 
the  winter  season. 

Still  another  important  influence  of  the  ocean  lies  in  its  own 

ability  to  transport  heat.  The  mean  latitudinal  heat  transport  by 

19  “"1 

the  ocean  is  indirectly  estimated  to  be  about  1(10)  cal  day  at 

the  central  latitude  of  the  northern  hemisphere,  and  hence  has  no 

controlling  on  the  form  of  the  atmospheric  general  circulation, 

19  -1 

since  a  reduction  by  that  amount  from  the  10.6(10)  cal  day 
required  poleward  heat  transport  would  not  change  the  circulation 
from  wave  to  symmetric  regime.  However,  important  aspects  of  the 
atmospheric  circulation,  within  the  wave  regime,  can  be  influenced 
by  the  oceanic  transport  of  heat.  As  a  hypothetical  example,  an 
anomalous  excess  of  insolation,  resulting  from  a  period  of  non¬ 
cloudiness,  could  be  stored  in  the  ocean  in  a  region  of  down- 
welling  current.  This  anomalous  heat,  when  brought  to  the  surface 
and  transferred  to  the  air,  could  affect  the  phase  of  the  ultra- 
long  atmospheric  waves.  Because  of  the  interaction  between  the 
ultra- long  waves  and  the  baroclinic  disturbances,  this  could 
result  in  a  shift  of  the  tracks  of  the  rain  producing  baroclinic 
disturbances,  east  or  west,  resulting  in  long- period  weather 
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anomalies  over  extensive  regions  of  the  globe.  The  atmosphere  s 
small  heat  capacity  and  large  velocity  gives  the  atmosphere 
a  short  memory  for  anomalies  of  insolation*  but  the  ocean  with 
its  large  heat  capacity  and  small  velocities,  has  a  long  memory 
for  anomalies  of  insolation. 


III. 


PHYSICS  OF  THE  STRATOSPHERE  AND  MESOSPHERE 


Extending  in  altitude  between  the  two  temperature  minima  of  the 
temperature-altitude  profile  illustrated  in  Fig.  4,  the  stratosphere 
and  mesosphere  form  a  region  from  perhaps  15  to  80  km  altitude, 
characterized  by  a  temperature  maximum  in  the  40  to  60  km  range, 
which  acts  somewhat  as  a  barrier  to  the  flow  of  energy,  in  forms 
other  than  radiative  and  conductive,  between  the  two  regions  of 
the  atmosphere  v/hich  bound  it.  Energy  within  the  region  - ~  derived 
mainly  by  absorption  of  solar  ultraviolet  radiation  by  ozone.  Wave 
motion  in  the  stratosphere  and  mesosphere  is  strongly  suppressed. 

TEMPERATURE  STRATIFICATION 

The  temperature  maximum  is  due  to  the  absorption  of  solar  UV 
in  the  2000-3000  A  region,  largely  by  the  Hartley  system  of  ozor.0 
for  which  the  absorption  cross  section  peaks  near  2600  a  (See  Fig.  19). 
The  transition  corresponding  to  the  photodissociation  of  ozone  is 


0^  +  h  v  (Hartley)  -  0^  -*  0*; 


-2  -1 

1  x  10  sec  above  40  km. 


(3.1) 


At  the  altitude  near  50  km  the  atmosphere  is  optically  thick 
for  this  wavelength  region,  that  is 


§(h)  =  1  for  h  a  50  km 

wherein  5(h)  =  ^(Hartley)  /  n(0,;hr)  dh» 

Jh  3 

MKfillNfi  m  BUNK 


(3.2) 
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mm 


n(0^;  h)  is  the  particle  number  density  of  ozone  at  altitude  h,  as 
given  by  Nagata  et  al.  (1957)  and  plotted  in  Fig.  45. 


FIGURE  45.  Ozone  Concentration  as  a  Function  of  Height 
(After  T.  Nagata  et  al.,  1967) 


Although  nfO^ ;  h)  peaks  at  h  =  20  km,  since  §(h)  =  1  already 
at  50  km,  the  peak  heating  occurs  at  the  relatively  higher  altitude. 
Evidently,  above  the  tropopause,  the  ozone  concentration  may  be 
explained  by  "photochemical  equilibrium",  that  is,  the  ozone  concen¬ 
tration  is  determined  by  the  balance  between  the  Hartley  absorption 
reaction  (3.1)  and  the  processes 

0  +  0  +  M  *  02  +  M;  k4; 
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jAbove  the  upper  mesosphere,  k^  =-  3  x 


10-33cm6sec' 


*] 


(3,3) 


0  +  0^  +  M  ^  0^  +  M;  ; 


k  **« 

k5  _ 


5.6  x  10-11  exp 


c^sec-1] 


(3.4) 


0  +  03  -  2  0*;  kg; 


jijelow  the  lower  mesosphere;  kg  8.2(10)  33  exp^^^cm^sec 


02  +  M  **  0*  +  M 


(3.5) 

(3.6) 


Under  steady-state  conditions  for  which 


(d/dt)  n(0)  =  (d/dt)  n(02)  =  (d/dt)  n(03)  =  0  (3.7) 

At  altitudes  of  50  km  and  above  the  concentration  of  ozone  is 
drastically  reduced  because  of  the  higher  intensity  of  solar  UV 
which  leads  not  only  to  the  ozone  photodissociation  (3.1),  but  also 
to  the  02  dissociation  reactions: 

02  +  h  v  (Schumann- Runge)  -  0(3P)  +  0(1D);  cig; 

jo9  <  10~17  cm2J  (3.8) 

02  +  h  v  (Herzberg)  -  2  0(3P);  a1Q; jo1Q  ^  10~23  cm2J  (3.9) 

The  details  of  the  photochemistry  of  atomic  oxygen  and  ozone, 
in  the  stratosphere  are  complicated  by  reactions  involving  water 
LHunt,  1966  and  Leavy,  1969], 

At  altitudes  below  20  km,  i.e.,  in  the  troposphere,  the  ozone 
concentration  is  higher  than  can  be  explained  in  the  terms  which 
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describe  it  for  the  stratosphere  and  mesosphere.  Evidently  this 
high  concentration  arises  from  low- Level  photochemistry  of  a  differ¬ 
ent  nature,  such  as  that  generating  smog,  etc. 

To  summarize,  at  altitudes  above  50  km  5(h)  <  1  because  there 
is  too  little  ozone  ro  make  the  atmosphere  optically  thick;  at 
altitudes  below  50  km,  since  most  of  the  solar  radiation  incident 
on  the  atmosphere  has  already  been  absorbed  at  higher  altitudes 
(5  >  1),  there  is  inadequate  solar  radiation  in  the  Hartley  band  to 
heat  up  the  atmosphere. 

Having  explained  the  temperature  maximum  at  the  stratopause , 
the  lower  bound  of  the  stratosphere  and  the  upper  bound  of  the 
mesosphere  will  next  be  described. 

The  lower  bound  of  the  stratosphere  is  given  by  the  tropopause, 
which  corresponds  to  the  temperature  minimum  between  the  tropopause, 
or  that  portion  of  the  atmosphere  which  is  in  direct  "contact"  with 
the  surface  of  the  earth,  and  the  thermopause. 

If  the  temperature  gradient  (with  increasing  height)  within 
the  troposphere  is  negative,  so  that  the  adiabatic  lapse  rate 

-dT/dh  <  g/Cp  (3.10) 

(where  g  =  acceleration  due  to  gravity  and  c^  =  specific  heat  at 
constant  pressure),  there  is  possibility  of  a  fluid  mechanical 
static  instability  by  which  a  parcel  of  air  of  slightly  different 
temperature  can  move  a  significant  vertical  distance,  leading  to 
overturning  of  the  atmosphere.  Such  an  overturning  in  the  tropo¬ 
sphere  is  particularly  visible  in  the  case  of  cumulus  clouds. 

The  vertical  motion  arising  from  static  instability  is  of 
vital  importance  in  keeping  water  droplets  or  ice  crystals  in 
suspension  in  the  troposphere.  Since  the  terminal  velocity  of  a 

spherical  water  droplet  of '  Ip  radius  predicted  by  the  Stokes’ 

-2 

formula  is  about  10  cm/sec,  modest  mean  vertical  motions  may  keep 
such  droplets  in  suspension. 
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The  clouds  made  up  of  ice  and  water  particles  fulfill  an 
essential  role  in  maintaining  the  thermal  balance  of  the  troposphere 
by  absorbing,  re-radiating  and  scattering  solar  and  terrestrial 
radiation,  as  well  as  by  the  large  thermal  inertial  provided  by 
their  heat  capacity. 

There  are  rare  occurrence^  of  other  types  of  clouds  at  higher 
altitudes- "nacreous  clouds",  which,  occur  rarely  between  20  and  30  km 
under  certain  geographic  and  climatic  conditions  associated  with 
strong  winds,  and  "noctilucent  clouds"  which  occur  occasionally 
between  75  and  90  km,  also  associated  with  winds. 

Under  meteorological  conditions  the  criterion  for  the  mainte¬ 
nance  of  turbulence  is  commonly  expressed  in  terms  of  the  dimension¬ 
less  Richardson’s  Number  Ri,  which  is  the  ratio  of  the  restoring 
buoyancy  force  per  unit  length  of  vertical  displacement  to  the 
square  of  the  wind  shear 

(g/TKST/az  +  g/c  ) 

Ri  =  - 5 - E_  (3.11) 

(du/dz  ) 

where  u  is  the  horizontal  wind  velocity.  Turbulence  persists  if 

Ri  <  Ric  (3.12) 

where  various  values  of  Ric  have  been  quoted,  ranging  from  2  down 
to  0.4  (Schlichting,  1960). 

Typical  values  of  Richardson  number,  as  computed  for  '  arious  al¬ 
titudes,  are  plotted  in  Fig.  46.  Also  indicated  are  the  v  ues  of 
temperature  and  of  wind  velocity  used  in  the  computations,  in  Fig. 
46(a),  Kochanski’s  wind  data  (Kochanski,  1964,  1966)  are  taken  with 
temperature  profiles  of  the  U.S.  Standard  Atmosphere  to  be  a  basis  for 
estimates  of  Reynolds  and  Richardson  numbers.  In  Fig.  46(b)  are 
Zimmermann’s  computations  of  Richardson  number  based  on  Sechrist 
et  al.  (1969)  data  of  simultaneous  measurements  of  wind  and  tempera¬ 
ture.  Zimmermann’s  data  show  that  turbulence  may  occur  low  in  the 
mesosphere  under  some  conditions . 


(a)  COMPUTATION  USING  KOCHANSKi  ( 1964,  1966)  WIND  DATA  AND  U.S. 

STANDARD  ATMOSPHERE,  1962,  TEMPERATURES  FOR  JANUARY,  LATITUDE  60°N 
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Conclusions  to  be  drawn  from  the  data  on  Richardson  Number  ar¬ 
rayed  in  Fig.  46  is  that  turbulence  in  the  troposphere  is  unambiguously 
present,  that  the  existence  of  turbulence  in  the  stratosphere  and 
mesosphere  is  unlikely,  and  chat  turbulence  in  the  lower  thermosphere 
is  near  the  critical  limit  for  suppression  by  the  temperature  gradient. 

That  the  tropopause  is  sometimes  abrupt  is  evidenced  by  the 
distinctive  flat  anvil-top  of  thunderclouds  delineating  an  altitude 
above  which  condensed  water  cannot  exist  stably,  and  where  the 
instabilities  discussed  above  no  longer  occur.  Above  the  tropopause 
there  may  be  strong,  steady,  prevailing  winds  driven  by  the  rotation 
of  the  earth  as  well  as  by  the  day/night  and  equator/pole  heating 
and  cooling. 

The  upper  limit  of  the  mesosphere  is  set  by  the  mesopause,  the 
minimum  in  temperature  below  the  thermosphere  which  is  heated  largely 
by  the  Schumann- Runge  continuum  leading  to  the  photodissociation 
of  0o ,  and  above  the  stratopause  arising  from  the  heating  ascribed 
to  the  Hartley  continuum  of  ozone.  The  inversion  due  to  the 
Schumann-Runge  heating  represents  a  smaller  heat  source  than  that 
of  the  stratopause,  so  that  the  conductive  heat  flow  in  the  meso¬ 
sphere  is  from  below  upwards  along  a  negative  temperature  gradient. 

ENERGY  BALANCE  AND  FLOW 

As  in  the  D-region  of  the  lower  thermosphere,  the  chemistry  of 
the  stratosphere  and  mesosphere  is  complex,  involving  not  only  the 
minor  species  ozone  and  its  decomposition  products  as  indicated  in 
Eqs.  3.1,  3.3,  3.4,  3.5,  3.6,  3.8,  and  3.9,  but  also  hydrogen  com¬ 
pounds  resulting  from  the  dissociation  of  water  vapor. 

Horizontal  transport  of  energy  may  be  effected  by  means  of  strong 
prevailing  horizontal  winds,  of  speeds  50-100  m/sec  (CIRA,  1965,  p. 

41  ff). 

Figure  47  shows  schematically  the  vertical  energy  transfer  from 
the  sun  and  within  the  various  upper  atmospheric  regions.  Other 
than  for  radiation  directly  from  the  sun,  from  region  to  region,  the 
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transport  is  generally  vertically  upward,  not  downward.  Since 
density  decreases  with  height  energy  transport  from  a  higher  region 
is  little  able  to  affect  a  lower  region.  The  higher  region,  due 
to  its  lower  density,  always  has  smaller  heat  capacity  than 
the  lower  region,  the  higher  region  cools  by  radia'  ve  emission 
rather  than  conductively  heats  the  lower  region. 

The  following  mechanisms  for  vertical  energy  transport  have 
been  considered  - 

(a)  radiative  transport;  given  the  prevailing 
temperatures  200  -  300°  K,  the  vibration-rotation 
levels  in  the  infrared  (2  -  30  u)  are  most  important 
to  this  process 

(b)  molecular  diffusivity 

(c)  turbulent  diffusivity 

(d)  transport  by  organized  macroscopic  flows  such 
as  vertical  winds  or  gravity  waves 

(e)  chemical  transport  (  e.g.,  atomic  and  molecular 
oxygen  are  transported  from  one  altitude  to 
another  where  they  recombine,  releasing  the 
bond-energy). 

Of  these  mechanisms,  chemical  transport  (e)  can  clearly  be 
considered  negligible  because  of  the  low  degree  of  dissociation. 
Molecular  diffusivity  (b)  can  be  calculated  by  the  standard  tech¬ 
niques  (Chapman  &  Cowling,  1339;  Hirschfelder,  Curtiss  &  Bird,  1954). 
Radiative  transport  (a)  is  possibly  the  most  important  single 
mechanism.  By  the  theory  of  chemical  reactions ,  quenching  and  the 
difficult  conversion  of  vibrational  into  translational  energy,  perhaps 
half  of  the  solar  excitation  absorbed  remains  in  the  form  of  vibra¬ 
tional  energy.  The  transfer  of  vibration  energy -to  form  infrared- 
active  species,  which  then  re-radiate,  tends  to  be  rapid  compared 
to  vibration-translation  deactivation.  By  these  simplifying  con¬ 
siderations,  the  radiative  transport  nay  be  calculated  by  the  gen¬ 
erally  recognized  techniques  of  Kondrat’ev.  (1965)  and  Goody  (1964) 
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after  the  appropriate  partition  of  energy  from  solar  input  into 
infrared  active  vibrational  and  rotational  energy  has  been  made. 

The  large-scale  mechanical  transport  of  energy,  either  by 
ordered  motions,  such  as  vertical  winds  or  various  other  types  of 
waves,  or  by  random  motion  (molecular  diffusion  or  turbulence). 
Because  of  static  stability  considerations  (high  Richardson’s  Number) 
random,  or  non-phase  coherent  motions,  appear  not  to  be  important 
for  vertical  transport  of  energy  in  the  stratosphere,  in  sharp 
contrast  to  the  situation  in  the  troposphere.  Some  phase  coherent 
(ordered)  motions  of  importance  in  the  stratosphere  and  mesosphere 
are  the  driven  motions  of  the  tides  (with  periods  of  fractions  of 
a  solar  day)  and  of  the  plentary  waves,  with  longer  periods.  These 
are  discussed  in  greater  length  in  the  following  section. 


MECHANICAL  TRANSPORT  BY  WAVE  MOTION 

One  may  consider  the  large-scale  mechanical  transport  of  energy 
to  be  accomplished  by  wave  action,  which  is  phase  coherent,  and  by 
turbulence,  which  is  not  phase  coherent.  Eckart  (1960)  nas  given  an 
excellent  description  of  the  several  modes  of  wave  action  indicating 
that  frequencies  greater  than  the  Brunt -Vaisala  frequency 

| the  waves  propagate  like  sound  waves  (hence  are 


H(M) 


termed  "acoustic  waves"),  frequencies  smaller  than  the  Brunt-Vaisala 
frequency  but  larger  than  the  Coriolis  frequency  (0=2  day-1)  are 
characterized  by  the  effect  of  gravity  on  buoyancy  dominanting  over 
the  spring-like  compressibility  effect  (hence  are  called  "gravity 
waves").  Some  of  the  gravity  waves  are  propagated  in  modes  with 
zero  amplitude  at  the  boundaries  (hence  they  are  termed  "internal 
gravity  waves"). 


Waves  of  another  type  are  those  called  by  Laplace  "oscillations 
of  the  second  kind”  and  include  the  tidal  oscillations  and  the 
planetary  waves.  The  tidal  oscillations  have  periods  corresponding 
to  1/x  of  the  solar  or  lunar  day  (X  =  1,  2,  3)  for  frequencies  in  the 
order  of  the  Coriolis  frequency;  planetary  waves,  which  have  also 
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been  described  in  one  form  by  Rossby,  have  smaller  frequencies  than 
the  Coriolis  frequency. 

The  major  period  of  the  atmospheric  tidal  oscillations  is  half 
a  solar  day,  rather  than  half  a  lunar  day,  as  would  be  expected  from 
the  application  of  Newton's  theory  of  gravitation.  Siebert  (1961) 
has  shown  that  solar  heating  plausibly  gives  rise  to  the  semi-diurnal 
tidal  oscillation,  with  the  diurnal  solar  heating  tide  suppressed  by 
anti-resonance  effects.  Adding  to  Siebert 's  calculations  for  the 
effect  of  water  vapor  absorption,  S.  T.  Butler  and  K.  A.  Small  (1963) 
have  shown  that  the  water  vapor  and  ozone  absorption  together 
can  adequately  account  for  the  observed  ground  level  atmospheric 
semi-diurnal  and  terdiurnal  tides. 

Theory  and  observation  described  by  Charney  and  Drazin  (1961) 
both  indicate  that  only  components  of  planetary  waves  (time  scale 
long  compared  to  a  day),  which  are  of  global  extent  can  propagate 
upward  out  of  the  troposphere  into  the  stratosphere  and  mesosphere. 

R.  E.  Dickinson  (1968)  has  described  the  following  characteristics 
of  the  terrestrial  planetary  waves: 

1.  Planetary  scale  disturbances  interact  with  the  vorticity 
of  the  earth's  rotation  and  propagate  through  the 
atmosphere  as  planetary  waves. 

2.  Planetary  waves  are  everywhere  evanescent  in  the  presence 
of  an  easterly  zonal  wind. 

3.  Strong  westerly  winds  trap  the  vertical  propagation  of 
planetary  scale  waves. 

4.  Trapping  of  vertical  propagation  of  planetary  waves  is 
absent  locally  near  the  equator. 

5.  Planetary  waves  are  guided  by  horizontal  shears.  Examples 
are  illustrated  in  Fig.  48  which  is  a  schematic  of  the 
winter-summer  zonal  winds  and  the  trajectories  of  planetary 
wave  rays  in  the  presence  of  the  winds.  The  maximum  westerly 
jet  in  middle  latitudes  at  the  stratopause  acts  as  a  reflecting 
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DAMPING  TIME  a  (DAYS)  ALTITUDE  (KM) 


r 


ZONAL  WIND  ( M/SEC )  ZONAL  WIND  ( M/SEC ) 


WAVE  NUMBER  ONE  WAVE  NUMBER  TWO 


DAMPING  OF  PLANETARY  WAVE  BY  ZONAL  WINDS 


FIGURE  48.  Schematic  of  Winter-Summer  Zonal  Winds  and  Trajectories  of 
Planetaiy  Wave  Guides  (After  R.  E.  Dickinson,  1968) 


barrier.  The  latitudinal  variation  of  winds  suggests  the 
formation  of  polar  and  equatorial  wave  guides  where  the  winds 
are  weaker.  Rays  in  the  polar  wave  guide  undergo  multiple 
reflections  between  the  pole  and  the  westerly  jet,  setting 
up  a  standing  wave  normal  mode  in  latitude.  Also,  in  the 
equatorial  wave  guide,  the  rays  are  deflected  onto  the  zero 
wind  line  where  they  are  absorbed,  giving  rise  to  rapid  at¬ 
tenuation  with  height. 

6.  Planetary  waves  in  weak  zonal  winds  are  strongly  damped  by 
radiative  processes.  In  Fig.  48  (after  R.  E.  Dickinson,  1968) 
is  shown  the  damping  time  for  wave  numbers  one  and  two 
versus  zonal  wind  for  a  number  of  values  of  the  parameter 
representing  fractional  damping  in  one  density  scale  height. 
The  shaded  region  of  parameter  space  indicates  where  the 
damping  is  severe  enough  to  cause  wave  amplitudes  to  decrease 
with  height. 

7.  Zonal  motions  may  be  coupled  to  both  solar  radiation  heating 
and  planetary  waves. 

Of  the  several  types  of  wave  motions  listed  above,  the  planetary 
waves  seem  to  dominate  the  circulation  of  the  troposphere;  the  planetary 
waves  and  tides  dominate  the  stratosphere  and  mesosphere.  In  the 
lower  thermosphere,  in  contrast,  the  gravity  waves  are  important. 

Gravity  waves  are  observed  within  the  troposphere,  stratosphere, 
mesosphere,  and  the  lower  thermosphere.  They  give  rise  to  the  wind 
fluctuations  which  are  seen  at  low  levels  as  gusts.  Observed  at 
higher  altitudes  as  deflections  of  smoke  trails,  they  have  impressive 
magnitude  in  the  lower  thermosphere,  where  Kochanski’s  (1966)  analysis 
shows  them  to  be  comparable  in  amplitude  to  the  tides. 

In  the  first  order,  in  the  lower  thermosphere  and  below,  by 
neglecting  the  effects  of  viscosity  and  assuming  the  conservation  of 
energy  as  the  wave  motion  is  propagated  upward,  one  may  note  that 
since  the  density  decays  upward  as  exp  (-z/H),  the  wave  velocity- 
must  increase  upward  as  exp  (z/2h‘).  The  tides  and  the  gravity  waves 
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accordingly  may  be  expected  to  have  velocity  values  three  orders  of 
magnitude  greater  at  100  km  than  at  the  earth’s  surface.  At  50 -km  al¬ 
titude,  on  the  other  hand,  gravity  wave  and  tidal  velocities  may  be 
but  little  more  than  an  order  of  magnitude  greater  than  at  the  surface. 

Mahoney  (1968)  has  descriced  results  of  tracking  falling  ROBIN 
spheres  by  two  different  radars  simultaneously,  to  show  as  in  Fig.  49 
that  the  wave  motion  common  to  the  two  tracks  of  the  same  sphere 
(tracks  1  and  3  in  the  figure)  is  of  finite,  but  small  amplitude  in 
the  stratosphere ,  and  correlates  well  with  the  motion  of  a  second 
ROBIN  sphere  falling  at  the  same  time  about  48  km  distant . 

Kochanski’s  (1964)  data  shown  in  Fig.  50  show  values  of  winds  at 
Fort  Churchill  (60°  N)  and  Wallops  Island  (30°  N)  and  elsewhere.  The 
large  amplitudes  of  wind  in  the  stratosphere  and  mesosphere  are  thought 
to  be  mainly  due  to  tidal  and  planetary  wave  effects. 

STRUCTURE  OF  ATMOSPHERIC  TURBULENCE 

The  study  of  irregular  motions  in  the  atmosphere,  by  means  of 
rocket  releases  of  chemicals  and  by  means  of  meteor  trails,  has  estab¬ 
lished  a  region  between  85  and  110  km  called  turbosphere.  In  this 
region  irregular  motions  of  both  turbulent  and  non -turbulent  nature 
have  been  observed,  as  described  in  Fig.  51.  The  curve  in  Fig.  51 
represents  the  diffusion  coefficient  calculated  using  the  altitude 
variation  of  the  atmospheric  density.,  temperature,  and  molecular 
weight  given  by  the  1965  COSPAR  international  reference  atmosphere 
for  a  10.7  solar  flux  index  85,  0400  h,  while  the  two  lines  with  bars 
indicate  the  theoretical  extremes  of  variation  with  season  and  solar 
cycle  (Golomb  and  MacLeod,  1966). 

Their  investigations  by  experimental  means  have  led  to  different 
interpretations  (Zimmerman,  S.  P.,  1966;  Justus,  C.  G.,  1966)  and 
cast  doubt  to  the  existence  of  natural  turbulence  (Bedinger  and 
Layzer,  1969). 


ALTITUDE  (KM) 


WIND  TOWARD  NORTH  (METERS/SEC) 


i  \i  i  -k,  Profiles  in  Stratosphere  and  Mesosphere  by 
FIGURE  49.  Honzontal  Veloc  ty  Prof  ?  Mahoney,  1968) 


RESULTANT 


FIGURE  50.  Wind  Near  38°  N  (After  Kochanslcl,  1964) 


Measurements  at  104  km  indicate  an  energy  dissipation  rate  of 

3  6 

e  =  5  x  10  cgs  and  a  molecular  viscosity  of  v  =  10  ,  yielding  an 

internal  scale  of  turbulence  (see  Fig.  51). 


k-1  =  (v3/e>1/4 

v 


3  x  10  cm 


The  maximum  horizontal,  scale  of  turbulent  inhomogeneities  may  be  esti¬ 
mated  of  the  order  3000  km,  and  the  vertical  scale  of  the  order  of  the 
scale  height  H.  This  indicates  that  the  spectrum  of  atmospheric 
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turbulence  is  very  broad,  compared  to  spectra  in  laboratory  turbulence, 
suggesting  that  one  ought  to  find  some  portions  of  the  spectrum  obeying 
the  Kolmogoroff  law  (Kolmogoroff,  1941a,  1941b): 


£ 

i! 


c2/3  k-S/3 


(3.13) 


Lower  altitude  flights  over  Australia  confirm  such  a  5/3  -  law,  see 
Fig.  52  (Reiter,  E.  R.  and  Burns,  A.,  1966).  However,  at  high  alti¬ 
tudes  (85  to  100  km),  the  presence  of  wind  shear  and  buoyancy,  or 
gravity  waves,  should  complicate  the  structure  of  turbulence. '  It  has 
been  contemplated  that  the  turbulence  may  be  generated  by  internal 
gravity  waves  (Hodges,  1967). 


FIGURE 51 .  Diffusion  Coefficient  (After  Golomb  and  MacLeod,  1966) 
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WAVELENGTH,  ft 

FIGURE  52.  Spectra  of  longitudinal  (U)  Component  of  Turbulence  Measured  by 
Project  TOPCAT  Over  Australia  (After  Reiter  and  Burns,  1966) 

The  characteristic  frequency  N,  given  by 

N2  =  2.(2.  Jti  + 

T\R  y  dz  / 

determines  the  mode  of  the  gravity  wave,  and  is  called  the  Brunt- 
VctisSla  frequency.  It  serves  also  to  derive  the  buoyancy  subrange  of 
the  turbulent  spectrum,  which  is  found  on  a  dimensional  argument  to  be 


Here  T  is  the  temperature,  R  the  gas  constant,  and  y  the  ratio  of 
specific  heat.  Measurements  in  the  troposphere  have  found  a  portion 
of  spectrum  obeying  the  k-3  law  (Tchen,  1969a),  s»e  Fig.  53.  Spectra 
in  the  altitude  of  85  to  160  km  have  also  been  compiled  by  Rosenberg 
(1967),  after  normalization  to  their  scale  heights,  see  Fig.  54.  The 
data  includes  coverage  of  both  upper  and  lower  layers;  that  for  upper 
layers  may  include  the  effects  of  gravity  waves .  According  to  the 
reduction  by  Tchen  (1969a),  these  spectra  follow  also  the  k-3  law. 

See  also  Zimmerman  (1969),  Figs.  55  and  56. 

It  is  evident  that  the  wind  shear  should  be  strong  enough  at  those 
altitudes  to  destroy  the  isotropy  and  to  supply  energy  from  the 
permanent  wind  shear  to  the  turbulent  motion.  One  may,  therefore,  ask 
whether  the  -turbulent  spectrum  under  such  a  wind  shear  would  alter  the 
Kolmogoroff  law  (Eq.  3.13)  or  the  buoyancy  law  (Eq.  3.14).  In  this 
respect,  Tchen  (1953,  1954)  has  found  a  spectrum 

F  ~  (e/u>s)  k"1  (3.15) 

in  the  shear  production  subrange,  with  support  from  experiments  ir 
wind  tunnel  turbulence,  see  Fig.  57  and  Table  10.  This  result  has 
been  used  by  Zimmerman  (1966)  to  derive  the  characteristic  parameters 
of  turbulence.  Recently,  by  means  of  a  cascade  model  of  shear  turbu¬ 
lence,  Tchen  (1969b)  stipulates  that  the  spectrum  (Eg.  3.15)  is  valid 
in  a  wave  number  region  smaller  than  a  critical  value  k£ 

k  <  k 

s 

where 

3  1/2 

kg  =  (ii^/e)  ,  u)s  =  |dU/dz  j  =  wind  shear  frequency 

and  that  k  should  be  sufficiently  small  that  the  vorticity  function 
k 

2 J  dh’k,2F 
0 
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fluctuation  remaining 


X(SH)  3  2  I  0  5  0.3 


WAVELENGTH  OF  HARMONIC 


FIGURE  54.  Spectrum  of  Residual  Fluctuation  Amplitude  Versus  Wavelength 
of  Harmonic  Removed  (After  Rosenberg,  1967) 


149 


FIGURE  57.  Longitudinal  Energy  Spectrum  in  a  Boundary  Layer  and  in  a  Pipe 
(Tchen,  1954) 


TABLE  10.  DATA  FOR  THE  ENERGY  SPECTRUM  IN  A  BOUNDARY 
LAYER  AND  IN  A  PIPE  (Tchen,  1954) 


Experi¬ 

mental 

points 

Type  of 
flow 

6;  cm 

Urn; 

cm/ sec 

Distance 
from 
wall  j 
6/x2 

Local 

mean 

velocity 

gradient 

u’ 

cm/ sec 

c 

+ 

Boundary 

layer 

7.6 

1524 

0.05 

Large 

119 

■ 

▲ 

Pipe 

12.3 

3048 

0.008 

Large 

256 

X 

Boundary 

layer 

7.6 

1524 

0.8 

Small 

32 

• 

Pipe 

12.3 

3048 

0.69 

Small 

113 
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and  consequently  the  inertial  transfer  should  not  intervene.  By  in¬ 
creasing  k,  srill  within  the  above  region  k  <  ks,  the  inertial  transfer 
begins  to  intervene.  From  that  point  on,  the  spectral  law  (Eq.  3.15) 
is  to  be  followed  by  a  new  spectral  law 

F  ~  ujg  k~3  (3.16) 

in  the  shear  turbulence. 

_  7 

We  note  now  the  difficulty  for  the  same  k  "  law  in  Eq.  3.14  and 
Eq.  3.15  to  represent  two  different  phenomena:  one  due  to  a  buoyancy 
force  and  the  other  due  to  a  wind  shear.  A  discrimination  between  the 
two  would  be  possible  by  a  careful  normalization  of  the  experimental 
data  to  N  or  u.'s,  respectively. 

We  conclude  that  it  is  extremely  hard  to  discriminate  and 
characterize  a  turbulent  motion.  We  believe  that  it  should  possess 
the  following  features:  an  irregular  nonlinear  motion  having  a  soec- 
trum  with  a  production,  a  mode  transfer,  and  a  dissipation.  The  mode 
transfer  is  a  nonlinear  mechanism  responsible  for  the  transfer  of 
energy  across  the  spectrum  and  must  be  present  in  general. 
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IV. 


PHYSICS  OF  THE  THERMOSPHERE.  AND  IONOSPHERE 


The  region  above  the  mesopause,  called  the  thermosphere,  par¬ 
ticularly  affects  radio  communications  and  decay  of  low  orbits  of 
earth  satellites . 


GENERAL  CHARACTERIZATION 

As  shown  in  Fig.  5,  a  profile  of  the  temperature  in  the  upper 
atmosphere  as  a  function  of  altitude  has  a  minima  at  15  and  85  km. 

The  region  above  the  85  km  temperature  minimum  is  called  the  thermo¬ 
sphere.  Its  upper  limit,  commonly  placed  at  about  500  km  altitude, 
depends  significantly  on  latitude. 

The  thermosphere  is  characterized  as  follows: 

(a)  Temperature  T  increases  with  height  h,  going  from  perhaps 
180°  K  at  h  =  80  km  to  perhaps  1500°  X  (±50%)  above  400  km. 

(b)  The  number  density  decreases  from  10~^  times  its  sea-level 

-1  2 

value  at  80  km  down  to  10  ~  times  its  sea-level  value  at 
500  km;  or  from  4  x  10^4  cm-3  at  80  km  to  5  x  10^  cm-3  at 
50  km. 

(c)  There  is  a  uniformly  mixed  region  from  sea  level  up  to  per¬ 
haps  100  km.  Above  100  km,  the  composition  changes  markedly 
with  altitude. 


(d)  For  more  than  200  km  above  about  90  km  there  is  a  relatively 
uniform  density  of  free  electrons,  ng  ~  10 5  -  10  6  cm-3. 
Because  of  the  drastic  fall-off  of  neutral  density,  the 
plasma  effects'  are  important  at  the  higher  altitudes. 


These  characteristics  (a)  -  (d)  may  be  grossly  explained  as 

£ollous-  MCEEHfi  FASH  BLANK 


For  particles  of  mass  M  and  temperature  T  in  a  gravitational 
field,  the  density  falls  off  with  height  according  to: 

-(h  -  h  )/H 

P  (h)  =  p  (h  )  e  ;  H  =  kT/Mg  ~  10  km  "scale  height."  (4.1) 

As  a  result  of  photochemistry  and  flow,  this  simple  relation 
(4.1)  does  not,  in  fact,  give  the  whole  story.  It  is,  however,  a 
first  approximation. 

The  precise/specific  reasons  why  the  atmosphere  is  uniformly 
mixed  up  to  perhaps  100  km  and  not  above  are  subjects  of  present 
scientific  controversy.  The  rotation  of  the  earth  combined  with 
differential  day/night  and  equator/pole  heating  gives  rise  to  large 
scale  atmospheric  winds.  These  move  in  a  horizontal  direction  be¬ 
cause  of  gravity  but  there  are  a  number  of  large  convection  cells 
which  overall  serve  to  produce  a  large  scale  low  frequency  stirring 
and  which  also  give  rise  to  various  high  frequency  disturbances  In¬ 
cluding  "internal  gravity  waves"  and  "turbulence"  (Hines  et  al,  1965). 

The  fluid  mechanical  morions  which  lead  to  these  irregular  mo¬ 
tions  are  large  enough  to  keep  the  atmosphere  uniformly  mixed  up  to 
perhaps  100  km,  where  the  gas  kinetic  mean  free  path  is  16  cm. 

The  principal  reason  for  the  temperature  increase  above  80  km 
and  for  the  composition  change  is  that  solar  ultraviolet  and  X-radia¬ 
tion  is  absorbed  in  these  upper  fringes  of  the  atmosphere  and  leading 
to  dissociation  and  ionization,  and  consequently  heating.  At  the  80- 
150  km  level  the  biggest  absorption  mechanism  by  far  is  the  Schumann- 
Runge  continuum  of 

02(X3I“)  +  hvSch_R  -  0(3P)  +  0(1D)  +  K.E.  (4.2) 

As  produced  by  solar  radiation,  the  approximate  energy  balance  of 
this  process  is  as  follows 
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wavelength  of  maximum  absorption  is  1500  A, 

i.e.,  hv  =  8.26  ev  (4.3a) 

energy  to  dissociation  G?(X3£*)  into  2  0(3P)  =  5.12  ev  (4.3b) 

^  y 

energy  to  excite  0(3P)  to  0(^D)  =  1.96  ev  (4.3c) 

balance  of  energy,  which  appears  as  kinetic 

energy  =  1.18  ev  (4.3d) 

The  energy  of  electronic  excitation  eventually  appears  a  kinetic 
energy  as  a  result  of  quenching  processes : 

0(1D)  +  N2  -  0(3P)  +  N2(vib)  +  K.E.  (4.4) 

Since  one  electron  volt  (kT  =  1  ev)  corresponds  to  a  temperature 
T  =  11,600°  K,  and  since  about  1%  of  the  atomic  oxygen  recombines 
over  night  to  be  re-dissociated  in  the  morning,  the  process  is  the 
single  most  substantial  source  of  heating  of  the  thermosphere.  It  is 
moreover  supplemented  by  kinetic  energy  of  macroscopic  mechanical  mo¬ 
tion  which  is  transported  upwards  in  the  atmosphere  by  "internal 
gravity  waves"  (Hines  et  al,  1965)  or  by  "internal  gravity  shocks" 
(Layzer,  1967). 

At  altitudes  above  150  km  the  photo -ionization  of  0  and  N2  con¬ 
tributes  significantly  to  the  photochemistry,  and  presumably  also  the. 
heating,  but  there  is  much  less  flux  of  solar  energy  at  the  higher 
energies  corresponding  to  wavelengths  less  than  11,200  A.  Also  there 

are  fewer  available  absorbing  atoms  and  molecules,  since  above  150  km 

—8  17 

there  only  remains  10  of  tne  earth’s  atmosphere,  or  10  cm 

column.  Thus  the  rise  in  temperature  above  150  km  is  slow. 

The  scale  height  H  of  Eq.  4.1  depends  on  molecular  weight  as 
well  as  on  temperature.  In  the  region  above  h  =  100  km  where  the  air 
is  no  longer  uniformly  mixed  there  will  be  a  diffusive  separation  of 
the  different  species,  with  the  lighter  species  (He,  0)  at  the  higher 
altitudes.  This  effect  supplements  the  photochemistry  to  cause  a 
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reduction  of  net  molecular  weight  with  increase  in  altitude  from  a 
value  29  at  80  km  to  a  value  18  at  500  km. 


Above  about  200  km,  the  fraction  of  free  electrons  becomes  ap¬ 
preciable  (Delcroix,  1965) 


number  density  of  free  electrons 
total  number  density 


>  10  ;  h  >  200  km  (4.5) 


under  these  conditions  the  motion  of  the  atmosphere  is  determined 
largely  by  the  motion  of  the  ionized,  rather  than  the  neutral,  con¬ 
stituents.  In  other  words,  the  geomagnetic  field  tends  to  become  a 
more  important  factor  than  the  prevailing  winds.  In  the  high  alti¬ 
tude  region  (h  >  200  km),  plasma  effects  dominate  over  the  effects  of 
continuum  fluid  dynamics. 

At  altitudes  above  80  km,  where  there  is  a  significant  degree  of 
dissociation,  the  partition  of  the  energy  of  the  atmosphere  in  dif¬ 
ferent  degrees  of  freedom  differs  drastically  both  from  that  of 
equilibrium  air  at  room  temperature,  and  that  of  equilibrium  air  at 
the  more  representative  temperature  of  1500°  K.  As  a  result  of  solar 
radiation,  more  energy  is  available  as  energy  of  dissociation  of  the 
C>2  molecules.  There  is  also  more  ionization  and  more  bound-state 
electronic  excitation  than  is  characteristic  of  local  thermodynamic 
equilibrium  (LTE). 

In  Fig.  58  are  shown  the  fractions  of  principal  neutral  species. 
The  data  of  Fig.  58  derived  from  CIRA-I965,  show  the  change  of  rela¬ 
tive  proportions  due  to  diffusion  without  other  mixing  above  about 
100  km.  In  Fig.  59  is  shown  the  variability  in  electron  density  en¬ 
countered  in  a  number  of  measurements. 

The  breakdown  of  energy  in  different  degrees  of  freedom,  based 
on  the  standard  temperature  and  the  composition  data  of  Figs.  59  and 
60,  is  shown  in  Fig.  60. 
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f  (  x )  =  n(  x )/ 


THI5  IS  A  MEASURE  OF  THE  DISCREPANCY 
BETWEEN  US  1962  AND  CIRA-1965 

PT  A  klfNADh 


tot 


tot 


US  62 
ClRA" 


SH 1-20-68-26 


FIGURE58.  Upper  Atmosphere  Chemistry:  Principal  Neutral  Species 
(Bauer,  to  be  published) 


for  the  electronic  energy  there  are  three  curves  in  Fig.  60:  a 
theoretical  curve  representing  LTE  and  constituting  a  plausible  lower 
bound,  and  two  other  curves  deduced  from  experimental  airglow  re¬ 
sults.  Specifically,  there  are  only  three  electronically  excited 
states  that  have  to  be  considered,  the  states  (a1^)  and  possibly 

(b^E*)  of  0o  which  are  respectively  0.98  and  1.60  ev  above  the  ground 
y  3  i 

state  (X  £  ),  and  the  state  0(  D)  which  is  1.96  ev  above  the  ground 
3  9 

state  0(  P).  Since  all  other  electronic  states  of  0,  C>2,  NT2,  NO, 
etc.  are  more  than  4  ev  above  the  respective  ground  states,  their 
overall  occupation  is  extremely  low  because  of  the  unfavorable 
Boltzmann  factor  exp  (-  EAT).  However,  they  may  be  important  for 
specific  processes. 

Thus  of  the  three  curves  for  electronic  energy  in  Fig.  59,  the 
experimental  data  refer  to  airglow  on  CL^a^A  )  and  0(^D) ,  while  the 
theoretical  LTE  curve  is  mainly  due  to  02(  a  Ag) .  The  experimental 
data  should  not  be  extended  in  altitude  range. 

Walker  (1968)  alleges  that  the  free  electrons  have  a  temperature 
larger  than  the  neutral  species  by  a  factor  of  two  at  h  =  10C  -  140 
km.  Walker  accounts  for  the  higher  temperature  to  be  a  result  of  the 
Schumann -Runge  absorption,  described  by  Eq.  4.2  by  which  a  large 
amount  of  0(^1))  is  produced  in  this  altitude  range.  The  0(^0)  trans¬ 
fers  energy  efficiently  into  the  vibration  of  N2  by  the  quenching 
reaction : 

-  OC^D)  +  N2  -  0(3F)  +  N2(v)  +  K.E.  (4.6) 

Despite  the  spin  change,  Zipf  (1968)  has  given  evidence  that  the  rare 
constant  of  process  (4.6)  is  perhaps  10-11  cm3/sec,  and  that  possibly 
half  of  the  energy  of  process  (4.6),  or  1  ev  transfers  into  the  vi¬ 
bration  of  Ng.  Since  the  interchange  between  kinetic  energy  of  free 
electrons  and  vibrational  energy  of  N2  is  rapid  (Schulz,  1962); 

Walker  (1968)  interprets  the  observed  enhanced  kinetic  energy  of  the 
electrons  in  terms  of  a  vibrational  temperature  of  N2  of  3000°  K. 


Walker’s  result  is  included  in  Fig.  60,  but  should  not  be  regarded  as 
generally  accepted . 

Also  shown  in  Fig.  60,  the  kinetic  energy  of  a  100  m  sec  1  wind, 
typical  of  these  altitudes,  is  a  minor  contributor  to  the  overall 
kinetic  energy. 

Until  recently,  the  terrestrial  atmosphere  was  thought  to  be 
mixed  to  a  uni.  orm  composition  to  a  uniform  height  of  approximately 
100  km.  In  the  altitude  region  80-100  km  the  uniform  mixing  was 
thought  to  be  affected  by  atmospheric  turbulence  which  terminated 
abruptly  at  a  "turbopause"  at  approximately  100  km.  Above  the  turbo¬ 
pause  the  various  chemical  species  were  supposed  to  diffuse  sepa¬ 
rately,  each  with  its  own  scale  height 

Hj  =  kT/M j g ;  j  =  C,  02,  NO,  Ng,  Ar,  He,....  (4.7) 

Experimental  work  cited  in  Table  11  gives  evidence  that  critical 
levels  in  the  upper  atmosphere  such  as  the  "turbopause,"  referring  to 
"cessation  of  turbulence,"  or  "homopause,"  referring  to  "cessation  of 
uniformity,"  may  vary  by  as  much  as  20  km  in  altitude  depending  on 
latitude,  season  and  possibly  even  time  of  day.  Evidence  from  a 
number  of  sources  indicates  the  existence  of  stratification  below 
the  turbopause;  e.g.,  sporadic-E,  wind  shears  and  different  ion  struc¬ 
ture  on  a  vertical  scale  of  order  1-3  km  (Kochanski,  1964).  Evidence 
for  compositional  changes  with  altitude  near  the  homopause  has  been 
reviewed  by  K.  S.  W.  Champion  (1967,  1968,  1969). 

The  evidence  concerning  the  state  of  motion  of  the  lower  thermo¬ 
sphere  may  be  summarized  as  follows: 

Below  80  km,  the  atmosphere  is  uniformly  mixed  at  all  times, 
while  above  120  km  it  is  diffusively  separated,  as  evidenced  by 
measurements  of  the  Ar-N2  concentration  ratio  (Kasprzak  et  al.,  1968). 
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GEOPHYSICAL  SERIES,  NO.  11  NOVEMBER  1960,  pp.  1649-1657,  AGU  TRANSLATION  BY  R.  MUDGE,  pp.  1099-1105. 

J.  E.  BLAMONT,  ET  AL,  ATMOSPHERIC  MOTIONS  STUDIED  BY  SODIUM  VAPOR  TRAILS  RELEASED  BY  ROCKETS,  COSPAR  INFORMATION  BULLETIN  13, 
p.  10,  JANUARY  1963. 


Over  all  altitudes  up  to  perhaps  200  km  there  is  evidence  for 
horizontal  winds.  In  the  region  60-160  km  the  magnitude  is  typically 
up  to  90  m/sec,  while  the  direction  varies,  according  to  Kochanski 
(1964),  whose  data  are  given  in  Fig.  50  in  the  discussion  on  the 
stratosphere  and  mesosphere,  and  Bedinger  and  Constantinides  (1969). 

Bedinger  and  Constantinides  (1969)  have  statistically  analyzed 
wind  data  derived  from  57  chemical  release  trails  over  Wallops  Island, 
Virginia,  released  in  evening  twilight  (23),  morning  twilight  (19), 
and  otherwise  during  the  night  (15).  The  height  interval  covered  by 
the  trails  ranged  from  63  km  to  220  km,  of  which  range  data  from  the 
heights  80  to  190  km  was  statistically  analyzed.  Sample  size  in  the 
interval  90  to  150  km  exceeds  30;  in  the  interval  90  to  130  km  sample 
size  exceeds  50.  The  seasonal  distribution  of  the  data  is  the  follow¬ 
ing:  Spring,  10  trails  on  7  occasions;  Summer,  20  trails  on  6 
occasions;  Fall,  11  trails  on  7  occasions;  and  Winter,  16  trails  on  5 
occasions . 

Figure  61  shows  the  average  wind  speed  in  the  range  of  altitude 
90  to  150  km,  ranging  in  value  from  less  than  50  to  above  90  m  sec 
Figure  62  shows  the  North  and  East  components  of  the  average  wind  and 
their  standard  deviations.  The  ratio  of  standard  de/iation  to  mean 
varies  from  1  to  2-  Figure  63  shows  the  average  value  of  the  vertical 
wind  shear  and  its  North  and  East  components. 

There  is  evidence  of  "internal  gravity  waves"  of  a  nature 
described  in  the  section  on  the  stratosphere  and  mesosphere.  In  the 
thermosphere,  the  gravity  wave  amplitude  increases  from  60  km  up  to 
about  180  km,  and  thdVi  falls  off,  as  describee  by  Kochanski  (1966), 
whose  results  are  sho'vr.  in  Fig.  64. 

From  rocket-borne  mass  spectrometer  results  (Narcisi  and  Bailey, 
1965)  and  ionosonde  measurements  of  sporadic-E,  there  is  clear  evidence 
of  horizontal  structure  in  the  80-100  km  region,  typically  with 
vertical  scale  of  1  to  3  km,  lifetime  on  the  order  of  one  hour  and 
horizontal  extent  of  about  100  km. 
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FIGURE  63.  Average  Magnitude  and  Components  of  Shear 
(After  Bedinger  &  Constantinides,  1969) 
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FIGURE 64.  Magnitude  of  Drift  and  Internal  Gravity  Waves 
{After  Kochanski,  1966) 

Spreading  of  rocket  release  vapor  trails  of  sodium,  tri-methyl 
aluminum  (TMA),  and  other  materials  show  a  marked  change  in  effective 
diffusion  coefficient  ac  altitudes  in  the  range  100-120  km.  The 
effective  diffusion  coefficient  at  the  lower  altitudes,  significantly 
larger  than  the  molecular  coefficient,  has  been  called  an  "eddy 
diffusivity . " 

The  probing  of  the  atmospheric  turbulence  by  chemical  releases 
poses  two  questions:  (1)  the  passage  of  rockets  or  chemical  dispensers 
unavoidably  induces  disturbances  which  should  be  eliminated  in  the 
interpretation  of  data,  and  (2)  the  diffusion  of  the  passive  conta¬ 
minant  by  the  turbulent  motion  in  the  atmosphere  does  reflect  the 
structure  of  the  atmospheric  turbulence,  but  the  two  motions  are 
obviously  not  identical.  Their  relation  should  be  investigated. 
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Radar  returns  from  the  ionization  trails  of  meteors  in  the  same 
altitude  regime  show  the  same  effect  as  an  onset  of  turbulence.  Be¬ 
cause  of  the  smaller  size,  the  inducing  of  turbulence  by  the  passage 
of  the  meteor  is  less  likely  than  in  the  case  of  rocket- releases. 

Experimental  data  shown  in  Fig.  51  exhibit  values  of  effective 
diffusivity  which  are  larger  than  molecular  diffusivity  at  altitudes 
below  about  120  km.  Since  many  of  these  data  were  obtained  by  measure¬ 
ments  of  the  up-trail,  which  differ  substantially  from  measurements 
made  on  the  down-trail,  the  reality  of  the  implied  effective  or  "eddy" 
diffusion  coefficients  shown  in  Fig.  51  are  controversial. 

Explanation  of  observational  data  should  give  reasons  for  a 
well-defined  turbopause  at  altitudes  of  about  100  km;  below  which  the 
atmosphere  is  turbulent  and  above  which  it  is  laminar.  Turbulence 
arises  when  a  fluid  is  subjected  to  a  srrain  too  severe  to  be  re¬ 
lieved  by  molecular  transport  in  the  form  of  a  velocity,  temperature, 
or  concentration  gradient.  Strong  destabilizing  concentration 
gradients  are  uncommon  in  the  atmosphere.  The  atomic  oxygen  gradient 
arising  in  the  0-0  2  cycle  described  above  is  in  fact  stabilizing,  not 
destabilizing,  so  that  it  will  not  induce  turbulence.  However, 
velocity  gradients  (wind  shear)  and  temperature  gradients  are  of  com¬ 
mon  occurrence  in  the  atmosphere. 

The  criteria  for  producrion  of  an  instability  leading  to  tur¬ 
bulence  by  velocity  and  temperature  gradients  may  be  described  in 
terms  of  a  Reynolds  number  Re  and  a  Richardson  number  Ri.  Instability 
due  to  wind  shear  is  implied  by  Reynolds  number,  in  terms  of  a 
transverse  coordinate. 

Re  >  Recrit  ~  3000  (4.8) 

wherein  Re  is  a  critical  value  of  Reynolds  number 

Instability  in  the  presence  of  a  density  or  temperature  gradient 
under  gravity  is  implied  by  Richardson  number: 

Ri  <  Ri  ~  0.04  -  1  (4.9) 

criu 

wherein  Ri  ..  is  a  critical  value  of  Richardson  number. 
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The  Reynolds  number  is  defined  (Schlichting,  1960,  p.13)  as 


follows 


inertial  force  _  p  u  5  u/Bx 


frictional  force 


U  3  u/3z" 


P  V  d/|i 


(4.10) 


where  u  =  local  velocity  (horizontal,  in  x-direction) 

z  =  vertical  coordinate  (perpendicular  to  the  flow) 
p  =  density  of  fluid 
(j  =  viscosity  of  fluid 

d  =  characteristic  dimension  of  fluid  motion 
V  =  free  stream  velocity  such  that  gradient  is  ~  V/d. 

The  Richardson  numtSer  Ri  is  defined  (cf.  Goody,  1964,  p.  367) 


Ri‘  =  (g/T)  (BT/Bz  +  g/cp) /( Bu/3z)‘ 


(4.11) 


Under  representative  conditions  for  typical  horizontal  wind 

-2  -1  -1 

speeds  100  m/sec,  and  wind  shears  of  10  -10  sec  ,the  effective 

Reynolds  number  is  greater  than  the  critical  value  Re^^  for  al¬ 
titudes  below  80-110  km,  implying  a  tendency  for  turbulence.  At 
higher  altitudes  the  Reynold’s  number  is  generally  much  lass  than 

Recrit ' 

In  the  stratosphere  and  mesosphere,  between  altitudes  20  and 
about  80  km  the  Richardson  number  is  always  far  greater  than  the 
critical  value,  implying  strong  stabilization  against  turbulence,  as 
shown  in  Fig.  46-  However,  because  the  positive  temperature  gradient 
in  the  lower  thermosphere  is  stabilizing,  not  destabilizing,  the 
tendency  of  the  lower  thermosphere  to  turbulence  is  marginal  below 
about  110  km.  Relatively  small  disturbances  in  the  lower  atmosphere, 
present  because  of  meteorological, , tidal,  and  diurnal  variation,  may 
be  amplified  in  propagating  upwards,  ir.  the  form  of  "internal  gravity 
waves"  (Eckart,  1960,  ahd  Hines  et  al.,  1965),  or  even  lead  to  the 
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thought  of  "in.ernal  gravity  shocks"  (Layzer,  1967).  Such  an  amnJi- 
fication  may  be  reasoned  in  an  elementary  fashion  by  the  conservation 
of  momentum 


2 

p  u  =  constant 

so  that,  since  density  decreases  with  altitude 


(4.12a) 


P 


~  a“8/H 


(4.12b) 


':he  velocity  will  consequently  increase  with  altitude,  according  to: 


u 


e+*/2H 


-4.12c) 


PLASMA  AND  HYDROMAGNETIC  EFFECTS 


The  significance  of  the  ionization  in  the  upper  atmosphere  is 
evidenced  by  the  reflection  of  radio  waves.  The  phenomenon  is  techno¬ 
logically  exploited  for  purposes  or  radio  communication. 

The  mapping  of  the  electric  and  magnetic  fields,  together  with 
the  distributions  of  ions  and  electrons  in  the  upper  atmosphere  have 
been  recently  the  subject  of  active  experimental  studies  of  ionized 
metallic  vapor  clouds. 


A  parameter  which  described  tne  relative  importance  of  collisions 
between  electrons,  charged  heavy  particles  and  neutral  particles  is 
the  collision  frequency.  Denoting  the  collision  frequency  v..  as  the 
frequency  of  collisions  between  particles  of  species  i  and  j,  where  i 
and  j  may  be  called  e  for  electrons,  i  for  ions,  and  n  for  neutral 
particles,  one  may  write  usinc  relations  of  Hochstim  (1969),  and 
Ratcliffe  and  Weekes  (1960) 


ei 


=  3.633 


1372 


3279 


In 1 1  +  ■ 


rrV2 

xe 


[(Ne(VTl)] 


.  6,  uor^)-3"  .. 


0 


.13) 
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vfen  =  2.33  (iO)'11  Nn  [1-1.21(10)"4  Te]  «  4.4  (10)"9 


(for  N2  as  neutral  gas) 


(555)  K„ 


(4.14) 


_q  N  +  N 

v^n  =  2.6  (IO)-  --n--  - —  by  Chapman  according  to  Hochstim  (1969)  (4.15) 


M 

where  M  is  the  mean  molecular  weight.  N  ,  N- ,  and  N  are  electron  and 
,  3  e  1  n 

ion  densities  (cn\  ),  and  T  =  =  T£  is  the  temperature  in  degrees 

Kelvin. 

Tltie  ratio  v  .j/v  is  plotted  in  Fig.  65  as  a  function  of  altitude, 
using  expressions  of  Chapman  (Hines  et  al.,  1965)  and  the  electron 
density  range  of  Fig.  64.  Hie  ratio  is  larger  than  unity  for  altitudes 
above  210  to  270  km.  Below  about  150  km,  in  the  D-  and  E- regions,  the 
electron  to  neutral  collisions  clearly  dominate.  Above  altitudes  from 
150  to  200  km,  the  electron  to  ion  collisions  must  also  be  taken  into 
account . 

Martyn  (1959)  estimates  the  value  of  v  .,  v  and  v.  to  be 

ei  en,  xn 

respectively  880,  2300,  and  9G  sec-J  at  the  FI  ledge,  and  900,  40,  and 
1  sec-1  at  the  F2  peak.  Thus  in  the  FI  ledge,  electron  collisions  are 
mainly  with  neutral  particles,  and  in  the  F2  iayer,  mainly  rfith  ions. 
Hochstim' s  (1365)  computations  corroborate  Martyn  in  a  general  way. 

Inthe  upper  F  region,  in  which  electron  collisions  are  mainly 
-3/2 

with  ions,  since  T  '  ,  as  temperature  increases  during  the  day, 

the  collision  frequency  decreases.  Conversely,  since  v  «  T,  in  the 
D-  and  E-regians,  wherein  electrons  collide  chiefly  with  neutral 
particles,  the  collision  frequency  will  increase  with  increasing 
temperature  during  the  day. 

In  the  case  of  a  largo  ratio  of  gyro  frequency  to  collision 
frequency,  the  ions  and  electrons  are  constrained  to  move  along  rather 
than  across  the  magnetic  field  lines.  The  gyro  frequencies  are 

electron  gyro  frequency  u)ge  =  cB/mc  (4.16a) 


ion  gyro  frequency  =  eB/K^c 


v4.16b) 


L7,3 


-v- 


The  ratios  «»Be/ven,  «iBi/vin  are  plotted  in  Fig.  66.  Figure  67  shows 
plot  of  theoretical  values  (solid  line)  of  collision  frequency  and 
experimental  determinations,  quoted  from  Hines  et  al.  (1965).  In 
terms  of  gyro  frequency  criteria,  the  electrons  follow  the  field  for 
h  >  75  km  while  the  ions  follow  the  field  only  for  h  >  125  km.  In 


FIGURE 65.  Plasma  Effects  (After  Bauer,  1969) 
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COLUSIONAL  FREQUENCY,  v,  sec 


FIGURE  67.  Collision  Frequency  y  of  Electrons 
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the  intermediate  region  (75  to  125  km)  the  electrons  and  ions  are  con¬ 
strained  by  different  forces,  suggesting  the  possibility  of  net  cur¬ 
rents,  which  are  observed  to  occur.  Such  occurrence  is  the  "dynamo 
region"  discussed  by  Hines  et  al.  (1965). 

One  energy  input  into  the  earth's  atmosphere  from  the  sun,  in 
particular  in  the  case  of  solar  storms,  flares  and  other  disturbances, 
comes  in  the  form  of  charged  particles.  The  bulk  of  the  energy  is 
carried  by  10  Kev  protons  arid  low-energy  electrons  having  roughly  the 
same  velocity  as  the  protons.  The  protons  moving  along  geomagnetic 
field  lines  enter  the  earth's  atmosphere  in  the  auroral  zones  of  the 
polar  regions  where  they  excite  auroras  and  heating,  leading  to 
significant  density  variations. 

A  different  but  important  effect  is  the  reflection  of  radio 
waves  from  the  ionosphere.  Even  a  small  electron  density  can  cause 
an  important  effect.  The  plasma  frequency  as  a  function  of  altitude 

u>  -  (n  e/m c)'2  (4.18) 

M  C 

is  plotted  in  Fig.  6  5.  Electromagnetic  signals  of  frequency  u)  <  ui^ 
are  reflected  on  normal  incidence  on  an  optically  thick  layer  of 
electron  density  ng,  or  plasma  frequency 

Effects  of  Coulomh  screening,  characterized  by  the  Debye  radius, 
given  ty 

hn  =  (kT/4nn  e2)^  =  6.9  (T/n  cm.  (4.17) 

V  S  Q 

are  not  important  in  the  thermosphere.  Typical  values  of  the  Debye 
radius  hp  are  of  the  order  0.3  cm,  which  is  less  than  the  mean  free 
path  above  80  km  and  yet  very  much  greater  than  the  Landau  radius 
rp  =  e2/kT. 

Time -dependent  phenomena  such  as  magneto-hydrodynamic  waves, 
discussed  by  Poeverlein  (1967),  may  have  significant  effects  on  the 
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overall  structure  and  motion  of  the  thermosphere  above  100-200  km,  at 
an  altitude  aependent  on  the  specific  phenomenon  under  investigation 
as  well  as  on  the  latitude  and  the  solar  activity. 

A  second  topic  under  which  the  upper  atmosphere  may  be  described 
and  its  regions  subdivided  is  that  of  chemistry.  The  solar  insolation 
especially  in  the  very  energetic  wavelengths  causes  ionization.  The 
chemical  reactions  which  take  place  in  the  stratosphere,  mesosphere, 
and  thermosphere  are  many  and  important  in  consideration  of  energy 
transformations . 

In  the  technologically  interesting  regions  of  ionization,  the 
chemistry  of  dissociation,  driven  by  the  solar  radiation,  and  not  the 
ionization  itself,  primarily  determines  the  constitution  and  charac¬ 
teristics  of  the  upper  atmosphere. 

The  species  that  can  be  ionized  may  be  arrayed  in  the  order  of 
ionization  potential  given  in  Table  12. 

The  wavelengths  corresponding  to  the  ionization  potentials  of 
all  atmospheric  constituents  listed  in  Table  12,  with  the  exception 
of  sodium  and  other  metals,  are  in  the  vacuum  ultraviolet.  The 
significant  photo ionization  occurs  above  the  sensible  atmosphere, 
dominantly  as  the  ionization  of  02  by  the  Schumann -Runge  continuum 
radiation  band.  Both  0 2  and  NO  may  be  ionized  by  Hydrogen  Lyman 
alpha  which  goes  through  to  60-80  km. 

Ionization  is  produced  by  solar  ultraviolet  and  to  a  lesser  ex¬ 
tent  by  solar  X-rays.  Johnson  (1961)  has  described  the  available 
energy  to  be: 


Ultraviolet 

(1350  -  800  A  | 
j9 ,25  -  14;S  ev  J 


X-ray,  at  solar  maximum 

50A  ±  factor  2  8  7 

250  ev  Total  flux  -  2.5  x  10“ 
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TSBLE  12 .  IONIZATION  POTENTIAL  OF  ATMOSPHERIC  CONSTITUENTS 


_ — — — - - - - - 

First  ionization 
Potential  (ev) 

Wavelength , 

A 

Comments 

Na 

5.12 

2420 

Metals  have  very  low  IP’s, 
i . e .  even  small  concentra¬ 
tions  can  be  major  sources 
of  ionization  at  lower  al¬ 
titudes 

NO 

9.25 

1340 

NO  is  a  minor  species  but^ 
with  a  very  low  IP :  '  is 

important  for  70-120  .or.. 

°2 

12.05 

1030 

C>2  is  important  for  100- 
160  km. 

n 

15.6 

910 

0+  is  important  from  150- 
500  km. 

H 

13.6 

9 1C 

H+  is  important  above  500 
km. 

N2 

15.5 

800 

IP  is  so  high  that  Ng 
is  not  important 

N 

14.55 

850 

N  isn’t  common  (Np  isn’t 
photodissociated  much  by 
solar  UV) ,  and  it  has  a 
high  IP . 

Ar 

15.7 

790 

Too  high  IP  to  be  impor¬ 
tant. 

l!_ 

24.5 

510 

Too  high  IP  to  be  impor¬ 
tant. 

The  ions  produced  at  various  altitudes,  according  to  Bortner  & 
Kummler  (1968),  Narcisi  &  Bailey  (1965),  Johnson  (1965),  are  listed 
below : 

H+ 

0+ 


N0+ 

Mg+,  SiO+,  H30+,  H502+,  etc 

The  chemical  reactions  involving  the  principal  constituents  of 
the  lower  thermosphere  are  complicated.  They  include  the  following 
processes : 

(a)  photo-ionization: 


hv  +  02  -  02  +  e” 

(4.18) 

hv  +  G  —  0^  +  e 

(4.19) 

hv  +  K2  -*  N2  +  e 

(4.20) 

charge  transfer:  (Exothermic  charge 
rapid) 

transfer  processes  are 

0+  +  NO  -  N0+  +  0 

(4.21) 

N2  +  °2  °2  +  N2 

(4.22) 

(c)  recombination: 


300 

km 

140 

km 

120 

km 

90 

km 

80 

km 

0  +  e 


0  4-  hv  radiative  -  slow 


(4.23) 


0+  +  e~  +  M  -  0  +  M  3-body  -  slow  at  (4.24) 

high  altitudes 
and  low  density 

NO+  +  e-  -  N’  +  0”  dissociative  -  rapid  (4.25) 

Although  the  air  constituent  NO  is  minor,  it  is  metastable  chemically. 
Also  because  NO  has  the  lowest  ionization  potential  of  any  of  the  air 
species,  N0+  accumulates  as  the  end-result  of  a  chain  of  reactions. 

At  the  higher  altitudes  of  the  lower  thermosphere  there  is  no 
NO,  lots  of  N2,  0,  and  some  Oj.  N^  charge  transfers  to  0^  and  0+. 
Since  0+  can  not  recombine  as  fast  as  it  accumulates. 

At  the  lower  altitudes  the  chemistry  is  complicated  (Bortner  and 
Kummler.  1968). 

NO  is  photo-ionized  by  Hydrogen  Lyman  alpha  (Craig,  1965). 

Mg,  SiO,  deposited  as  a  result  of  meteor  ablation  around  90-100 
km  are  easily  ionized. 

Hydrated  water  ions  H+  (H20)p»  n  =  1,  2,...  are  present. 

Because  of  the  relatively  high  densities  of  the  lower  altitudes, 
significant  numbers  of  negative  ions,  0~,  CO^,  N0~,  especially  at 
night,  are  formed  by  dissociative  attachment  processes  of  type 

e  +  AB  -*  A  +  B 

Excited  states  due  to  chemi-ionization  are  also  present. 

The  several  subregions  of  the  mesosphere  and  lower  thermosphere 
have  been  categorized  according  to  the  degree  of  ionization  by  geo¬ 
physicists  interested  in  the  propagation  of  radio  waves,  as  the  C-, 

D- ,  E-  and  F-regions.  The  C-region  is  in  one  lower  mesosphere;  the 
D-region  is  within  the  upper  mesosphere;  the  E-region  in  the  lower 
thermosphere,  and  the  F-region  is  in  the  upper  thermosphere.  The 
ionospheric  F-region  has  a  lower  limit  at  the  rop  of  the  E-region  at 
a  height  of  about  2^0  km  and  no  generally  accepted  upper  height  limit. 
Each  of  these  regions  is  distinguished  from  the  others  in  its  chemical 
nature . 
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The  bulk  of  the  atmosphere  in  the  F-region  is  neutral  gas.  Since 
mixing  processes  other  than  diffusion  are  of  little  significance,  the 
composition  of  the  F-region  is  determined  by  ionization  production 
and  loss  processes  and  by  diffusive  separation  of  the  chemical  con¬ 
stituents.  In  the  lower  F-region,  the  principal  constituents  are 
molecular  nitrogen  dominating  atomic  oxygen.  With  increasing  alti¬ 
tude,  atomic  hydrogen  becomes  the  dominant  constituent.  The  scale 
height  at  the  base  of  the  F-region  is  about  30  km. 

The  formation  of  the  F- layer  below  600  km  is  governed  by  the 
continuity  equation 

||  +  div  (Nv)  =  q  -  L(N)  (4.26) 

wherein  N  =  electron  density 

v  =  local  mean  electron  velocity 
q  =  production  race  of  electron 
L  =  loss  rate  of  electrons  and 

I»(N)  =  a  N2  +  e  N  (4.27) 

wherein  a  =  recombination  coefficient 
e  =  attachment  coefficient 

For  equal  densities,  the  flux  of  ionization  transport  becomes 
(Hines  et  al. ,  1965,  Eq.  4.3) 

div  (Nv)  =  div  ^N(vgm  +  vT  +  vp)J  (4.28) 

by  which  three  types  of  movement  are  recognized  explicitly,  i.e., 
movement  caused  by  electromagnetic  forces,  by  temperature  changes  of 
the  atmosphere,  and  by  plasma  diffusion.  The  solution  of  the  con¬ 
tinuity  equation  is  complicated  if  all  three  types  of  ionization 
movement  are  simultaneously  taken  into  account. 
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Below  about  600  km,  the  F-region’s  neutral  atmosphere  consists 
primarily  of  0  and  Ng >  with  09>  N,  H  and  He  as  minor  constituents.  The 
type  of  electron  loss  depends  on  the  positive  ions  present.  Elec¬ 
trons  may  combine  directly  with  atomic  ions  or  dissociative] y  recom¬ 
bine  with  Og,  Ng  and  N0+ .  The  combination  with  atomic  ions  is  sc 
slow  that  below  the  Fg  peak  the  rate  is  negligible  compared  to  that 
for  dissociative  recombination. 

The  most  important  loss  reactions  in  the  F-region  are  those  for 
0  (Hines  et  al,  1965,  Eq.  4.4  and  4.5) 


0+  +  00  -  0+  +  0 


(4.29a) 


Oj  +  e  -  0*  +  0** 


(4.29b) 


and  those  for  N0 


0+  +  Ng  -  NO+  +  N 


(4.30a) 


NO  +  e  -  N*  +  0 


(4.30b) 


The  rate  of  the  general  reaction  of  (4.29a)  and  (4.30a)  is  (Hines 
et  aL,  1965,  Eq.  4.6) 


=  -  Kx  n(0+)  n(XY) 


(4.31) 


The  rate  of  the  general  reaction  of  (4.29b)  and  (4.30b)  is 


XT  =  -  Kg  Ne  "(XY+) 


(4.32) 


In  (4.31)  and  (4.32),  and  Kg  are  constant  rate  coefficients;  Ng 
is  the  electron ■ member  density;  n(XY)  is  the  density  of  neutral  di¬ 
atomic  molecules  (Og  or  Ng) ;  n(XY'*')  is  the  density  of  positive  di¬ 
atomic  ions  (0g  or  N0+).  The  attachment  coefficient  (6)  of  (4.27) 
equals 


describe  it  for  the  stratosphere  and  r^esosphere.  Evidently  this 
high  concentration  arises  from  low- Level  photochemistry  of  a  differ¬ 
ent  nature,  such  as  that  generating  smog,  etc. 

To  summarize,  at  altitudes  above  50  km  |(h)  <  1  because  there 
is  too  little  ozone  ro  make  the  atmosphere  optically  thick;  at 
altitudes  below  50  km,  since  most  of  the  solar  radiation  incident 
on  the  atmosphere  has  already  been  absorbed  at  higher  altitudes 
(|  >  1),  there  is  inadequate  solar  radiation  in  the  Hartley  band  to 
heat  up  the  atmosphere. 

Having  explained  the  temperature  maximum  at  the  stratopause, 
the  lower  bound  of  the  stratosphere  and  the  upper  bound  of  the 
mesosphere  will  next  be  described. 

The  lower  bound  of  the  stratosphere  is  given  by  the  tropopause, 
which  corresponds  to  the  temperature  minimum  between  the  tropopause, 
or  that  portion  of  the  atmosphere  which  is  in  direct  "contact"  with 
the  surface  of  the  earth,  and  the  thermopause. 

If  the  temperature  gradient  (with  increasing  height)  within 
the  troposphere  is  negative,  so  that  the  adiabatic  lapse  rate 

-dT/dh  <  g/Cp  (3.10) 

(where  g  =  acceleration  due  to  gravity  and  cp  =  specific  heat  at 
constant  pressure),  there  is  possibility  of  a  fluid  mechanical 
static  instability  by  which  a  parcel  of  air  of  slightly  different 
temperature  can  move  a  significant  vertical  distance,  leading  to 
overturning  of  the  atmosphere.  Such  an  overturning  in  the  tropo¬ 
sphere  is  particularly  visible  in  the  case  of  cumulus  clouds. 

The  vertical  motion  arising  from  static  instability  is  of 
vital  importance  in  keeping  water  droplets  or  ice  crystals  in 
suspension  in  the  troposphere.  Since  the  terminal  velocity  of  a 

spherical  water  droplet  of  Ip.  radius  predicted  by  the  Stokes' 

_2 

formula  is  about  10  cm/sec,  modest  mean  vertical  motions  may  keep 
such  droplets  in  suspension. 


8  =  Kx  n(XY) 


(4.33) 


and  the  recombination  coefficient  (a)  of  (4.27)  is 


a  -  K, 


(4.34) 


The  net  loss  rate  for  che  reactions  4.29  and  4.30  is  (Hines  et  al. . 
1965,  Eq.  4.6c)  ’ 


a  B  M 
8  +  a  N 


(4.35) 


where  a  is  constant  and  0  is  height  dependent.  If  $  »  a  N)  as  in 
the  lower  F- region,  then 


L(N)  s:  a  N  ;  0  »  a  N 

If  8  «  a  N,  as  in  the  upper  F- region ,  then 


(4.36) 


L(N)  -  6  N;;P  «  aN 

'high 


(4.37) 


and  the  net  loss  is  controlled  by  an  attachment-like  process. 

Although  an  appreciable  fraction  of  the  solar  ultraviolet  radia¬ 
tion  is  consumed  in  producing  N*  ions,  mass  spectrometry  indicates 
N2  to  be  relatively  scarce  in  the  entire  F-region.  This  scarcity  of 
N2  ions  implies  that  either  the  dissociative  recombination  of  N+  is 
so  rapid  that  it  makes  no  contribution  to  observed  ionization,  or 
that  the  reaction  (Hines  et  al,  1965,  Eq.  4.9) 


N+  +  0  ->  N0+  +  N 


(4.38) 


preferentially  prevents  the  rapid  removal  of  all  electrons,  thereby 
making  a  substantial  contribution  to  F-region  ionization. 
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Measurements  of  the  F  region  made  at  night,  when  q  =  0,  suggest 
that  the  attachment  coefficient  has  the  empirical  value 

,300-lu 
-T  _4  '  qn  ' 

6  (sec  )  =  10  e  ^  (4.39) 

It  is  apparent  that  loss  processes  for  day  and  night  in  the  F-region 
are  different,  possibly  due  to  a  temperature  dependence  of  the  rate 
coefficients,  to  a  progressive  change  in  abundance  of  constituents, 
or  to  the  different  height  distributions  of  the  principal  positive 
ions. 

The  temperature  of  the  F-region  is  not  accurately  known.  The 
electron  temperature  Tg  and  the  neutral  gas  temperature  Tn  are  prob¬ 
ably  equal  at  night,  but  unequal  during  the  day. 

In  the  ionized  E-region,  which  ranges  in  height  from  an  upper 
limit  at  about  140  km  to  a  lower  limit  at  about  90  km,  a  well  defined 
layer  of  ionization  is  found  with  its  peak  near  100  km. 

Above  about  100  km,  mixing  processes  characteristic  of  the  lower 
levels  appear  to  quench,  with  the  consequence  that  diffusive  separa¬ 
tion  of  the  chemical  constituents  takes  place.  Also  at  90  to  100  km, 
there  is  an  important  chemical  change  resulting  from  the  photo  disso¬ 
ciation  of  C>2  into  0,  which  becomes  complete  at  about  160  km. 

The  ionizing  radiations  responsible  for  the  E-layer  must  be  of 
such  nature  that  they  penetrate  the  regions  above  without  serious  loss, 

then  are  absorbed  within  the  E -region.  Early  attention  was  drawn  to 

the  possibility  that  the  pertinent  radiation  may  be  specific  to  the 

ionization  of  O^j  a  molecule  known  to  be  virtually  absent  above  about 

160  km.  Another  source  would  be  x-rays  of  wave  length  10-D00A, 

shorter  than  the  broadband  100-1200A  for  whicn  Ng  has  appreciable 

cross  sections .  The  ionization  produced  n!  seems  to  have  a  large 

—6  —7  “1 

recombination  coefficient  of  10  or  10  cm  sec  ,  whose  effect  is  to 
remove  ions  nearly  as  rapidly  as  they  are  produced. 
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The  dominant  electron  loss  process  in  the  E-region  is  one  of 
dissociative  recombination  with  rate  coefficient  aQ  suggested  by 
Bates  and  Massey  (1947).  The  reaction  is  of  the  form 

( XY)+  +  e  -  (XY)*  -  X*  +  Y*  (4.40) 

in  which  a  positive  diatomic  species  (such  as  Ng ,  Og  or  N0+  in  the  B- 
and  E-regions)  combines  with  an  electron  to  form  one  or  two  excited 
atoms . 

The  estimates  of  dissociative  recombination  coefficient  are 

ao  (N2>  =  3(10)~7  cm3  sec~1  (4.41a) 

ao  (Op  =  1.7(10)~7  cm3  sec"1  (4.41b) 

aQ  (N0+)  =  2(10)"8  cm3  sec"1  (4.41c) 

The  collision  frequency  for  electrons  in  the  E-region  is  given  in 
Fig.  67. 

The  D-region  extends  from  the  lowest  level  of  atmospheric  ioniza¬ 
tion  (as  low  ao  50  km)  to  about  90  km.  It  includes  the  level  of  the 
mesopause,  at  80  to  SO  km  height,  where  the  temperature  reaches  a 
minimum  variously  estimated  from  130°  K  to  250°  K.  In  the  D-region 
and  below,  the  atmosphere  is  subject  to  convective  overturning  which 
leads  to  a  mixing  of  atmospheric  constituents .  As  a  result  of  the 
mixing,  the  relative  concentrations  of  most  chemical  species  are 
virtually  independent  of  height. 

In  addition  to  Mg  and  Og,  ir>  the  approximate  proportions  of  4:1, 
various  trace  constituents  play  roles  in  important  processes.  For 
example,  C09  and  0^  play  an  important  part  along  with  Cg  in  maintain¬ 
ing  radiative  energy  balance  and  so  in  controlling  the  temperature 
distribution.  The  ionization  balance  may  be  similarly  influenced  by 
the  trace  constituent  NO. 
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MOTION  OF  PLASMAS  IN  THE  UPPER  ATMOSPHERE 

The  plasma  ir.  the  upper  atmosphere  consists  of  electrons,  ions, 
and  neutral  part  ieles  which  interact  among  themselves  and  move  in 
certain  configurations  of  external  electric  and  magnetic  fields.  In 
addition,  induced  fields  are  also  excited  by  the  motion  of  the  ionized 
components . 

In  the  following  section  an  orbital  theory  of  a  single  charged 
particle  in  ar.  electric  field  and  a  magnetic  field  is  developed,  and 
some  elementary  features  of  the  motion,  together  with  the  characteristic 
parameters  involved,  are  reviewed. 

Since  the  orbital  theory  does  not  reproduce  the  collective  effects 
in  the  motion  of  an  assembly  of  particles,  it  is  necessary  to  develop  a 
hydrodynamic  theory  to  describe  the  motion  and  plasma  phenomena  in  the 
upper  atmosphere,  e.g.,  aurora,  airglow,  instability,  wave  propagation, 
etc.  A  hydrodynamic  theory  describing  the  motion  of  charged  particles 
requires  the  knowledge  of  transport  coefficients:  conductivity  and 
diffusion  tensors. 

In  a  later  section,  a  theory  of  transport  coefficients  for  both 
the  weakly  ionized  gas,  valid  at  altitude  up  to  200  km,  and  the 
strongly  ionized  gas,  valid  at  altitude  above  200  km  is  presented.  Not 
included  in  the  theory,  however,  is  the  effect  of  fluctuations  in  the 
structure  of  such  transport  properties  which  exhibit  anomalous  trans¬ 
port  phenomena.  Discussions  of  the  hydrodynamic  equations  for  a  plasma 
including  such  transport  properties  are  also  omitted. 

ELEMENTARY  FEATURES  OF  THE  ORBITAL  MOTION  OF  AN  IONIZED  GAS  IN  AN 
ELECTRIC  AND  MAGNETIC  FIELD 

The  following  paragraphs  discuss  in  turn  interaction  processes  be¬ 
tween  particles,  trajectories  of  charged  particles,  and  particle  drift. 

Interaction  Processes 

Consider  first  the  short-range  collisions,  or  binary  collisions, 
of  which  the  simplest  type  is  the  elastic  collision.  A  particle  of 
charge  e  and  ness  m  placed  at  the  origin  of  coordinates,  exerts  a 
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For  a  plasma;  distant  encounters  are  also  important,  which  extend 
beyond  the  limit  r  .  The  theory  of  Debye-HUckel  for  electrolytes  in¬ 
troduces  a  shielding  distance  Xp,  called  Debye  distance.  When  a  single 
particle  in  a  vacuum  exerts  a  potential  e/r,  the  same  particle  in  a 
plasma  does  not  extend  its  influence  to  such  a  long  distance,  on 
account  of  the  cloud  of  particles  surrounding  it  and  shielding  its 
electric  tield.  On  a  dimensional  argument,  the  balance  between  the 
kinetic  energy  and  the  potential  energy  is  written 

^7  v2  =  |  V  tp  (4.45) 

where  v  is  the  velocity  and  cp  is  the  potential.  ji  the  thermal  speed 
is  taken  as  representative  of  the  average  velocity,  and  the  potential 
e/X^  is  truncated  at  the  Debye  distance,  and  the  same  equation  of 
energy  balance  is  applied  for  n  particles  having  charges  n  e  in  a  Debye 
sphere,  one  finds,  after  omitting  7, 


.  -a  KT  ,  ne  e 
v  m  m  XD 


(4.46) 


yielding 

K,  =  r-  =  (4rrne2/KT)^  (4.47) 

u  a.d 

j. 

The  numerical  coefficient  (4n)2  comes  from  an  analytical  theory;  kp  is 
called  the  Debye  wave  number. 

Many  problems  of  plasmas  exhibit  features  which  do  not  depend  on 
collisions  of  the  scales  rQ  and  X^,  e.g.,  in  oscillations  of  wave¬ 
length  larger  than  X^  and  rQ.  Under  such  a  circumstance,  the  motion 
of  an  electron  is  characterized  by  a  self-consistent  electric  field  E 
and  is  determined  by  the  dynamical  equation 


(4.48) 
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and  the  Maxwell  equation 
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The  equation  possesses  simple  solutions  in  several  special  cases  as 
described  below.  When  B_  vanishes,  and  E  is  constant,  the  particle 
moves  with  a  constant  acceleration  eE/m.  When  E  vanishes,  tbe 
acceleration  e  v  x  B/mc  is  perpendicular  to  v  and  brings  no  change  to 
the  absolute  value  of  v: 


at  <*/>  - 0 


(4.53) 


entailing 


2  2  2 

v  s  v,,  +  v  =  constant 


Since  B  does  not  change  the  parallel  component  of  velocity,  it  will  leave 


v  =  constant 
jl 
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FIGURE 68.  Characteristics  of  Plasmas  (After  Tchen,  1968) 


191 


If  E  vanishes  and  B  is  constant  in  time  and  space,  the  motion  in 
the  plane  perpendicular  to  B  will  be  a  gyration  of  angular  frequency, 
called  cyclotron  frequency  (or  gyro  frequency) 


ibc  =  eB/mc 


(4.54) 


for  a  singly  charged  particle,  or 


u3c  =  Z  e  B/mc 


(4.55) 


for  a  multiply  charged  (Z  times)  particle.  The  gyration  perpendicular 
to  B  has  a  radius 


“c 


(4.56) 


The  component  v,  will  not  be  affected  by  B.  When  the  two  motions  are 
combined,  the  resultant  particle  path  will  be  a  helix  of  constant 
pitch  around  the  magnetic  line  of  force. 

Particle  Drifts  in  An  Electric  Field  and  Gravitational  Field 

If  E  and  B  are  constant  in  time  and  space,  one  finds  a  solution 
of  4.42,  constant  in  time  and  space. 


v 


E  x  B 


called  the  drift  velocity.  The  complete  solution 


(4.57a) 


v  =  V  +  v*  (4. 57b) 
~  ~d  ~ 

satisfying  4.42  consists  of  a  constant  part  and  a  variable  part  v ’ . 
The  latter  is  determined  by  the  following  equation  obtained  by  a  sub¬ 
stitution  of  4.57a  and  4.57b  into  4.42: 
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Iv  is  to  be  remarked  that  the  role  of  an  electric  force  eE  can  be 
replaced  by  any  other  force  (e.g.,  a  gravity  -  mg)  to  produce  a 
similar  feature  as  described  above.  In  view  of  the  change  of  charge  e 
into  m  the  gravitational  drift 

Vj  =  g^/u>c  (4.59) 


will  depend  on  the  sign  of  the  charge  of  the  particle. 


Particle  Drifts  in  An  Inhomogeneous  Magnetic  Field 


Consider  a  particle  moving  in  an  inhomogeneous  magnetic  field 
parallel  to  the  z-axis,  with  a  strength  varying  along  the  x-axis.  The 
particle  which  gyrates  in  the  xy  plane  will  change  its  radius  of  gyra¬ 
tion  over  the  orbit,  so  that  a  drift  will  result,  as  shown  in  Fig.  70. 
By  using  the  first  order  approximation  in  V^/v^,  Alfven  (1950) 
developed  the  formula 


(4.60) 


where  V  B  is  the  gradient  of  B  in  the 

i 


plane  perpendicular  to  B  and  vB. 


INHOMOGENEOUS  « 
MAGNETIC  FIELD  U 


GRADIENT  OF  B 


+  PARTICLE 


DRIFT 


-  PARTICLE 


DRIFT 


FIGURE  70. 


Drift  by  an  Inhomogeneous  Magnetic  Field 
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If  a  particle  moves  with  a  velocity  v^  along  a  curved  magnetic 
field  with  a  radius  of  curvature  R,  we  expect  a  similar  drift 


(4.61) 


/Ruc 


according  to  4.59  with  the  gravity  force  mg  replaced  by  the  centri- 
2 

fugal  force  mvy'R.  By  neglecting  the  current  for  a  strong  B,  we  have 


1 

R 


v  B 
x 

B 


The  two  drifts  4.60  and  4.61  from  the  inhomogeneous  magnetic 
field  amount  to 


+  %  v  2) 
x 


(4.62) 


They  are  in  the  same  direction  of  B  x  vB. 

The  orbital  motion  of  a  single  particle  in  an  external  electric 
and  magnetic  field  described  in  this  section  is  not  a  valid  description 
of  the  motion  of  an  assembly  of  particles.  Such  an  assembly  involves 
new  effects,  such  as  the  norm?.!  and  shear  stresses  due  to  the  inter¬ 
action  between  the  particles  of  like  kinds  or  unlike  kinds,  and  the 
effect  of  the  self-consistent  electric  field  induced  by  the  motion  of 
the  particles.  A  hydrodynamic  theory  to  include  these  effects  is 
elaborated  in  the  following  section. 


TRANSPORT  PROPERTIES  OF  PLASMAS  IN  THE  UPPER  ATMOSPHERE 

In  the  following  description  of  aggregates  of  charged  particles, 
the  topics  discussed  in  turn  are:  weakly  ionized  plasmas  such  as  the 
ionosphere  below  200  km,  arbitrary  degrees  of  ionization  such  as  the 
higher  ionosphere,  and  the  expression  of  some  features  in  geomagnetic 
polar  coordinates. 
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The  important  role  played  by  the  development  of  electric  and 
magnetic  field  on  the  motion  of  a  plasma  is  a  fundamental  feature  of 
plasma  dynamics,  as  distinguished  from  the  hydrodynamics  of  a  neutral 
gas.  The  diffusion  and  conductivity  properties  of  plasmas  are 
important  in  studying  the  origin  and  structure  of  plasma  inhomogeneities 
in  the  ionosphere,  and  their  scattering  of  radio  waves.  For  atmospheric 
conditions  of  dimension  inhomogeneities  longer  than  the  mean  free  path, 
it  is  necessary  to  consider  the  hydrodynamic  equations  of  motions  of  all 
constituents:  electrons,  ions,  and  neutral  molecules,  together  with  the 
Maxwell  equations  of  the  fields. 


The  equations  of  motion  of  the  three  components ,  denoted  by  indices 
a  and  b  are : 


d _ .  a_ 

dt  )t 


+  v 


~a 


.V 


(4.63b) 


where  the  number  density  n  and  the  velocity  v  satisfy  the  equation  of 

cl 

continuity 


-jt  +  v-(na  =  P  -  L 


(4.64) 


The  right-hand  *ide  represents  the  production  and  loss  processes  which 
will  be  subsequeucly  neglected.  Other  notations  are:  ma  is  the  mass, 
pa  the  pressure,  g  the  gravitational  acceleration,  eg  the  charge,  IE 
the  electric  field,  B  the  magnetic  field,  and  v  b  the  collision  fre¬ 
quency  between  the  particles  of  kinds  a  and  b- 
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In  the  low  altitude  of  the  atmosphere  (i.e.,  up  to  heights  of 
about  250  km),  the  magnetic  field  can  be  considered  as  constant  in 

2 

time  and  position,  as  the  electron  pressure  is  small  compared  to  B  /8tt. 
Consequently,  the  Maxwell  equation 

curl  E  =  0  (4.65) 

«v 

is  valid.  In  view  of  the  small  Debye  length,  the  condition  of  quasi¬ 
neutrality  holds 

n.  =  n  =  n  (4.66) 

-*•  0 

entailing 

div  ^  =  0  (4.67a) 

where 

j  =  en  (v.  -  v  )  (4.67b) 

is  the  electric  current.  The  divergence-free  condition  of  the  current 
is  a  consequence  of  the  equation  of  continuity  for  the  electrons  and 
ions. 

Further  the  pressures  satisfy  the  equation  of  state 

P-,  =  nKT  ,  K  =  Bolzmann  constant  (4.68) 

a  cl 

For  constant  Tg  and  v^  (n  =  neutral  molecules),  one  may  complete  the 
system  of  equations  for  the  determination  of  the  variables :  n  ,  v  , 
p ,  E,  j ■  It  will  be  convenient  to  express  the  gravity 

g  =  o>  n..  -g 


197 


198 


j 


By  further  introd.  Jing  the  notations 


v--=v  -  v  ,  a  =  e ,  i 

~a  ~a  ~n  ’ 


(4.73a) 


O.  eaB 

£  _  a  ==  - S - 

=  m  c  v 

a  van  a  an 


(4.73b) 


u  =  (—  E  -  t  )/Q 
~a  m  ~n  ~a  a 
3 


(4.73c) 


E  =  E  +  v_  x  B/c 
~n  ~  ~n  ~ 


(4.73d) 


the  asymptotic  momentum  equation  4.71  is  reduced  to 


v*  x  e„  -  ~  v*  +  u  =0 
~a  ~B  X  ~a  ~a 
a 


(4.74a) 


giving  the  components  of  v*  parallel  and  perpendicular  t 

"  ^  3 


V*7  =  K  U  ,• 

a"  a  all 


(4.74b) 


1  o 

- = - 7T  (K  u  xe^+K  u) 

2  +  K  2  a  ~a  ~B  a  ~a  x 


(4.74c) 


or  in  terms  of  E  and  F  : 

~n  ~a 


naeava//  "  aa//  En//  "  6 a//  e=  ^a/7 

a 


n  e  v*  =  a  n  E  -6  —  F 

a  a  ax  aP  aP  e  ~ax 

Q 


(4.75a) 


+  CTaH  (^n  x  ~B}x  ~  6aH  if  (Ea  x 
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a 


(4.75b) 


6 


aP 


aP 


3  6.  y 

m  0  aP 
3  3 


6 

an  an 


n  e2 
a  a 

ma°a 


K 


aH 


KaP  = 


1  +  K  2 
a 


KaH  =  ~ 


K  *' 

a 


1  +  tc  2 
a 


The  generalized  Ohms  law  (4.75a)  for  the  particles  a  can  be  re¬ 
written  as 


n  v*  =  —  5  .E  -  D  .7N 


a~a 


~n 


»a 


(4.76) 


where 


VN  =  n 


+ 


e 

H 


and  the  diffusion  tensor  D  is  shown  related  to  the  conductivity  tensor 

«$3 

cL  by  the  relation 


D 

«a 


KT 

a 

:a  efcn 

3 


(4.77) 


called  the  Einstein  relation. 

In  a  cartesian  coordinate  system  where  the  z-axis  (//)  is  parallel 
to  B  (as  in  Fig.  71),  and  have  the  form 


(Vig  = 


aP 


aH 


D 


aH 


-a 


aH 


aP 


0 

0 

<7/ 


-D 


aH 


aP 


D 


a!! 


The  columns  and  the  rows  are  numerated  by  the  indices  i  and  j  respectively. 
oap  is  called  the  Pedersen,  oaH  the  Hall,  and  aa»  the  direct  conductivity. 

cr u  is  the  conductivity  parallel  to  the  magnetic  field  and  is  that 
which  would  exist  for  all  directions  in  absence  of  the  magnetic 
field. 

Op  is  the  conductivity  parallel  to  that  component  of  electric 

field  that  is  normal  to  the  magnetic  field.  Op  is  often  called 
the  "Pedersen  conductivity." 

aH  is  the  conductivity  for  a  direction  perpendicular  to  both  elec¬ 
tric  and  magnetic  fields.  ou  is  often  called  the  "Hall  con- 

n 

ductivity . " 

From  the  above  transport  relations  for  each  constituent,  and  from 
their  difference  for  ions  and  electrons,  the  relation  for  the  current 
density  ^  nay  be  derived. 

2  =  g  *  In  "  e  5  '  VN  (4.80a) 


where 


(4.80b) 
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Fe,  AND  B  ARE  LOCALLY  IN  THE  PLANE  OF  THE  PAPER 
/.'  IS  PARALLEL  TO  B 

P  IS  NORMAL  TO  //  IN  THE  PLANE  OF  Fe  AND  B 
H  IS  NORMAL  TO  //,  P  AND  TO  Fe,  B 

FIGURE  71 .  The  // ,  P,  H  Coordinate  System 
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Figure  72  shows  the  tensors  versus  altitude.  According  to 
4.76  and  neglecting  vNa  and  v^ 


— —  a  E„ 
naea  a//  11 


e 

=  — K 
mana  a// 


7/ 


or 


va  B 

a// 


cE 


7/ 


(4.82) 


and  similarly  for  other  components.  Thus,  Fig.  72a  shows  the  formula 
4.82.  Furthermore,  according  to  4.81,  the  quantities 


A//S 


B 

nec  0 H 


- 


Ke>// 


etc 


are  shown  in  Fig.  72b. 

Arbitrary  Degree  of  Ionization 

The  above  transport  properties  have  been  derived  for  the  case  of 
a  weakly  ionized  plasma,  v/here  the  collision  between  the  charged 
particles  themselves  is  negligible  compared  to  the  collision  between  a 
charged  particle  and  a  neutral  particle.  This  is  the  case  a  ititudes 
lower  than  200  km.  At  higher  altitudes,  the  collision  between  the  ion 
and  electron  is  not  negligible  and  has  to  be  included  in  the  study  of 
transport  properties  and  plasma  motions.  The  plasma  has  then  a  high 
degree  of  ionization.  Thus  for  an  arbitrary  degree  of  ionization,  the 
transport  properties  are  found  to  be  Gurevich  and  Tsedilina,  67): 


V/  = 


2 

e  n 


m  (m  +  v  •  ) 
ev  ven  ei ' 


(4.83a) 


eR 


% 


—  Fv  •  +  V  (1  +  )1 

L  ei  en  1  in'J 


(4.83b) 
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ct  u=  -  n  ji  +  n>  v  ./m.v-  +  n?/v?  i 

eH  IP  A  ei-  e  ei  1  in  v  mj 


CTi //  =  -®J2 - 


IP  .  V .  V  .  +  V 
x  zn  ex  en 


0.  =  — £_2 —  (v2  +  v  v  .  +  fl2) 

ip-v-  A  en  en  ei  e 
x  in 

e2n  ni  2  2 

CTiH  =  -^-TJ  (ven  +  n’iveivir/n'e  +  V 
ivinH 


(4.83c) 


(4.83d) 


(4 . 83e) 


(4. 83f ) 
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(4.84c) 


(4.84d) 


(4.84e) 


(4. 84f ) 
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The  Einstein  relation  (Eq.  4.77)  is  invalid  for  a  plasma  with  an 
arbitrary  degree  of  ionization. 

The  following  conclusions  and  observations  can  be  drawn: 

1.  v.'ith  the  exception  of  D^,  all  transport  properties  trans¬ 
verse  to  the  magnetic  field  decrease  with  decreasing 
collision  frequency,  while  all  parallel  components  increase. 
The  large  parallel  conductivity  in  the  upper  atmosphere 
suggests  that  disturbances  produced  in  the  sporadic  E- layer 
may  influence  a  plasma  inhomogeneity  in  the  higher  altitude 
by  the  joining  lines  cf  force. 


2. 


3. 


The  current  density  nev*  formulated  above  in  Eq.  4.44  is 
useful  in  transforming  the  equation  of  continuity  into  an 
equation  of  diffusion,  by  a  substitution  for  nv^  Such  an 
equation  is  an  equation  of  diffusion  of  the  second  order,  if 
the  self-consistent  electrostatic  field  E  is  neglected.  The 
inclusion  of  the  effect  of  such  a  field  will  transform  into 
a  diffusion  equation  of  fourth  order,  which  is  responsible 
for  a  spreading  of  a  plasma  inhomogeneity  in  the  atmosphere 
along  with  the  diffusion  process. 


In  the  presence  of  an  external  electric  field  E^, 
velocities  of  the  electrons  and  ions  are 


the  drift 


V  =  v  - 


ae 
—  E 


en 


~eff ,e 


(4.85a) 


V.  =  V  +  —  E  or  . 
~xo  ~n  _  ~eff,x 
en  7 


(4.85b) 


depending  on  the  effective  electric  fields  J'or  electrons  and 
ions: 
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i  in  g 

E  =E+ivxB - — 

. 6  ^-O  'STi  _ 

5  c  e 


i  m±9 

E££.  =  E  -4  —  v  xB  +  -=— 
~«ff,i  ~o  ~n  ~ 

c  e 


(4.86a) 


(4.86b) 


4.  The  motion  of  plasmas  in  the  atmosphere  will  depend  on  the 
■■er/ime  of  the  plasmas,  i.e.,  whether  they  contain  fluctua¬ 
tions  or  not.  The  presence  of  the  micro- fluctuations  will 
change  the  transport  properties,  and  hence  the  motion  and  the 
structure  of  plasma  inhomogeneities  in  the  atmosphere.  A 
preliminary  attempt  of  the  investigation  of  such  nonlinear 
effects  has  been  made  by  Tchen  (1969). 

5.  Experiments  with  visible  artificial  ion  clouds  as  plasma 
inhomogeneities  have  been  performed  as  a  means  of  studying 
the  interaction  processes  between  the  interplanetary  plasma 
and  the  artificial  plasma  inhomogeneities  (Haerendel,  Lttst, 
and  Rieger,  1967).  Wie  experiments  help  in  the  understanding 
of  the  interplanetary  plasma  and  its  properties. 

Geomagnetic  Polar  Coordinates 

For  many  purposes  it  is  advantageous  to  use  a  system  of  magnetic 
polar  coordinates  described  in  the  following  paragraphs : 

In  the  F- region  and  above,  is  very  much  larger  than  Op  and 
which  become  vanishingly  small. 

In  the  D-  and  E-regions,  a  horizontal  current  system  nay  be 
imagined  to  be  flowing  in  a  relatively  thin  spherical  shell  with 
possible  conjugate  sources  and  sinks  in  the  Northern  and  Southern 
Hemispheres.  In  this  special  case,  except  at  highest  latitudes, 
vertical  current  density  may  be  ignored.  Consider  a  coordinate  system 
(x,  y,  z)  whose  axes  point  to  the  south,  east,  and  zenith,  respectively. 
In  this  new  coordinate  system,  the  tensor  $  of  (4.78)  becomes,  when  I 
is  the  magnetic  dip  angle,  positive  in  the  Northern  Hemisphere.  (Hines 
et  a!.,  1965). 
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o  '= 


.  2 


o  sin  I 
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°// 


(o//  -  Op)  sin  I  cos  I 


Ojj  sin  I 

(o/7 

-  °p 

°// 

o 

cos  I 

-I^J  cos  I 

0 

cos2  : 

.  2 


(4.87) 


Since  vertical  current  density  J  vanishes,  one  may  eliminate  E.  , 
for  example,  from  the  three  component  equation  (4.80a),  written  in 
the  x,y,z  coordinate  system,  leading  to  a  two  dimensional  relation¬ 
ship  between  horizontal  components  of  the  total  electric  field  and 

J  ,  J  of  the  current  density  in  the  form 
x  y 


J  = 

o  , 

r» 

+  a 

E 

(4.88a) 
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XX 

X 

xy 
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J  = 

:  o 
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+  a. 
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(4.88b) 
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CTXX  “ 

cxy  = 
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F  ap  <7/ 
r-l 

r  o2  vg 

sin  I 

(4.88c) 

(4.88d) 

CTyy  " 
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2_x 

+  a3  cos  I) 

(4-88e) 

r  = 

o 

cp  COS  I 

+  <7/ 

sin2I) 

(4.88f) 

CJ3  = 

2 

.  a!l 
°P  *  Op 

• 

(4.88g) 

o  ,  o ,  a  may  be  regarded  as  components  cf  a  two-dimensional  con- 
xx  xy  yy 

ductivity  tensor  relating  horizontal  current  density  to  horizontal 
components  of  the  total  electric  field. 


In  the  special  case  where  E  is  independent  of  height,  the  height 
integrated  conductivity  tensor,  with  components  ,faxx  ^2>  ?axy  dz> 

Jo  dz  relates  the  height  integrated  horizontal  component  of  the 
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current  density  fj  dz  to  the  horizontal  component  of  the  to 
field.  These  components  which  play  an  important  part  in  th 
the  dynamo  current  system  are  shown,  for  typical  mid-day  condition, 
Fig.  73. 


FIGURE  73.  Height  Integrated  Conductivity  Versus  Latitude 
(After  J.  A.  Fejer,  1953) 


2r 


*E 


Over  much  of  the  earth  (i.e.,  in  latitudes  higher  than  10  deg), 
the  height  integrated  conductivity  Ja  dz  by  which  meridional  currents 
result  from  east-west  variations  of  potential,  is  seen  to  predominate. 
At  the  magnetic  dip  equator,  however,  Jcrxy  dz  vanishes,  and  large 
values  are  assumed  by  the  height  integrated  conductivities  Jo  dz, 
by  which  meridional  fields  produce  meridional  currents,  and  foyy  dz, 
by  which  east-west  fields  produce  east-west  currents. 

Although  in  the  foregoing  treatment  the  gravitational  and  pressure 
terms  of  4.63a  are  neglected,  it  may  be  that  these  terms  can  assume 
importance  in  the  F-region  and  above.  Symmetry  about  the  magnetic 
equator  and  the  zero  value  of  the  current  density  parallel  to  B  in  the 
F-region  and  above  are  assumed. 

The  pressure  and  gravitational  terms  in  the  H  component  of  4.80a 
reduce  to 


-Er 


V  a:  v.  a: 
en  on 


(4.89) 


In  vector  notation  the  common  transverse  drift  velocity  of  4.89  is 
E  x  B  E  x  B 

— — ip  when  measured  with  respect  to  the  neutral  gas,  and  — — —  when 
B  B^ 

measured  with  respect  to  the  coordinate  system  in  which  E  is  measured. 

In  the  //  component  (parallel  to  B)  of  4.89,  and  the  corresponding 
equation  for  ions,  the  pressure  and  gravitational  terms  cannot  be 
neglected.  The  condition  which  follows  from  J  =  0  is 


V  =  V.  =  V. 
e/i  iff  // 


(4.90) 


leading  to  (Hines  et  al. ,  1965) 


m  w  v„ 
e  e  // 


=  F  =  -e  E 

e// 


7/  +  meg  sin  1  "  |  V>e 


(4.91a) 


mi  wi  v//  =  Fi/X  e  E//  +  m±g  sin  I  -  -  v//P. 


(4.91b) 
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Addition  of  4.91a  and  4.31b  and  substitution  of 


Pe  =  Pi  =  KTN 


(4.9^) 


when  K  is  Boltzmann’s  constant  and  T  is  the  temperature,  assumed  equal 
for  ions  and  electrons,  results  in  the  velocity  of  "ambipolar  diffu¬ 
sion"  along  tne  lines  of  force,  to  be 


sin  I 


+  mi)g  + 


X  9(TN)1 
N  "dz  J 


(4.93) 


where  z  is  the  height  and  (TN)  is  assumed  independent  of  latitude. 
The  upward  component  vz  of  velocity  v^  is  given  by 


vz  =  "  v//  sin  1 


(4.94) 


Moreover,  isothermal  conditions  are  assumed  and  we  use  the  dimension¬ 
less  height 

h  =  |  (4.95a) 

where  the  ion  scale  height  is 

„  _  XT  (4.95b) 

n>±  g 

and  if  mg  and  meu)e  are  neglected  in  comparison  with  m^  and  m^uu,  then 
4.93  and  4.94  become 


(4-96) 


Equation  4.96  indicates  the  vertical  component  of  the  velocity  of  ambi¬ 
polar  diffusion  due  to  gravity  and  pressure  along  the  magnetic  lines 
of  force. 
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p  ,  4  89  and  4.94  indicated  drift  velocity  with  respect 

Equation  •  move  with  neutral  particle  velo- 

the  neutral  medium  which  the  assumption  is  made  for  the 

city  ^  with  respect  to  t  e  =  Q  then  4.89,  4.94,  and 

F-region  that  neutral  particle  velocity  ^  , 

4.96  indicate  velocities  with  respect  to  earth. 
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